{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "cell_id": "00000-b62acb5e-bf38-4990-b0b2-de1a4ccf4944",
    "output_cleared": false
   },
   "source": [
    "# SPECT simulation \n",
    "\n",
    "Analysis after the simulation spect_main.py"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:09.858500Z",
     "start_time": "2024-01-23T19:56:09.436562Z"
    },
    "cell_id": "00001-ec633223-4fbf-496a-990e-5c5826265081",
    "execution_millis": 2,
    "execution_start": 1605009142403,
    "output_cleared": false,
    "source_hash": "206465ed"
   },
   "outputs": [],
   "source": [
    "# %matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from matplotlib import collections  as mc\n",
    "import numpy as np\n",
    "from ipywidgets import interact\n",
    "import os\n",
    "from pathlib import Path\n",
    "import uproot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:09.858824Z",
     "start_time": "2024-01-23T19:56:09.858154Z"
    },
    "cell_id": "00002-5be4f8b1-56ff-4026-8ab0-0934df8e2854",
    "execution_millis": 35,
    "execution_start": 1605009142412,
    "output_cleared": false,
    "source_hash": "bd411e75"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The Current Working Directory (CWD) is: \n",
      " /Users/dsarrut/src/gate/dqprm/2025/gate-exercices/6_spect\n"
     ]
    }
   ],
   "source": [
    "# The following command display the current working directory (where jupyter has been launched)\n",
    "cwd = os.getcwd()\n",
    "print('The Current Working Directory (CWD) is: \\n', cwd)\n",
    "f = uproot.open('output/spect.root')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:09.858943Z",
     "start_time": "2024-01-23T19:56:09.858247Z"
    },
    "cell_id": "00003-c0a724eb-32bf-4cef-94c4-9a9d5913e516",
    "execution_millis": 4,
    "execution_start": 1605009142450,
    "output_cleared": false,
    "source_hash": "365ca338"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "List of tree : \n",
      " ['Hits_spect_crystal;1', 'Singles_spect_crystal;1', 'spectrum;1', 'scatter1;1', 'peak113;1', 'scatter2;1', 'scatter3;1', 'peak208;1', 'scatter4;1']\n"
     ]
    }
   ],
   "source": [
    "print(\"List of tree : \\n\", f.keys())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:09.872113Z",
     "start_time": "2024-01-23T19:56:09.858295Z"
    },
    "cell_id": "00004-740bc229-33ef-4fb7-ae32-98d6b3c28afb",
    "execution_millis": 61,
    "execution_start": 1605009142457,
    "output_cleared": false,
    "source_hash": "a77506f"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of hits: 10978.0\n",
      "Hits keys: \n",
      "PostPosition_X  PostPosition_Y  PostPosition_Z  TotalEnergyDeposit  PreStepUniqueVolumeID  GlobalTime  LocalTime  StepLength  TrackLength  "
     ]
    }
   ],
   "source": [
    "hits = f['Hits_spect_crystal']\n",
    "print('Number of hits:', hits.num_entries)\n",
    "print('Hits keys: ')\n",
    "for i in hits.keys():\n",
    "    print(i, ' ', end='')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "cell_id": "00005-8b47ee00-7104-4efa-9801-0b2fe98e45d5",
    "output_cleared": false
   },
   "source": [
    "# Plot 1: position of hits, hits times"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:10.093886Z",
     "start_time": "2024-01-23T19:56:09.867143Z"
    },
    "cell_id": "00006-a1153054-743a-4484-b3d7-256edea767c7",
    "execution_millis": 3252,
    "execution_start": 1605009142519,
    "output_cleared": false,
    "source_hash": "7f8e5400"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'Time in ns')"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/sAAAGJCAYAAADVDSPZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACQLElEQVR4nOzdd1wT9/8H8FdYYQdBIFBR0VrFbVERJ1YqKrUOOhxV3K0FrdK6vnWPUu3Q1kWn2FZqta6qrRZx1YoLRVylaq074CIMZd/vD3+5GlkJJGTwej4e99DcfXJ5X8hd8r7PkgiCIICIiIiIiIiIzIaFoQMgIiIiIiIiIt1isk9ERERERERkZpjsExEREREREZkZJvtEREREREREZobJPhEREREREZGZYbJPREREREREZGaY7BMRERERERGZGSb7RERERERERGaGyT4RERERERGRmWGybyQkEgnmzp2rUdn69etjxIgReo2nqvr06YOxY8dW2+t16NABU6dO1fl+CwsLMXXqVPj4+MDCwgL9+/cvs2z9+vXx0ksvVbjP/fv3QyKRYP/+/boL9Cn//vsvJBIJYmNj9fYaKrGxsZBIJDhx4oTeX4vInBny2l4d16XKyM7OhoeHB9atW1ep52vz3fokfVzXRowYgfr164uP7927BwcHB/z66686ew2imojXzpJKu3bOnTsXEonEgFGVHkdBQQF8fHywatUqA0Zlvpjsl0L1Ja9abG1t8dxzzyEyMhJpaWnVEsPhw4cxd+5cZGRkVMvr6dKff/6J33//HdOmTRPX/fXXX5g6dSpat24NJycneHl5ITQ0tNQfUiNGjFB7/x0dHdGgQQO88sor2LRpE4qLi0s8Z9q0aVi5ciUUCoVOj+Xbb7/FRx99hFdeeQVr167F5MmTdbp/lbi4OCxbtkwv+9aVVatWVcvNAyJzUFGyGBQUhObNm5e7j/Pnz2Pu3Ln4999/dRaXqZ3Hn332GZycnDBo0CBDh6Jzbm5uGDNmDGbNmmXoUIiMBq+dumFK105ra2tERUVh0aJFyM3NNXQ4ZsfK0AEYs/nz58PX1xe5ubk4dOgQVq9ejV9//RVnz56Fvb29Tl/r0aNHsLL6789x+PBhzJs3DyNGjICLi4ta2dTUVFhYGO99mo8++gg9evTAs88+K677+uuv8c033yAsLAxvv/02lEolvvjiC3To0AG7du1CcHCw2j6kUim+/vprAI/fm6tXr2L79u145ZVXEBQUhG3btsHZ2Vks369fPzg7O2PVqlWYP3++zo5l7969eOaZZ7B06VKd7bNr16549OgRbGxsxHVxcXE4e/YsJk2apJPXqFevHh49egRra2ud7A94/EVXu3Zto29VQmSqnr62nz9/HvPmzUNQUJBajXBVlHUel3ZdMrSCggJ89tlnmDx5MiwtLQ0djl689dZb+Pzzz7F371688MILhg6HyCTx2qnOFK+dI0eOxPTp0xEXF4dRo0YZOhyzYrwZoxHo3bs33njjDYwZMwaxsbGYNGkSrly5gm3btun8tWxtbdWS/fJIpVKdJnG6lJ6ejp07d+K1115TWz948GBcv34dX3/9NcaNG4cpU6bg6NGjcHV1LbWJpZWVFd544w288cYbGDt2LBYuXIjTp08jOjoa+/fvL9FFwMLCAq+88gq+++47CIKg0+N5+mZLVVlYWMDW1lavN2xULVJM5SJf0+Tk5Bg6BDJChry2V8d1SVs7duzAnTt3SnyfmBM/Pz80b97cpGoMiYwNr53qTPHa6eLigp49e/JaqAfG88k0Aaq77leuXAHwuD/3ggUL0LBhQ0ilUtSvXx//+9//kJeXp/a8EydOICQkBLVr14adnR18fX1L3LV6sl/h3LlzMWXKFACAr6+v2Jxd1RyptL5J//zzD1599VW4urrC3t4eHTp0wM6dO9XKqPoVbdiwAYsWLUKdOnVga2uLHj164NKlS2plL168iLCwMMjlctja2qJOnToYNGgQlEplue/Rzp07UVhYWKKm3t/fH46Ojmrr3Nzc0KVLF1y4cKHcfT5p+vTp6NmzJzZu3Ii///5bbduLL76Iq1evIjk5ucL95OTk4N1334WPjw+kUikaN26Mjz/+WLxRoOrzvm/fPpw7d078G2jSJ+vQoUNo3749bG1t0aBBA3z33Xdq25/u3xUUFISdO3fi6tWr4us8eSd6+fLlaNasGezt7VGrVi20bdsWcXFx5cZQWp/9ESNGwNHRETdv3kT//v3h6OgId3d3vPfeeygqKip3f/Xr18e5c+dw4MABMcagoCC1Mnl5eYiKioK7uzscHBwwYMAA3Llzp8S+fvvtN3Tp0gUODg5wcnJCaGgozp07V+7rq2RkZGDy5MmoX78+pFIp6tSpg+HDh+Pu3bsAgPz8fMyePRv+/v6QyWRwcHBAly5dsG/fvlLfn48//hgrV65EgwYNYG9vj549e+L69esQBAELFixAnTp1YGdnh379+uH+/fuVOhbV+3758mX06dMHTk5OGDp0KADgjz/+wKuvvoq6detCKpXCx8cHkydPxqNHjzR6P8i8PHltj42NxauvvgoA6N69e4lrkCbfK6Xtv6zzuLR+p6rmsikpKejWrRvs7e3x7LPP4ueffwYAHDhwAAEBAbCzs0Pjxo2xZ8+eEq958+ZNjBo1Cp6enpBKpWjWrBm+/fZbjd6PrVu3on79+mjYsGGJbRs3bkTTpk1ha2uL5s2bY8uWLSX6xJfl1KlT6N27N5ydneHo6IgePXrgyJEjpZZ9+PAh3nzzTbi5ucHZ2RnDhw/HgwcP1Mps27YNoaGh8Pb2hlQqRcOGDbFgwYIKr6sqL774IrZv367TG9VENQmvnerKu3Y+TdNcBnj8m6dbt25wcnKCs7Mz2rVrp/Z7tKq/aV588UUcOnSo1N9bVHlsxq+Fy5cvA3icpALAmDFjsHbtWrzyyit49913cfToUURHR+PChQvYsmULgMc1wz179oS7uzumT58OFxcX/Pvvv9i8eXOZrzNw4ED8/fff+PHHH7F06VLUrl0bAODu7l5q+bS0NHTs2BEPHz7ExIkT4ebmhrVr1+Lll1/Gzz//jAEDBqiV//DDD2FhYYH33nsPSqUSS5YswdChQ3H06FEAjxOmkJAQ5OXlYcKECZDL5bh58yZ27NiBjIwMyGSyMmM/fPgw3NzcUK9ePY3eU4VCIR6fpoYNG4bff/8d8fHxeO6558T1/v7+AB6PGdCmTZsyny8IAl5++WXs27cPo0ePRuvWrbF7925MmTIFN2/exNKlS+Hu7o7vv/8eixYtQnZ2NqKjowE8roUpz6VLl/DKK69g9OjRCA8Px7fffosRI0bA398fzZo1K/U577//PpRKJW7cuCF2F1DdGPnqq68wceJEvPLKK3jnnXeQm5uLlJQUHD16FEOGDNH8Tft/RUVFCAkJQUBAAD7++GPs2bMHn3zyCRo2bIjx48eX+bxly5ZhwoQJcHR0xPvvvw8A8PT0VCszYcIE1KpVC3PmzMG///6LZcuWITIyEj/99JNY5vvvv0d4eDhCQkKwePFiPHz4EKtXr0bnzp1x6tSpcn+oZ2dnizeHRo0aheeffx53797FL7/8ghs3bqB27drIzMzE119/jcGDB2Ps2LHIysrCN998g5CQEBw7dgytW7dW2+e6deuQn5+PCRMm4P79+1iyZAlee+01vPDCC9i/fz+mTZuGS5cuYfny5XjvvffUvmi1OZbCwkKEhISgc+fO+Pjjj8VuQBs3bsTDhw8xfvx4uLm54dixY1i+fDlu3LiBjRs3VvTnJBOgVCrFm1FPKigoKPd5Xbt2xcSJE/H555/jf//7n3jt8fPzq9T3CqDZefy0Bw8e4KWXXsKgQYPw6quvYvXq1Rg0aBDWrVuHSZMm4a233sKQIUPEsU2uX78OJycnAI+/mzp06ACJRILIyEi4u7vjt99+w+jRo5GZmVlht6XDhw/j+eefL7F+586deP3119GiRQtER0fjwYMHGD16NJ555ply9wcA586dQ5cuXeDs7IypU6fC2toaX3zxBYKCgsQf4E+KjIyEi4sL5s6di9TUVKxevRpXr14Vf+ADj5MLR0dHREVFwdHREXv37sXs2bORmZmJjz76qMKY/P39sXTpUpw7d67CvshENQWvnbq/dpZGk1wGeHydGzVqFJo1a4YZM2bAxcUFp06dwq5du8Tfo1X9TePv7w9BEHD48GGNBrwmDQlUwpo1awQAwp49e4Q7d+4I169fF9avXy+4ubkJdnZ2wo0bN4Tk5GQBgDBmzBi157733nsCAGHv3r2CIAjCli1bBADC8ePHy31NAMKcOXPExx999JEAQLhy5UqJsvXq1RPCw8PFx5MmTRIACH/88Ye4LisrS/D19RXq168vFBUVCYIgCPv27RMACH5+fkJeXp5Y9rPPPhMACGfOnBEEQRBOnTolABA2btyo0fv1pM6dOwv+/v4alT148KAgkUiEWbNmqa0PDw8XHBwcynyeKr7JkyeX2GZjYyOMHz++3NfdunWrAEBYuHCh2vpXXnlFkEgkwqVLl8R13bp1E5o1a6bJ4Qj16tUTAAgHDx4U16WnpwtSqVR49913xXWqv8O+ffvEdaGhoUK9evVK7LNfv34av/6Trly5IgAQ1qxZI64LDw8XAAjz589XK9umTRuN/mbNmjUTunXrVmK96nwJDg4WiouLxfWTJ08WLC0thYyMDEEQHn8mXVxchLFjx6o9X6FQCDKZrMT6p82ePVsAIGzevLnENtXrFhYWqn22BUEQHjx4IHh6egqjRo0S16neH3d3dzE+QRCEGTNmCACEVq1aCQUFBeL6wYMHCzY2NkJubq7Wx6J636dPn14i7ocPH5ZYFx0dLUgkEuHq1avlvh9k3FTnRXnL0+f209f2jRs3lrhWCILm3yulKes8Lu261K1bNwGAEBcXJ67766+/BACChYWFcOTIEXH97t27S1xzRo8eLXh5eQl3795Ve61BgwYJMpms1M+/SkFBgSCRSNSunSotWrQQ6tSpI2RlZYnr9u/fLwAocR19+ru1f//+go2NjXD58mVx3a1btwQnJyeha9eu4jrV38/f31/Iz88X1y9ZskQAIGzbtk1cV9pxvPnmm4K9vb14zRCEx9eC0q7zhw8fFgAIP/30U+lvBlENwmun/q6dc+bMEZ5M/TTNZTIyMgQnJychICBAePTokVrZJ3/3afqb5uk4VG7duiUAEBYvXlzm8ZH22Iy/HMHBwXB3d4ePjw8GDRoER0dHbNmyBc8884w4VU5UVJTac959910AEJvQq/p779ixo8K7kZX166+/on379ujcubO4ztHREePGjcO///6L8+fPq5UfOXKk2kAiXbp0AfC4KwAAseZ+9+7dePjwoVax3Lt3D7Vq1aqwXHp6OoYMGQJfX1+tp8xT1XpnZWWV2FarVq1S7wQ/6ddff4WlpSUmTpyotv7dd9+FIAj47bfftIrnSU2bNhXfT+Bxa4zGjRuL7622XFxccOPGDRw/frzSMT3trbfeUnvcpUuXSsf3pHHjxqlNpdKlSxcUFRXh6tWrAID4+HhkZGRg8ODBuHv3rrhYWloiICCgRFP7p23atAmtWrUq0VIFgPi6lpaW4me7uLgY9+/fR2FhIdq2bYuTJ0+WeN6rr76q1lJFVav3xhtvqI2hERAQgPz8fNy8ebPSx1Jaywk7Ozvx/zk5Obh79y46duwIQRBw6tSpct8PMg0rV65EfHx8iaVly5aV3md1fK+oODo6qo3m3LhxY7i4uMDPz0+tFlz1f9W1RBAEbNq0CX379oUgCGrnSUhICJRKZannpMr9+/chCEKJ75Nbt27hzJkzGD58uFrXsG7duqFFixblHktRURF+//139O/fHw0aNBDXe3l5YciQITh06BAyMzPVnjNu3Di1vsDjx4+HlZWV2nR5T57HWVlZuHv3Lrp06YKHDx/ir7/+KjcmAOIxVvTdRVST8Nqp22tnaTTNZeLj45GVlYXp06fD1tZWreyTv/uq+puG10L9YDP+cqxcuRLPPfccrKys4OnpicaNG4sDcFy9ehUWFhZqI84DgFwuh4uLi5jgdOvWDWFhYZg3bx6WLl2KoKAg9O/fH0OGDIFUKtVJnFevXi3R9BD4r8n51atX1ZoG1q1bV62c6uRS9UP09fVFVFQUPv30U6xbtw5dunTByy+/jDfeeKPcJvwqQgX9DnNycvDSSy8hKysLhw4dKtGXvyLZ2dkAIDZ3evq1K5pD9OrVq/D29i7x/Cffr8p6+r0FHr+/T/fx1NS0adOwZ88etG/fHs8++yx69uyJIUOGoFOnTpXan62tbYnuIFWJ70kVfa4uXrwIAGWOOP3k7AqluXz5MsLCwiqMY+3atfjkk0/w119/qX2Z+/r6Vhiz6vPt4+NT6vrKHouVlRXq1KlToty1a9cwe/Zs/PLLLyX+BhWNj0GmoX379mjbtm2J9ZrcmCxLdXyvqNSpU6fENVUmk1V4jty5cwcZGRn48ssv8eWXX5a67/T09Apf/+nvE9X1+envXtW68n4E37lzBw8fPkTjxo1LbPPz80NxcTGuX7+u1uWqUaNGauUcHR3h5eWlNqXXuXPnMHPmTOzdu7fEzQJNzmPVMRrD/NdExoLXTt1eO0ujaS6j6sZcUTejqv6m4bVQP5jsl6OsC82TKvpASiQS/Pzzzzhy5Ai2b9+O3bt3Y9SoUfjkk09w5MgRrRNdXShrhPYnLwyffPIJRowYgW3btuH333/HxIkTER0djSNHjpSatKi4ubmVmzjm5+dj4MCBSElJwe7duyvVP/Hs2bMASv+xl5GRofUYALqkyXurDT8/P6SmpmLHjh3YtWsXNm3ahFWrVmH27NmYN2+ezuLThYqOvbi4GMDjvu5yubxEOU1noyjPDz/8gBEjRqB///6YMmUKPDw8YGlpiejoaPHLSpOYdX0sUqm0xEi9RUVFePHFF3H//n1MmzYNTZo0gYODA27evIkRI0aIr0H0tOr8XqnqOfLGG28gPDy81LLl1dC5urpCIpHo5EakPmVkZKBbt25wdnbG/Pnz0bBhQ9ja2uLkyZOYNm2aRuex6hgN+d1FVBPw2lk6XSTXuvhNw2uhfjDZr6R69eqhuLgYFy9eVBu0LS0tDRkZGSUGqOvQoQM6dOiARYsWIS4uDkOHDsX69esxZsyYUvevzYlXr149pKamllivaj6o6WB5T2vRogVatGiBmTNn4vDhw+jUqRNiYmKwcOHCMp/TpEkTbNq0qdRtxcXFGD58OBISErBhwwZ069atUnF9//33kEgkePHFF9XW37x5E/n5+RUOolevXj3s2bMHWVlZarX7VX2/qqK8v7eDgwNef/11vP766+LNkkWLFmHGjBklmlMZKkZNqEaF9fDwKDFbg6bPV93oKcvPP/+MBg0aYPPmzWrxzpkzR+vXqygWoPLHAgBnzpzB33//jbVr12L48OHi+vj4eJ3ESKavonNO2+8VTfapK+7u7nByckJRUVGlzhErKys0bNhQnP1GRXV9fnoGmbLWPR2Tvb19md+XFhYWJWrdLl68iO7du4uPs7Ozcfv2bfTp0wfA45G47927h82bN6Nr165iuafjLo+qbEXfXUSkGV47S147S6NpLqP6zXP27NlSK9oA3fym4bVQP9hnv5JUX/TLli1TW//pp58CAEJDQwE8vkv1dK2uakTw0qa1UHFwcADwuNZAk1iOHTuGxMREcV1OTg6+/PJL1K9fH02bNq1wH0/KzMxEYWGh2roWLVrAwsKi3JgBIDAwEA8ePCi1D/iECRPw008/YdWqVRg4cKBWMal8+OGH+P333/H666+XaF6ZlJQEAOjYsWO5++jTpw+KioqwYsUKtfVLly6FRCJB7969KxVbVTg4OJTaxOnevXtqj21sbNC0aVMIgqD3/mZPc3Bw0OjzWJaQkBA4Ozvjgw8+KDX20qbpe1JYWBhOnz6tNjqsiuocU90xf/KcO3r0qNq5oQtVPRag9FgFQcBnn32mu0DJpJX1PVDZ7xXVPqtyHmvK0tISYWFh2LRpU6k36TQ5RwIDA3HixAm1dd7e3mjevDm+++47sUsX8HgqqzNnzlQYU8+ePbFt2za1ZvhpaWmIi4tD586dS3TB+fLLL9XO8dWrV6OwsFD8nijtPM7Pz8eqVasqPD6VpKQkyGSyMmdsISLt8NpZ8tpZGk1zmZ49e8LJyQnR0dHIzc1VK1ve7y9tf9MkJSVBIpEgMDBQ4+dQxVizX0mtWrVCeHg4vvzyS7EZ37Fjx7B27Vr0799frAlYu3YtVq1ahQEDBqBhw4bIysrCV199BWdnZ/EkK41qGrn3338fgwYNgrW1Nfr27StewJ40ffp0/Pjjj+jduzcmTpwIV1dXrF27FleuXMGmTZtKNB+uyN69exEZGYlXX30Vzz33HAoLC/H999+LF6DyhIaGwsrKCnv27MG4cePE9cuWLcOqVasQGBgIe3t7/PDDD2rPGzBggNqxFRYWimVyc3Nx9epV/PLLL0hJSUH37t1L7ccUHx+PunXrljvtHgD07dsX3bt3x/vvv49///0XrVq1wu+//45t27Zh0qRJGs1Lqmv+/v746aefEBUVhXbt2sHR0RF9+/ZFz549IZfL0alTJ3h6euLChQtYsWIFQkNDSx2zQN8xrl69GgsXLsSzzz4LDw+PMvusl8bZ2RmrV6/GsGHD8Pzzz2PQoEFwd3fHtWvXsHPnTnTq1KnEDZgnTZkyBT///DNeffVVjBo1Cv7+/rh//z5++eUXxMTEoFWrVnjppZewefNmDBgwAKGhobhy5QpiYmLQtGlTtcSgqqp6LMDjVjANGzbEe++9h5s3b8LZ2RmbNm0y+mbLVH1at24NS0tLLF68GEqlElKpFC+88ALi4uIq9b0CVP081saHH36Iffv2ISAgAGPHjkXTpk1x//59nDx5Env27KlwLuV+/frh+++/x99//602zeoHH3yAfv36oVOnThg5ciQePHiAFStWoHnz5hWe5wsXLkR8fDw6d+6Mt99+G1ZWVvjiiy+Ql5eHJUuWlCifn5+PHj164LXXXkNqaipWrVqFzp074+WXXwbw+OZyrVq1EB4ejokTJ0IikeD777/XqutWfHw8+vbty36qRDrCa2fp186naZrLODs7Y+nSpRgzZgzatWuHIUOGoFatWjh9+jQePnyItWvX6uQ3TXx8PDp16iROcU46ou/h/k2RatqPiqbmKCgoEObNmyf4+voK1tbWgo+PjzBjxgy1qXZOnjwpDB48WKhbt64glUoFDw8P4aWXXhJOnDihti88NT2QIAjCggULhGeeeUawsLBQm4bv6SlGBEEQLl++LLzyyiuCi4uLYGtrK7Rv317YsWOHWhnV9CBPT6n39DRt//zzjzBq1CihYcOGgq2treDq6ip0795d2LNnTwXv3GMvv/yy0KNHD7V1qunHylqenGLw6bL29vZC/fr1hbCwMOHnn38WpxJ8UlFRkeDl5SXMnDlToxizsrKEyZMnC97e3oK1tbXQqFEj4aOPPlKbQkQQtJ96LzQ0tMT6bt26qU3XUto0LdnZ2cKQIUMEFxcXtemjvvjiC6Fr166Cm5ubIJVKhYYNGwpTpkwRlEplubGUNfVeaVMaljUFytMUCoUQGhoqODk5CQDEYyrrfCntOFXrQ0JCBJlMJtja2goNGzYURowYUeKcKM29e/eEyMhI4ZlnnhFsbGyEOnXqCOHh4eL0NMXFxcIHH3wg1KtXT5BKpUKbNm2EHTt2lJjySvX+fPTRR6XG/PQ5Ut4xVnQs5U0lef78eSE4OFhwdHQUateuLYwdO1Y4ffp0ib8dmZ6KvkdKu7aUdm3/6quvhAYNGgiWlpbi+aTp90ppyjqPy5o+qrTrX1nXOgBCRESE2rq0tDQhIiJC8PHxEaytrQW5XC706NFD+PLLLyuMNS8vT6hdu7awYMGCEtvWr18vNGnSRJBKpULz5s2FX375RQgLCxOaNGlSIqanv1tPnjwphISECI6OjoK9vb3QvXt34fDhw2plVH+/AwcOCOPGjRNq1aolODo6CkOHDhXu3bunVvbPP/8UOnToINjZ2Qne3t7C1KlTxem0nnw/S5t678KFCwL+f6pfIuK1U0Uf187Sfu9pksuo/PLLL0LHjh0FOzs7wdnZWWjfvr3w448/its1/U1TWhwZGRmCjY2N8PXXX1d4fKQdiSBUcuQwojL88ccfCAoKwl9//VWiqb2+bN26FUOGDMHly5fh5eVVLa9JRET6tWDBAqxZswYXL16scIDR1q1bw93d3aTGvZg0aRIOHjwoNl8lItIFba6dxmDZsmVYsmQJLl++rDaFH1Ud++yTznXp0gU9e/YstUmkvixevBiRkZFM9ImIzMjkyZORnZ2N9evXi+sKCgpKjCuzf/9+nD59GkFBQdUcYeXdu3cPX3/9NRYuXMhEn4h0qrRrp7EqKCjAp59+ipkzZzLR1wPW7BMREZHJ+PfffxEcHIw33ngD3t7e+OuvvxATEwOZTIazZ8+yvycREdH/4wB9REREZDJq1aoFf39/fP3117hz5w4cHBwQGhqKDz/8kIk+ERHRE1izT0RERERERGRmTKbPfnR0NNq1awcnJyd4eHigf//+SE1NVSuTm5uLiIgIuLm5wdHREWFhYUhLS1Mrc+3aNYSGhsLe3h4eHh6YMmVKib5/RERERERERKbMZJL9AwcOICIiAkeOHEF8fDwKCgrQs2dP5OTkiGUmT56M7du3Y+PGjThw4ABu3bqFgQMHituLiooQGhqK/Px8HD58GGvXrkVsbCxmz55tiEMiIiIioirQpDIoKCgIEolEbXnrrbfUymhSGbR//348//zzkEqlePbZZxEbG6vvwyMiqhKTbcZ/584deHh44MCBA+jatSuUSiXc3d0RFxeHV155BQDw119/wc/PD4mJiejQoQN+++03vPTSS7h16xY8PT0BADExMZg2bRru3LkDGxubCl+3uLgYt27dgpOTE0fPJaJKEQQBWVlZ8Pb2hoWFydxz1QivkURUFdpeH3v16oVBgwahXbt2KCwsxP/+9z+cPXsW58+fh4ODA4DHyf5zzz2H+fPni8+zt7eHs7MzgMeVQa1bt4ZcLsdHH32E27dvY/jw4Rg7diw++OADAMCVK1fQvHlzvPXWWxgzZgwSEhIwadIk7Ny5EyEhIRodG6+PRFQVlfr9KJioixcvCgCEM2fOCIIgCAkJCQIA4cGDB2rl6tatK3z66aeCIAjCrFmzhFatWqlt/+effwQAwsmTJ0t9ndzcXEGpVIrL+fPnBQBcuHDhUuXl+vXrOr82Gtr169cN/r5y4cLF9JfKXh/T09MFAMKBAwfEdd26dRPeeeedMp/z66+/ChYWFoJCoRDXrV69WnB2dhby8vIEQRCEqVOnCs2aNVN73uuvvy6EhIRoHBuvj1y4cNHFos310SRH4y8uLsakSZPQqVMnNG/eHACgUChgY2MDFxcXtbKenp5QKBRiGVWN/pPbVdtKEx0djXnz5pVYf/36dfGOMBGRNjIzM+Hj4wMnJydDh6JzqmPiNZKIKqOq10elUgkAcHV1VVu/bt06/PDDD5DL5ejbty9mzZoFe3t7AEBiYiJatGih9hsxJCQE48ePx7lz59CmTRskJiYiODhYbZ8hISGYNGlSmbHk5eUhLy9PfCz8f2NaXh+JqDIqc300yWQ/IiICZ8+exaFDh/T+WjNmzEBUVJT4WPUmOzs780JNRFVijs04VcfEayQRVUVlro+lVQYBwJAhQ1CvXj14e3sjJSUF06ZNQ2pqKjZv3gxAs8qgsspkZmbi0aNHsLOzKxFPWRVGvD4SUVVoc300uWQ/MjISO3bswMGDB1GnTh1xvVwuR35+PjIyMtRq99PS0iCXy8Uyx44dU9ufarR+VZmnSaVSSKVSHR8FEREREelSWZVB48aNE//fokULeHl5oUePHrh8+TIaNmyot3jKqjAiIqouJjMylCAIiIyMxJYtW7B37174+vqqbff394e1tTUSEhLEdampqbh27RoCAwMBAIGBgThz5gzS09PFMvHx8XB2dkbTpk2r50CIiIiISKdUlUH79u1TqwwqTUBAAADg0qVLAB5X+Dw9VfPTlUFllXF2di61Vh94XGGkqsVnbT4RGYLJJPsRERH44YcfEBcXBycnJygUCigUCjx69AgAIJPJMHr0aERFRWHfvn1ISkrCyJEjERgYiA4dOgAAevbsiaZNm2LYsGE4ffo0du/ejZkzZyIiIoK190REREQmpqLKoNIkJycDALy8vABoVhkUGBioVqGkKqOqUCIiMkYm04x/9erVAB5Pn/KkNWvWYMSIEQCApUuXwsLCAmFhYcjLy0NISAhWrVollrW0tMSOHTswfvx4BAYGwsHBAeHh4WpTsRARERGRaYiIiEBcXBy2bdsmVgYBjyuB7OzscPnyZcTFxaFPnz5wc3NDSkoKJk+ejK5du6Jly5YA1CuDlixZAoVCUaIy6K233sKKFSswdepUjBo1Cnv37sWGDRuwc+dOgx07EVFFJIJqaFDSSGZmJmQyGZRKJZtjEVGlmPN1xJyPjYj0T9trSFkDVakqg65fv4433ngDZ8+eRU5ODnx8fDBgwADMnDlTbf9Xr17F+PHjsX//frEy6MMPP4SV1X/1Yvv378fkyZNx/vx51KlTB7NmzRIrnPRxbERET6rMNcRkavaJiIiIiJ5UUZ2Vj48PDhw4UOF+6tWrh19//bXcMkFBQTh16pRW8RERGZLJ9NknIiIiIiIiIs0w2SciIiIiIiIyM0z2iYiIiIiIiMwMk30iMiopNzLw6e+pSLmRYehQTMrBgwfRt29feHt7QyKRYOvWrSXKXLhwAS+//DJkMhkcHBzQrl07XLt2Tdyem5uLiIgIuLm5wdHREWFhYSXmlSYiIiIi08Bkn4iMRsqNDLz2RSKW77uE175IZMKvhZycHLRq1QorV64sdfvly5fRuXNnNGnSBPv370dKSgpmzZoFW1tbsczkyZOxfft2bNy4EQcOHMCtW7cwcODA6joEIiIiItIhjsZPREZjz/k05BUWw8pCgrzCYuw5n4aWdVwMHZZJ6N27N3r37l3m9vfffx99+vTBkiVLxHUNGzYU/69UKvHNN98gLi4OL7zwAoDHU1f5+fnhyJEj6NChg/6CJyJ0X9sdALAvfJ+BIyGqGVTnHMDzjswXa/aJyGgEN/WE1MoChcUCpFYWCG7qaeiQzEJxcTF27tyJ5557DiEhIfDw8EBAQIBaU/+kpCQUFBQgODhYXNekSRPUrVsXiYmJZe47Ly8PmZmZagsRERERGR6TfSIyGi3ruGDDm4GY0P1ZbHgzkLX6OpKeno7s7Gx8+OGH6NWrF37//XcMGDAAAwcOFOefVigUsLGxgYuLi9pzPT09oVAoytx3dHQ0ZDKZuPj4+OjzUIiIiIhIQ2zGT0RGpWUdFyb5OlZcXAwA6NevHyZPngwAaN26NQ4fPoyYmBh069at0vueMWMGoqKixMeZmZlM+ImIiIiMAJN9IiIzV7t2bVhZWaFp06Zq6/38/HDo0CEAgFwuR35+PjIyMtRq99PS0iCXy8vct1QqhVQq1UvcRERERFR5bMZPRGTmbGxs0K5dO6Smpqqt//vvv1GvXj0AgL+/P6ytrZGQkCBuT01NxbVr1xAYGFit8RIRERFR1bFmn4jIDGRnZ+PSpUvi4ytXriA5ORmurq6oW7cupkyZgtdffx1du3ZF9+7dsWvXLmzfvh379+8HAMhkMowePRpRUVFwdXWFs7MzJkyYgMDAQI7ET0RERGSCmOwTEZmBEydOoHv3/6YRUvWjDw8PR2xsLAYMGICYmBhER0dj4sSJaNy4MTZt2oTOnTuLz1m6dCksLCwQFhaGvLw8hISEYNWqVdV+LERERERUdUz2iYjMQFBQEARBKLfMqFGjMGrUqDK329raYuXKlVi5cqWuwyMiIiKiasY++0RERERERERmhsk+ERERERERkZlhsk9ERERERERkZpjsExEREREREZkZJvtEREREREREZobJPhEREREREZGZYbJPREREREREZGaY7BMRERERERGZGSb7RERERERERGaGyT4RERERERGRmbEydABEREREVH26r+0u/n9f+D4DRkJERPrEmn0iIiIiIiIiM8Nkn4iIiIiIiMjMMNknIiIiIiIiMjNM9omIiIiIiIjMDJN9IiIiIiIiIjPDZJ+IiIiIiIjIzDDZJyIiIiIiIjIzTPaJiIiIiIiIzAyTfSIiIiIiIiIzw2SfiIiIiIiIyMww2SciIiIiIiIyM0z2iYiIiIioxuq+tju6r+1u6DCIdM6kkv2DBw+ib9++8Pb2hkQiwdatW9W2jxgxAhKJRG3p1auXWpn79+9j6NChcHZ2houLC0aPHo3s7OxqPAoiIiIiIiIi/TKpZD8nJwetWrXCypUryyzTq1cv3L59W1x+/PFHte1Dhw7FuXPnEB8fjx07duDgwYMYN26cvkMnIiIiIiIiqjZWhg5AG71790bv3r3LLSOVSiGXy0vdduHCBezatQvHjx9H27ZtAQDLly9Hnz598PHHH8Pb21vnMRMRERERERFVN5Oq2dfE/v374eHhgcaNG2P8+PG4d++euC0xMREuLi5iog8AwcHBsLCwwNGjR0vdX15eHjIzM9UWIiJjU1E3pye99dZbkEgkWLZsmdp6dnMiIiIiMh9mlez36tUL3333HRISErB48WIcOHAAvXv3RlFREQBAoVDAw8ND7TlWVlZwdXWFQqEodZ/R0dGQyWTi4uPjo/fjICLSlibdnABgy5YtOHLkSKktmdjNiYiIiMh8mFQz/ooMGjRI/H+LFi3QsmVLNGzYEPv370ePHj0qtc8ZM2YgKipKfJyZmcmEn4iMjibdnG7evIkJEyZg9+7dCA0NVdvGbk5ERERE5sWsavaf1qBBA9SuXRuXLl0CAMjlcqSnp6uVKSwsxP3798vs5y+VSuHs7Ky2EBGZmuLiYgwbNgxTpkxBs2bNSmyvTDcngF2diIiIiIyVWSf7N27cwL179+Dl5QUACAwMREZGBpKSksQye/fuRXFxMQICAgwVJhGR3i1evBhWVlaYOHFiqdsr080JYFcnIiIiImNlUsl+dnY2kpOTkZycDAC4cuUKkpOTce3aNWRnZ2PKlCk4cuQI/v33XyQkJKBfv3549tlnERISAgDw8/NDr169MHbsWBw7dgx//vknIiMjMWjQIDZRJSKzlZSUhM8++wyxsbGQSCQ63feMGTOgVCrF5fr16zrdPxERERFVjkkl+ydOnECbNm3Qpk0bAEBUVBTatGmD2bNnw9LSEikpKXj55Zfx3HPPYfTo0fD398cff/wBqVQq7mPdunVo0qQJevTogT59+qBz58748ssvDXVIRER698cffyA9PR1169aFlZUVrKyscPXqVbz77ruoX78+gMp1cwLY1YmIiIjIWJnUAH1BQUEQBKHM7bt3765wH66uroiLi9NlWERERm3YsGEIDg5WWxcSEoJhw4Zh5MiRANS7Ofn7+wNgNyciIiIiU2ZSyT4REZUuOztbHIwU+K+bk6urK+rWrQs3Nze18tbW1pDL5WjcuDEA9W5OMTExKCgoYDcnIiIiIhNmUs34iYiodOV1c9IUuzkRERERmQ8m+0REZkDVzenpJTY2ttTy//77LyZNmqS2TtXNKSsrC0qlEt9++y0cHR31HzwRUSVFR0ejXbt2cHJygoeHB/r374/U1FS1Mrm5uYiIiICbmxscHR0RFhaGtLQ0tTLXrl1DaGgo7O3t4eHhgSlTpqCwsFCtzP79+/H8889DKpXi2WefLfP6SkRkLJjsExGRWeq+tru4EJF5OnDgACIiInDkyBHEx8ejoKAAPXv2RE5Ojlhm8uTJ2L59OzZu3IgDBw7g1q1bGDhwoLi9qKgIoaGhyM/Px+HDh7F27VrExsaqtYy6cuUKQkND0b17dyQnJ2PSpEkYM2aMRuNFEREZCvvsExEREZFJ2rVrl9rj2NhYeHh4ICkpCV27doVSqcQ333yDuLg4vPDCCwCANWvWwM/PD0eOHEGHDh3w+++/4/z589izZw88PT3RunVrLFiwANOmTcPcuXNhY2ODmJgY+Pr64pNPPgHweJyTQ4cOYenSpeIUz0RExoY1+0RERERkFpRKJYDH3ZIAICkpCQUFBWozkjRp0gR169ZFYmIiACAxMREtWrSAp6enWCYkJASZmZk4d+6cWKa0WU1U+yhNXl4eMjMz1RYiourEZJ+IiIiITF5xcTEmTZqETp06oXnz5gAAhUIBGxsbuLi4qJX19PSEQqEQyzyZ6Ku2q7aVVyYzMxOPHj0qNZ7o6GjIZDJx8fHxqfIxEhFpg8k+EREREZm8iIgInD17FuvXrzd0KACAGTNmQKlUisv169cNHRIR1TDss09ERERUjicHedwXvs+AkWhOFbOpxFtVkZGR2LFjBw4ePIg6deqI6+VyOfLz85GRkaFWu5+Wlga5XC6WOXbsmNr+VKP1P1nm6RH809LS4OzsDDs7u1JjkkqlkEqlVT42IqLKYs0+EREREZkkQRAQGRmJLVu2YO/evfD19VXb7u/vD2trayQkJIjrUlNTce3aNQQGBgIAAgMDcebMGaSnp4tl4uPj4ezsjKZNm4plntyHqoxqH0RExog1+0RERERkkiIiIhAXF4dt27bByclJ7GMvk8lgZ2cHmUyG0aNHIyoqCq6urnB2dsaECRMQGBiIDh06AAB69uyJpk2bYtiwYViyZAkUCgVmzpyJiIgIsWb+rbfewooVKzB16lSMGjUKe/fuxYYNG7Bz506DHTsRUUVYs09EREREJmn16tVQKpUICgqCl5eXuPz0009imaVLl+Kll15CWFgYunbtCrlcjs2bN4vbLS0tsWPHDlhaWiIwMBBvvPEGhg8fjvnz54tlfH19sXPnTsTHx6NVq1b45JNP8PXXX3PaPSIyaqzZJyIiompjiv3fyXgJglBhGVtbW6xcuRIrV64ss0y9evXw66+/lrufoKAgnDp1SusYiYgMhTX7RERERERERGaGyT4REdH/6762u1rNMxEREZGpYrJPREREREREZGaY7BMREZHJYSsMIiKi8jHZJyIiIiIiIjIzHI2fqoVCmYtzt5Ro5i2DXGZr6HCIiIiIqAZiiyCqSZjsk94plLkY891xKJS5kMts8fXwdkz4iYhIbzi9n+Go3nu+70REhsdm/KR3524poVDmwsXORqzhJyIiIiIiIv1hzT7pnarpvqpmv5m3zNAhERFVK0PVNLOWlYiIqOZisk96p2q6zz77RERERERE1YPJPlULucyWST4RkREpbZAqtgAwP7puVcLxEIiITAf77BMRERERERGZGdbsU4U4bR4RUcXYP550oTI15/zsERFRaZjsU7k4bR4REREREZHpYTN+KhenzSMiourUfW33UscTICIiIu0w2adyqZruZzzK57R5REbs4MGD6Nu3L7y9vSGRSLB161ZxW0FBAaZNm4YWLVrAwcEB3t7eGD58OG7duqW2j/v372Po0KFwdnaGi4sLRo8ejezs7Go7BlWSx0SPTAk/t0REZKyY7FO5VE33F4e1ZBN+IiOWk5ODVq1aYeXKlSW2PXz4ECdPnsSsWbNw8uRJbN68GampqXj55ZfVyg0dOhTnzp1DfHw8duzYgYMHD2LcuHHVdQhEREREpEPss08V4rR5RMavd+/e6N27d6nbZDIZ4uPj1datWLEC7du3x7Vr11C3bl1cuHABu3btwvHjx9G2bVsAwPLly9GnTx98/PHH8Pb21vsx1GQcYE1/zOG95XR3RERUGazZJyKqgZRKJSQSCVxcXAAAiYmJcHFxERN9AAgODoaFhQWOHj1a5n7y8vKQmZmpthARERGR4THZJyKqYXJzczFt2jQMHjwYzs7OAACFQgEPDw+1clZWVnB1dYVCoShzX9HR0ZDJZOLi4+Oj19ipbOw3TkRERE9isk9EVIMUFBTgtddegyAIWL16dZX3N2PGDCiVSnG5fv26DqIkIiIioqpin30iohpClehfvXoVe/fuFWv1AUAulyM9PV2tfGFhIe7fvw+5XF7mPqVSKaRSqd5iJsNjf/GaxxzGOSAiItbsExHVCKpE/+LFi9izZw/c3NzUtgcGBiIjIwNJSUniur1796K4uBgBAQHVHS4RERERVRFr9omIzEB2djYuXbokPr5y5QqSk5Ph6uoKLy8vvPLKKzh58iR27NiBoqIisR++q6srbGxs4Ofnh169emHs2LGIiYlBQUEBIiMjMWjQIKMZib+02kZjrIFkTTiZIo73QERkfpjsExGZgRMnTqB79/9+rEdFRQEAwsPDMXfuXPzyyy8AgNatW6s9b9++fQgKCgIArFu3DpGRkejRowcsLCwQFhaGzz//vFriJyIiIiLdMqlk/+DBg/joo4+QlJSE27dvY8uWLejfv7+4XRAEzJkzB1999RUyMjLQqVMnrF69Go0aNRLL3L9/HxMmTMD27dvFH7OfffYZHB0dDXBERES6ERQUBEEQytxe3jYVV1dXxMXF6TIsk1UdtfOsSTVObJlBRETmwqT67Ofk5KBVq1ZYuXJlqduXLFmCzz//HDExMTh69CgcHBwQEhKC3NxcsczQoUNx7tw5xMfHY8eOHTh48CDGjRtXXYdAJkqhzEXChTQolLkVFyYiIo1wukDzpfrb8u9LRGQ4JlWz37t3b/Tu3bvUbYIgYNmyZZg5cyb69esHAPjuu+/g6emJrVu3YtCgQbhw4QJ27dqF48ePo23btgCA5cuXo0+fPvj444+Npl8qqVMoc3HulhLNvGWQy2wN8vpjvjsOhTIXcpktvh7eziBxEBERERERaUrrZP/evXuYPXs29u3bh/T0dBQXF6ttv3//vs6C08aVK1egUCgQHBwsrpPJZAgICEBiYiIGDRqExMREuLi4iIk+AAQHB8PCwgJHjx7FgAEDSuw3Ly8PeXl54uPMzEz9HgipMYZE+9wtJRTKXLjY2Yg3Hkwh2Tf0TRIiTWlb81dTagrLG3yQTc1rhpryWSciIv3QOtkfNmwYLl26hNGjR8PT0xMSiUQfcWlNNbK0p6en2npPT09xm0KhgIeHh9p2KysruLq6imWeFh0djXnz5ukhYtKEMSTaqmRZdcOhmbesWl+/MozhJgkRERERERmO1sn+H3/8gUOHDqFVq1b6iMfozJgxQxzVGnhcs+/j42PAiGoWY0i0VcmyKdWSG8NNEiKqGYxx+kOqHvzbExEZN62T/SZNmuDRo0f6iKVK5HI5ACAtLQ1eXl7i+rS0NHGqKblcjvT0dLXnFRYW4v79++LznyaVSiGVSvUTNFXIWBJtuczWpJJlY7hJQkREREREhqP1aPyrVq3C+++/jwMHDuDevXvIzMxUWwzF19cXcrkcCQkJ4rrMzEwcPXoUgYGBAIDAwEBkZGQgKSlJLLN3714UFxcjICCg2mMmzchltujh52lSybahqW6SLA5rySb8RFXEEcWJiIjIFGlds+/i4oLMzEy88MILausFQYBEIkFRUZHOgntadnY2Ll26JD6+cuUKkpOT4erqirp162LSpElYuHAhGjVqBF9fX8yaNQve3t7o378/AMDPzw+9evXC2LFjERMTg4KCAkRGRmLQoEEciZ/Mjqm1RiAiIiIiIt3ROtkfOnQorK2tERcXV+0D9J04cQLdu/9Xu6LqSx8eHo7Y2FhMnToVOTk5GDduHDIyMtC5c2fs2rULtrb/JTzr1q1DZGQkevToAQsLC4SFheHzzz+vtmMgIiLTZ2p9lTUdvZ+j/BMREZkPrZP9s2fP4tSpU2jcuLE+4ilXUFAQBEEoc7tEIsH8+fMxf/78Msu4uroiLi5OH+ERaYzT4hERERERkT5p3We/bdu2uH79uj5iIaoRVNPiTduUIk6PV5MolLlIuJBW446biIiIiKg6aV2zP2HCBLzzzjuYMmUKWrRoAWtra7XtLVu21FlwROaoJk+Lp7rRoZolgIMHkrnjwH41E7tDEBGRMdA62X/99dcBAKNGjRLXSSSSahmgj8gc1ORp8WryjQ4iIiIiouqkdbJ/5coVfcRBVC2Moa+8qkbb0HEYQk2+0UFkatgqgdhCgYjItGmd7NerV08fcRDpnTE1Ia+p0+LV5BsdRERERETVSetkHwBu3bqFQ4cOIT09HcXFxWrbJk6cqJPAiHSNTciNQ0290UHVhzXSRERERJVI9mNjY/Hmm2/CxsYGbm5ukEgk4jaJRMJkn4wWm5ATEREREVFNoXWyP2vWLMyePRszZsyAhYXWM/cRGQybkBORqfdBNpZWC8YSB2D6f1NdMqa/CxERGZ7Wyf7Dhw8xaNAgJvpkktiEnIiIiIiIagKtk/3Ro0dj48aNmD59uj7iISIiIjJJ2tassyaeiIj0SetkPzo6Gi+99BJ27dqFFi1awNraWm37p59+qrPgiIiIiIiIiEh7lUr2d+/ejcaNGwNAiQH6iIiIiIiIiMiwtE72P/nkE3z77bcYMWKEHsIhIiIyTmxybfyq42+keo2aNhhgee8tzw0iIuOk9Sh7UqkUnTp10kcsRERERERERKQDWif777zzDpYvX66PWIiIiCql+9ru4lLediJTx8+yuoMHD6Jv377w9vaGRCLB1q1b1baPGDECEolEbenVq5damfv372Po0KFwdnaGi4sLRo8ejezsbLUyKSkp6NKlC2xtbeHj44MlS5bo+9CIiKpM62b8x44dw969e7Fjxw40a9asxAB9mzdv1llwRERERERlycnJQatWrTBq1CgMHDiw1DK9evXCmjVrxMdSqVRt+9ChQ3H79m3Ex8ejoKAAI0eOxLhx4xAXFwcAyMzMRM+ePREcHIyYmBicOXMGo0aNgouLC8aNG6e/gyMiqiKta/ZdXFwwcOBAdOvWDbVr14ZMJlNbiIio+lVUuyUIAmbPng0vLy/Y2dkhODgYFy9eVCujSe0WEWnHmGrijSkWXenduzcWLlyIAQMGlFlGKpVCLpeLS61atcRtFy5cwK5du/D1118jICAAnTt3xvLly7F+/XrcunULALBu3Trk5+fj22+/RbNmzTBo0CBMnDiRM1ARkdHTumb/yTujRERkHCqq3VqyZAk+//xzrF27Fr6+vpg1axZCQkJw/vx52NraAqi4douIyBTt378fHh4eqFWrFl544QUsXLgQbm5uAIDExES4uLigbdu2Yvng4GBYWFjg6NGjGDBgABITE9G1a1fY2NiIZUJCQrB48WI8ePBA7ebBk/Ly8pCXlyc+zszM1NMREhGVTutkn4iIjE/v3r3Ru3fvUrcJgoBly5Zh5syZ6NevHwDgu+++g6enJ7Zu3YpBgwaJtVvHjx8Xf/QuX74cffr0wccffwxvb+9qOxbAvEc8N3TNqr5ev7T9muPfz1gY+nNkKnr16oWBAwfC19cXly9fxv/+9z/07t0biYmJsLS0hEKhgIeHh9pzrKys4OrqCoVCAQBQKBTw9fVVK+Pp6SluKyvZj46Oxrx58/RwVEREmtG6GT8REZmWK1euQKFQIDg4WFwnk8kQEBCAxMREABXXbpUlLy8PmZmZagsRkbEYNGgQXn75ZbRo0QL9+/fHjh07cPz4cezfv1/vrz1jxgwolUpxuX79ut5fk4joSazZJyIyc6raKVVNlIqnp6dazVVFtVulMYeaK1OvITX1+ImqU4MGDVC7dm1cunQJPXr0gFwuR3p6ulqZwsJC3L9/H3K5HAAgl8uRlpamVkb1WFWmNFKptMRggERE1YnJPhERVdqMGTMQFRUlPs7MzISPj48BIzJtTNyrh6bdRPj3MD83btzAvXv34OXlBQAIDAxERkYGkpKS4O/vDwDYu3cviouLERAQIJZ5//33UVBQIM5CFR8fj8aNG5fZhJ+IyBiwGT8RkZlT1TyVVjP1ZM1VRbVbpZFKpXB2dlZbiIiqS3Z2NpKTk5GcnAzgcbel5ORkXLt2DdnZ2ZgyZQqOHDmCf//9FwkJCejXrx+effZZhISEAAD8/PzQq1cvjB07FseOHcOff/6JyMhIDBo0SByrZMiQIbCxscHo0aNx7tw5/PTTT/jss8/UbnQSERmjStXsJyQkICEhAenp6SguLlbb9u233+okMCIi0g1fX1/I5XIkJCSgdevWAB7XwB89ehTjx48HoFntFlFVsJac9OHEiRPo3v2/z5YqAQ8PD8fq1auRkpKCtWvXIiMjA97e3ujZsycWLFig1rx+3bp1iIyMRI8ePWBhYYGwsDB8/vnn4naZTIbff/8dERER8Pf3R+3atTF79myMGzeu+g6UiKgStE72582bh/nz56Nt27bw8vKCRCLRR1xERKSF7OxsXLp0SXysqt1ydXVF3bp1MWnSJCxcuBCNGjUSp97z9vZG//79AajXbsXExKCgoKBE7RYRkbEJCgqCIAhlbt+9e3eF+3B1da1witGWLVvijz/+0Do+IiJD0jrZj4mJQWxsLIYNG6aPeIiIqBLKq92KjY3F1KlTkZOTg3HjxiEjIwOdO3fGrl27YGtrKz6nototIkPR11SMT7Y24DSBRERkbrRO9vPz89GxY0d9xEJktBTKXJy7pUQzbxnkMtuKn0BUzSqq3ZJIJJg/fz7mz59fZhlNareIiIiIyDRoneyPGTMGcXFxmDVrlj7iITI6CmUuxnx3HAplLuQyW3w9vB0TfiIiHWJ/fiIiIt3TOtnPzc3Fl19+iT179qBly5biFCQqn376qc6CIzIG524poVDmwsXORqzhZ7JPRERERETGTOtkPyUlRRzN+ezZs2rbOFgfmSNV031VzX4zb5mhQyIiIgNjawQiIjJ2Wif7+/ZxABsyDtXVj17VdL86XotjAxDR05hU6gbfRyIiqmm0TvaJjIE++tGXl2jLZbZ6T745NgARERGR4XCGDjI3GiX7AwcORGxsLJydnTFw4MByy27evFkngRGV59wtJW5l5EJqZYFbGVXvR28MiTbHBiAiIiIiIl3RKNmXyWRif3yZjP2VyfDcnaR4mF+IBw+LIbWygLuTtEr7M4ZEm2MDEBERERGRrmiU7K9Zs6bU/xPpW1lN6+9k5cHO2hK17GyQW1iEO1l5VXodY0i0q3NsACIyX+ybTkRERAD77JMRK69pfTNvGZ6pZQeFMhfP1LKrcnJuLIl2dYwNQERERERE5o/JPhmt8prW6yM5Z6JNRERERETmgsk+Ga2KmtbrIjnnVHdERERERGSOLAwdgC7NnTsXEolEbWnSpIm4PTc3FxEREXBzc4OjoyPCwsKQlpZmwIipPKra+8VhLXU2Or5CmYuEC2lQKHPFbgLTNqWI3QWIiIiIiIjMgdnV7Ddr1gx79uwRH1tZ/XeIkydPxs6dO7Fx40bIZDJERkZi4MCB+PPPPw0RKmlAl03rnx4DYESgr8Yj8LMFABERERERmZJKJfsJCQlISEhAeno6iouL1bZ9++23OgmssqysrCCXy0usVyqV+OabbxAXF4cXXngBwOOZBfz8/HDkyBF06NChukOlavb0GACAoNEI/Ck3MjAh7hSUuQWQ2Vpj+ZA2aFnHpVpjr068sUE1BUetJ0PhZ4+IiKqD1sn+vHnzMH/+fLRt2xZeXl6QSCT6iKvSLl68CG9vb9ja2iIwMBDR0dGoW7cukpKSUFBQgODgYLFskyZNULduXSQmJpaZ7Ofl5SEv779p3TIzM/V+DKQfT48B0LmROzo3ci83sVUocxEZdxLX7z+CAED5qACRcSex4c2OZpkIlzcDAhERERERmQ6tk/2YmBjExsZi2LBh+oinSgICAhAbG4vGjRvj9u3bmDdvHrp06YKzZ89CoVDAxsYGLi4uas/x9PSEQqEoc5/R0dGYN2+eniM3X9VVS6zJ65Q1gn95cZ27pURWbiEsLSQoLBZgKZEgK7ew3Cb/pqy8GRCIiIiIiMh0aJ3s5+fno2PHjvqIpcp69+4t/r9ly5YICAhAvXr1sGHDBtjZ2VVqnzNmzEBUVJT4ODMzEz4+PlWOtSaorlrilBsZiIw7iazcQjxTy67c19F2DIBm3jI8U8sO1+89xMOCIjjYWOGZWnZlNvk3dRXNgEBERERERKZB62R/zJgxiIuLw6xZs/QRj065uLjgueeew6VLl/Diiy8iPz8fGRkZarX7aWlppfbxV5FKpZBKpdUQrfmpjlpihTIXE+JO4fqDR7CxtMDNB490+jpPtgZwd5LiTlaeWfdlL6v1AxERERERmRatk/3c3Fx8+eWX2LNnD1q2bAlra2u17Z9++qnOgquq7OxsXL58GcOGDYO/vz+sra2RkJCAsLAwAEBqaiquXbuGwMBAA0dqnqqjlvjcLSWy8gphbSlBflExPG2lpb5OVboT6HJGAFNQ046XiIiIiMgcaZ3sp6SkoHXr1gCAs2fPqm0z9GB97733Hvr27Yt69erh1q1bmDNnDiwtLTF48GDIZDKMHj0aUVFRcHV1hbOzMyZMmIDAwECOxK8n1VFL3MxbBm8XW9zKAJykVlg+pE2J1+Ggc9rhaPxERERERKZP62R/3759+ohDJ27cuIHBgwfj3r17cHd3R+fOnXHkyBG4u7sDAJYuXQoLCwuEhYUhLy8PISEhWLVqlYGjNm/6riXW5IYCB53THG+MkLHhFGVkLvhZJiKi6qZ1sv+kGzduAADq1Kmjk2Cqav369eVut7W1xcqVK7Fy5cpqioiqQ0U3FDjonOZ4Y4SIiIiIyDxYaPuE4uJizJ8/HzKZDPXq1UO9evXg4uKCBQsWoLi4WB8xElWJqoZ6cVhL1lRXQHVjJONRPm+MEBERERGZMK1r9t9//3188803+PDDD9GpUycAwKFDhzB37lzk5uZi0aJFOg+SqKpq+qBzmvbD52j8RERERETmQetkf+3atfj666/x8ssvi+tatmyJZ555Bm+//TaTfSIjo20//Jp+Y4SIdIf91I0f/0ZEROZL62b89+/fR5MmTUqsb9KkCe7fv6+ToIhId0rrh081T1FREWbNmgVfX1/Y2dmhYcOGWLBgAQRBEMsIgoDZs2fDy8sLdnZ2CA4OxsWLFw0YtfHrvra7uNQENelYiYiITJ3WyX6rVq2wYsWKEutXrFiBVq1a6SQoItId9sMnAFi8eDFWr16NFStW4MKFC1i8eDGWLFmC5cuXi2WWLFmCzz//HDExMTh69CgcHBwQEhKC3NxcA0ZORFRz1bQbikSkW1o341+yZAlCQ0OxZ88eBAYGAgASExNx/fp1/PrrrzoPkKgqOGc8++HTY4cPH0a/fv0QGhoKAKhfvz5+/PFHHDt2DMDjWv1ly5Zh5syZ6NevHwDgu+++g6enJ7Zu3YpBgwYZLHai6sBkioiIzI3WNfvdunXD33//jQEDBiAjIwMZGRkYOHAgUlNT0aVLF33ESFQpqr7q0zaliH3WzYFCmYuEC2laHY9cZosefp5M9Guwjh07IiEhAX///TcA4PTp0zh06BB69+4NALhy5QoUCgWCg4PF58hkMgQEBCAxMbHM/ebl5SEzM1NtISIiIiLD07pmHwC8vb05EB8ZPXOcM151A+NWRi6cpFZYPqQNWtZxMXRYZAKmT5+OzMxMNGnSBJaWligqKsKiRYswdOhQAIBCoQAAeHp6qj3P09NT3Faa6OhozJs3T3+BExEREVGlaJTsp6SkoHnz5rCwsEBKSkq5ZVu2bKmTwMwJm5Ibhur9Vo1Cbw591c/dUuJWRi6y8wrw4GE+IuNOYsObHfm5ogpt2LAB69atQ1xcHJo1a4bk5GRMmjQJ3t7eCA8Pr/R+Z8yYgaioKPFxZmYmfHx8dBEyEREREVWBRsl+69atoVAo4OHhgdatW0MikaiN4KwikUhQVFSk8yBNmbbTnpFujQisD0CCzo1qm8X73sxbBiepFR48zIeNpQWycgvNosUC6d+UKVMwffp0se99ixYtcPXqVURHRyM8PBxyuRwAkJaWBi8vL/F5aWlpaN26dZn7lUqlkEqleo2diIiIiLSnUbJ/5coVuLu7i/8nzZljU3JT8PRNls6Nahs6JJ2Qy2yxfEgbRMadRFZuIZ6pZWcWLRZI/x4+fAgLC/VhWiwtLVFcXAwA8PX1hVwuR0JCgpjcZ2Zm4ujRoxg/fnx1h0tEREREVaRRsl+vXj3x/1evXkXHjh1hZaX+1MLCQhw+fFitLJlnU3JN6LPrgib7NuebLC3ruGDDmx3ZNYS00rdvXyxatAh169ZFs2bNcOrUKXz66acYNWoUgMctsyZNmoSFCxeiUaNG8PX1xaxZs+Dt7Y3+/fsbNngiIiIi0prWA/R1794dt2/fhoeHh9p6pVKJ7t27sxn/U2ritGf67Lqg6b7N/SaLXGZbIz5LpDvLly/HrFmz8PbbbyM9PR3e3t548803MXv2bLHM1KlTkZOTg3HjxiEjIwOdO3fGrl27YGvLzxoZB06PR0REpDmtk31BECCRSEqsv3fvHhwcHHQSlLmpaYmZPmvVNd23IW6ycCBGMmZOTk5YtmwZli1bVmYZiUSC+fPnY/78+dUXGBERERHphcbJ/sCBAwE8/jE4YsQItQGZioqKkJKSgo4dO+o+QjI5+qxVL23fZSXZFd1k0WVyzoEYiYiIiIjImGic7MtkjxM2QRDg5OQEOzs7cZuNjQ06dOiAsWPH6j5CMjn6rFV/et8AtEqyVQm+u5MU/9tyRmfJuTmPEUBERERERKZH42R/zZo1AID69evjvffeY5N9Kpc+uy48ue+EC2kaJ9lP1r472FghM7cAbg5SnSTn5j5GABERERERmRat++zPmTNHH3EQibRpXq9Ns/4na9/vP8yHk60VMh7l6yQ5r4kDMRIRERERkfHSKNl//vnnkZCQgFq1aqFNmzalDtCncvLkSZ0FRzWPtn3ftWnW38xbBjdHG9zOyIWXiy0+eqUV7mTl6Sw5r2kDMRIRERERkfHSKNnv16+fOCAf51smfdK07/vTtfdPNuu/lZELqZUFbmWU9nyJ+K+Hky1a1nEpc59ERERERESmSqNk/8mm+2zGT/qkSd/38mr/3Z2keJhfiAcPiyG1soC703998lNuZOD6/Ydws7fGvew8tRsBm5OuY+728xAEAfVqO3A0fSIiA+u+truhQyAiIjJpWvfZv379OiQSCerUqQMAOHbsGOLi4tC0aVOMGzdO5wFSzaJJ3/fyav//VmTDykICB3trPCoowrrEq0i+qcStB4+QlVcIAHiYX4R6bvbijYR9f6UhamOKuP+rd3M4mj4REREREZk0rZP9IUOGYNy4cRg2bBgUCgWCg4PRvHlzrFu3DgqFArNnz9ZHnFSDVNT3vbTa/5QbGfjqj3/w+7nbyCsEsvOKAAA/Jd0odR9pGY9w4fbjhH7zyZtq2yQSCUfTJyIiIiIik2ah7RPOnj2L9u3bAwA2bNiAFi1a4PDhw1i3bh1iY2N1HR9RCara/8VhLfH18Ha4cFuJgasOY/vpx4m+JnKLBIz9LgkpNzIw8Pln8OSQk6/6++glbjJPCmWuOAUkEREREZGx0DrZLygoEAfr27NnD15++WUAQJMmTXD79m3dRkf0hCeTKrnMFj38PJGelYuIuJMoLBa03l9hsYBZ287CzVGKiT2eReeGtVHX1RbbTt8UxwQgKo9q/Ihpm1L4mSEiIiIio6J1M/5mzZohJiYGoaGhiI+Px4IFCwAAt27dgpubm84DJPNRldHunxyUz83RBmM7N4SbozXeWZ+Mh/nFlY7p9HUlXolJhIONJZxtrZGVWwQ3B2m5MwEQqWg6ewQRERERUXXTumZ/8eLF+OKLLxAUFITBgwejVatWAIBffvlFbN5P9LSq1oAeungX1+49hK2VBS6mZWPm1jMY+90JZOZq2G6/HPmFxYAAKB8VwMnWChmP8sucCUB1LLpsts1m4KZLdeOqos8MERHpx8GDB9G3b194e3tDIpFg69atatsFQcDs2bPh5eUFOzs7BAcH4+LFi2pl7t+/j6FDh8LZ2RkuLi4YPXo0srOz1cqkpKSgS5cusLW1hY+PD5YsWaLvQyMiqjKta/aDgoJw9+5dZGZmolatWuL6cePGwd7eXqfBkfmoSg2oQpmLrw5dRk5+IbJyCyEAyCsqhqB9y/0yPXhUABsrCyx9vRWKBZTZ+qC8af+0pVDm4tDFu/jq0GXcy86v8v6o+mkyewQREelPTk4OWrVqhVGjRmHgwIElti9ZsgSff/451q5dC19fX8yaNQshISE4f/48bG0fX7OHDh2K27dvIz4+HgUFBRg5ciTGjRuHuLg4AEBmZiZ69uyJ4OBgxMTE4MyZMxg1ahRcXFw4ExURGTWtk30AsLS0RGFhIQ4dOgQAaNy4MerXr6/LuMjMlDaCfmlKa+p/7pYS97Lz4eNiD0VmLvKLilGJLvoVyi8sxj93cjC6S4Myy+iq2bbqpsG1ew+Rk1/4+NjYDNwkVTR7BBER6U/v3r3Ru3fvUrcJgoBly5Zh5syZ6NevHwDgu+++g6enJ7Zu3YpBgwbhwoUL2LVrF44fP462bdsCAJYvX44+ffrg448/hre3N9atW4f8/Hx8++23sLGxQbNmzZCcnIxPP/2UyT4RGTWtm/Hn5ORg1KhR8PLyQteuXdG1a1d4e3tj9OjRePjwoT5iJDPw9Aj6ZdWah685hsk/JSN8zTGxWbsq8c/OL4SnsxRSK60/thpLvpahFs/Tzet11WxbddPA1cEGEkhw/yGbgRMREenSlStXxGmiVWQyGQICApCYmAgASExMhIuLi5joA0BwcDAsLCxw9OhRsUzXrl1hY2MjlgkJCUFqaioePHhQ5uvn5eUhMzNTbSEiqk5aZ01RUVE4cOAAtm/fjoyMDGRkZGDbtm04cOAA3n33XX3ESGZCNYJ+WbWghy7exeX0bOTkFeJyejYOXbwrPu/r4e0Q3MQD1+8/wqOCyg/IV5HtZ24j5UZGiTEGUm5kIOFCGgBUeNNCE+INjLxCNPRwxOyXmrIJPxERkQ4pFAoAgKenp9p6T09PcZtCoYCHh4fadisrK7i6uqqVKW0fT75GaaKjoyGTycTFx4dT+xJR9dK6Gf+mTZvw888/IygoSFzXp08f2NnZ4bXXXsPq1at1GR/VKAIECAAk///vf231L9xWYv2JG9USxfRNpzGqU0Oxuf7NB48wIe4UcvILxRsPPfw8K95ROdjXm4iIyLzNmDEDUVFR4uPMzEwm/ERUrbRO9h8+fFji7iYAeHh4sBk/VUnnRu541sMRtzNy4eVii86N3MVtm0/erLY4zt/Oxke//wUJJFBk5kJmZ4WsvEK42ut2ejX29SYiItIfuVwOAEhLS4OXl5e4Pi0tDa1btxbLpKenqz2vsLAQ9+/fF58vl8uRlpamVkb1WFWmNFKpFFKptMrHQURUWVo34w8MDMScOXOQm/tfP+ZHjx5h3rx5CAwM1GlwVLPIZbZYOzIAS19vjbUjA9QS4YHPP1OtsaRl5uFOdh4e5hfCykICdyep2E/f3UnKqfKIiIiMnK+vL+RyORISEsR1mZmZOHr0qPibNTAwEBkZGUhKShLL7N27F8XFxQgICBDLHDx4EAUFBWKZ+Ph4NG7cWG1mKiIiY6N1zf5nn32GkJAQ1KlTB61atQIAnD59Gra2tti9e7fOA6Sapazabj8vGRylFsjO019//dIUC0B6Vh7e61kf9Ws7wN1Jiv9tOaOTqfeIiIioarKzs3Hp0iXx8ZUrV5CcnAxXV1fUrVsXkyZNwsKFC9GoUSNx6j1vb2/0798fAODn54devXph7NixiImJQUFBASIjIzFo0CB4e3sDAIYMGYJ58+Zh9OjRmDZtGs6ePYvPPvsMS5cuNcQhExFpTOtkv3nz5rh48SLWrVuHv/76CwAwePBgDB06FHZ2djoPkAh4PHp9YfXm+aL8IgEbkm5g7cj2Opt6j4iIiKruxIkT6N69u/hY1Uc+PDwcsbGxmDp1KnJycjBu3DhkZGSgc+fO2LVrF2xt//vuXrduHSIjI9GjRw9YWFggLCwMn3/+ubhdJpPh999/R0REBPz9/VG7dm3Mnj2b0+4RkdHTOtkHAHt7e4wdO1bXsRCVyd1JiuJioeKCOiS1lMDK0gLuTlLcy85TG0xPVbPPqfKIiIgMJygoCIJQ9u8DiUSC+fPnY/78+WWWcXV1RVxcXLmv07JlS/zxxx+VjpOIyBAqleynpqZi+fLluHDhAoDHTaAiIyPRpEkTnQZHpPK3IgtF1Zzs5xUJqO1kg+y8QjGx5yj6RERERERkCrQeoG/Tpk1o3rw5kpKS0KpVK7Rq1QonT55EixYtsGnTJn3EqBcrV65E/fr1YWtri4CAABw7dszQIVG5JDBEK/5X/OtgcVhLtb75cpktevh5Vluir1DmckBAIiIiIiLSitY1+1OnTsWMGTNKNIeaM2cOpk6dirCwMJ0Fpy8//fQToqKiEBMTg4CAACxbtgwhISFITU2Fh4eHocOjUjwnd4QFgKJqft2fjl9HDz8Pg9XgK5S5GPPdcZMfEFA1vgFbQxARERERVQ+ta/Zv376N4cOHl1j/xhtv4Pbt2zoJSt8+/fRTjB07FiNHjkTTpk0RExMDe3t7fPvtt4YOjcpwJysP1lZaf1yrTJGZi/E/JOHnEzcMUrNe2oCApkZ1w2LaphTxxgUZxs2bN/HGG2/Azc0NdnZ2aNGiBU6cOCFuFwQBs2fPhpeXF+zs7BAcHIyLFy8aMGIiIiIiqiyts6egoKBSByg5dOgQunTpopOg9Ck/Px9JSUkIDg4W11lYWCA4OBiJiYklyufl5SEzM1NtoepnIQFyCwwzHP/tjFzM/uUswtccq/ZEVVUTnvEo32QHBDSHGxbm4MGDB+jUqROsra3x22+/4fz58/jkk0/U5ohesmQJPv/8c8TExODo0aNwcHBASEgIcnN5g4aIiIjI1GjdjP/ll1/GtGnTkJSUhA4dOgAAjhw5go0bN2LevHn45Zdf1Moam7t376KoqAienp5q6z09PcWpBJ8UHR2NefPmVVd4VIZT1zIM8roWAIoB5BYU4XJ6Ng5dvItX2tap8Hm6arZuDgMCcgYD47B48WL4+PhgzZo14jpfX1/x/4IgYNmyZZg5cyb69esHAPjuu+/g6emJrVu3YtCgQaXuNy8vD3l5eeJj3hAlIiIiMg5aJ/tvv/02AGDVqlVYtWpVqduAx1OdFBVVdw9r3ZsxY4Y4Zyvw+Iesj4+PASOqmdrUdTHI63ZqWBuHr9wFBAkECADKnhFAleC7O0nxvy1ndNbPXi6zNckkX8UcbliYg19++QUhISF49dVXceDAATzzzDN4++23xWlUr1y5AoVCodbqSSaTISAgAImJiWUm+7whSkRERGSctE72i4sN05RaV2rXrg1LS0ukpaWprU9LS4NcLi9RXiqVQiqVVld4Jk9fA7Hdy87X2b40JXeWYkrvxpjycy5uZ+TCy8UWnRu5l1r2yYH07G0skZ1XBFf7/5qt1/QE19RvWJiDf/75B6tXr0ZUVBT+97//4fjx45g4cSJsbGwQHh4OhUIBAKW2elJtKw1viBIREREZJ62TfVNnY2MDf39/JCQkoH///gAe38BISEhAZGSkYYMzcfocOf76g0c62Y+mpFYW+HJ4W7Ss44K1IwMqvIHxZL/0ezl5cLa1Nul+9mR+iouL0bZtW3zwwQcAgDZt2uDs2bOIiYlBeHh4pffLG6JERERExqnGJfsAEBUVhfDwcLRt2xbt27fHsmXLkJOTg5EjRxo6NJNW2kBslUn2FcpcHLp4F4CAzo3c/7+/d/Um++3q10LLOi4ANKuVfrJf+jO17PDBgBa4k5XHZutkNLy8vNC0aVO1dX5+fti0aRMAiC2b0tLS4OXlJZZJS0tD69atqy1OIiIiItKNGpnsv/7667hz5w5mz54NhUKB1q1bY9euXSWar5J2dDEQm0KZi/A1x3A5PRsCBDzr4Yi1IwPwnKeTHiIu2/X7D8Xj0AT7pZOx69SpE1JTU9XW/f3336hXrx6Ax4P1yeVyJCQkiMl9ZmYmjh49ivHjx1d3uERERERURRon+7du3YK3t7c+Y6lWkZGRbLavY7pIeM/dUuJ2xiMIggBAgtsZj1sIhLb0xjd//oNbGXkV7kMXbjx4hEMX7+CVtpr3PdZXv3R9jYNANcvkyZPRsWNHfPDBB3jttddw7NgxfPnll/jyyy8BPB5UddKkSVi4cCEaNWoEX19fzJo1C97e3mKXJyIiIiIyHRaaFmzWrBni4uL0GQuZAbnMFj38PCudlDbzlqG2k/T/x70X4OViC3cnKc7dUqLLs6UPjqcPxQIASKrt9cqiGgdh2qYUcTwEospo164dtmzZgh9//BHNmzfHggULsGzZMgwdOlQsM3XqVEyYMAHjxo1Du3btkJ2djV27dsHWljeZiIiIiEyNxjX7ixYtwptvvoktW7bgiy++gKurqz7johrM2tIC9jZWcJRaoXPD2hj/QxIe5hfByday2mIQALg5Wlfb65VFV+MgEAHASy+9hJdeeqnM7RKJBPPnz8f8+fOrMSoiIiIi0geNa/bffvttpKSk4N69e2jatCm2b9+uz7iohjp3S4l72Xmo7WCD9KxcfPvnv7iZkYucvELcyymAjWX11bafupZR5jaFMhcJF9L0XtOuarrPkf2JiIiIiEgbWg3Q5+vri71792LFihUYOHAg/Pz8YGWlvouTJ0/qNECqWVTJ7bV7DwEAEgkgCEBBkYDaTtYoLgYUmfpvym5jKUFw09IHbNTnFINP48B/RERERERUGVqPxn/16lVs3rwZtWrVQr9+/Uok+0Rl0WSgOVVye+jiHaw6cBlX7+XAUiKBu6MUozv5ooG7A97deBr3cwr0Fqergw0+ebWlOPXe06q7ab2+Bv4jIiIiIiLzpVWm/tVXX+Hdd99FcHAwzp07B3f36hswjUybNrXhcpktXmnrg86N3HHo4h0oHxVgQ9J1rD5wGXKZLT55tRUm/HgK2XlFOo/T3sYS/+vdBN2blD0Noy6mGCQiIiIiItInjZP9Xr164dixY1ixYgWGDx+uz5jIDFWmNlyV9CdcSMO97HzxuX/8fU8viT4AFBQV4zm5U4XlRgTWByBB50a1WetORERERERGR+Nkv6ioCCkpKahTp44+4yEzVZXa8Kef+0wt/SXXhcUC/lZkl9mE/+kWCp0b1dZbLERERERERJWlcbIfHx+vzzjIzGk70NzT/ftV/fgBCZ6TO8LZ1gqZuYU6j/Px9BRCmds5FR4REREREZkCjq5H1UbTgeYUylyErzmG2xmP4OVih7Uj2wMAYhP/FWvU+zSXY/2JGzqLzdZKAktLC9SpZYfOjcoei4L99TWjyWCMRERERESkP0z2yegcungXl9OzIQgCLqdn49DFu6jlYK1Wo96+niskKK8OXjNWFhI842KL5UOex52svAqTU06FV7HqnJqQiIiIiIhKx2SfjJAAAQIAyf//K4iJ9c0Hj+Bka4Uuz9XGn//cxc0HjyCzs0ZTL2fEX0jX6lUcbCwwunMDBDf11CjRV+FUeOVjVwciIiIiIsNjsk9Gp3Mjdzzr4YjbGbnwcrFF50bukMts8cGAFpgQdwpZeYX4JP5vfPRKKzFJB4Bztw9BocyDBYAioexaf3sbC1hbWMDHzR7BTT3xvy1nWAutQ+zqQERERERkeEz2yejIZbZYOzKgRFP5O1l5yMkvhKv94xrjO1l56OHnKT5v8/jO4nPSs3Kx7sg1HPg7HRmP8pFbIMBSAjxTyw4rnmiyX1YtNPucVx67OhARERERGR6TfTJKpTWVr6jG+MnnyGW2aPmKC1JuZGDrqZtI+CsdD3LyYWdjBQ8nW7Wp9Z7eJ/ucVx27OhARERERGRaTfTIZqqb8e86nIbipp0bT9/1vyxlcu/cQD/OLUKeWHe5l56n1IS+tFjrhQhr7nBMRERERkUljsk8mQ6HMxaSfTuF2Ri52nLmFuDGB5Sbhqib6rvY2yMl/iPs5+ajrZl9uiwCAfc6JiIiIiMj0Mdknk7Ez5Rb+ufMQAPDPnYfYmXILo7s0KLP8k0n7sx6OGNu5ITo3ql1hLT37nBMRERERkaljsk8m4+aD3HIfq6j66T/jYocPBrQodVq9igbgY59zIiIiIiIyZUz2yWRGnu//vDfWHbuK/MJi2FhZoP/z3iXKpNzIwKsxh5FX+HjivYbuDlg3pkOJRJ8D8BERERERkTljsl/D6TLx1fdNg5Z1XLDxrUBxgL4nR9RX2XM+TUz0AeB2KQPslTXdHhERERERkblgsl/D6Srxra7a8pZ1XEpN8lWCm3rii4OXxYTfq5QB9jgAHxERERERmTsm+zWcrhJfY6ktf1z731Hssx/a0rtEHByAj4iIiIiIzB2T/RpOV4mvMdWWV1T7D3AAPiIiIiIiMm9M9kkniS9ry4mIiIiIiIwHk33SGVOvLU+5kVHu4H9EREREZHq6r+1u6BCIDMLC0AEQ6YNCmYuEC2lQKHM1Kp9yIwOvfZGI5fsu4bUvEpFyI6NaXpeIiIiIiEgfmOyT2VHNDDBtU4o4Q0BFHk/ZVwwrCwnyCoux53xatbwukaF8+OGHkEgkmDRpkrguNzcXERERcHNzg6OjI8LCwpCWpv25QERERESGx2SfROZSK13azAAVCW7qCamVBQqLBUitLBDc1LNaXpfIEI4fP44vvvgCLVu2VFs/efJkbN++HRs3bsSBAwdw69YtDBw40EBREhEREVFVsM8+AfivVlo1mv7Xw9uZbP/7yswM0LKOCza8GVilPvvGNCMBUVmys7MxdOhQfPXVV1i4cKG4XqlU4ptvvkFcXBxeeOEFAMCaNWvg5+eHI0eOoEOHDoYKmYiIiIgqgck+ASi9VtpUk/3KzgygyZR9+nhdouoUERGB0NBQBAcHqyX7SUlJKCgoQHBwsLiuSZMmqFu3LhITE8tM9vPy8pCXlyc+zszM1F/wRERERKQxJvsEwPxqpQ01M4Cpz0hA5m39+vU4efIkjh8/XmKbQqGAjY0NXFxc1NZ7enpCoVCUuc/o6GjMmzdP16ESERERURUx2ScArJUmMnfXr1/HO++8g/j4eNja6u78njFjBqKiosTHmZmZ8PHx0dn+iYiIiKhymOyTiLXSROYrKSkJ6enpeP7558V1RUVFOHjwIFasWIHdu3cjPz8fGRkZarX7aWlpkMvlZe5XKpVCKpXqM3QiIiIiqgSOxk81mrnMQFBVfB/MX48ePXDmzBkkJyeLS9u2bTF06FDx/9bW1khISBCfk5qaimvXriEwMNCAkRMRERFRZbBmn2osc5qBoCr4PtQMTk5OaN68udo6BwcHuLm5ietHjx6NqKgouLq6wtnZGRMmTEBgYCBH4iciIiIyQUz2qcY6d0uJGw8eQQLgxoNHJj0DQVWY00wMVDVLly6FhYUFwsLCkJeXh5CQEKxatcrQYRERERFRJTDZpxrLQgIoHxZAACD5/8c1kbnNxECa279/v9pjW1tbrFy5EitXrjRMQERERESkM2bVZ79+/fqQSCRqy4cffqhWJiUlBV26dIGtrS18fHywZMkSA0VLlaWr/uWnrmVAAGBpAQj//7gmUjXdXxzWkk34iYiIiIjMhFkl+wAwf/583L59W1wmTJggbsvMzETPnj1Rr149JCUl4aOPPsLcuXPx5ZdfGjBi0oaqf/m0TSliP/PKCm7qCVtrCxQLgK21BYKbeuowUtMil9mih58nE30iIjI7c+fOLVEZ1KRJE3F7bm4uIiIi4ObmBkdHR4SFhSEtLU1tH9euXUNoaCjs7e3h4eGBKVOmoLCwsLoPhYhIK2bXjN/JyanMaaLWrVuH/Px8fPvtt7CxsUGzZs2QnJyMTz/9FOPGjavmSKkydNm/vGUdF2x4MxB7zqchuKknWtZx0W2wREREZBSaNWuGPXv2iI+trP77CTx58mTs3LkTGzduhEwmQ2RkJAYOHIg///wTwONpSkNDQyGXy3H48GHcvn0bw4cPh7W1NT744INqPxYiIk2ZXc3+hx9+CDc3N7Rp0wYfffSR2l3XxMREdO3aFTY2NuK6kJAQpKam4sGDB6XuLy8vD5mZmWpLTWJsU7Kp+pdnPMrXSf/ylnVcENWzMRN9IiIiM2ZlZQW5XC4utWvXBgAolUp88803+PTTT/HCCy/A398fa9asweHDh3HkyBEAwO+//47z58/jhx9+QOvWrdG7d28sWLAAK1euRH5+fpmvWdN/QxKR4ZlVsj9x4kSsX78e+/btw5tvvokPPvgAU6dOFbcrFAp4eqo31VY9VigUpe4zOjoaMplMXHx8fPR3AEZGl03mdYX9y4mIiEhbFy9ehLe3Nxo0aIChQ4fi2rVrAICkpCQUFBQgODhYLNukSRPUrVsXiYmJAB5XFrVo0ULtN2RISAgyMzNx7ty5Ml+zJv+GJCLjYPTJ/vTp00v0s3p6+euvvwAAUVFRCAoKQsuWLfHWW2/hk08+wfLly5GXl1fp158xYwaUSqW4XL9+XVeHZvRKazJvDNi/nIiIiDQVEBCA2NhY7Nq1C6tXr8aVK1fQpUsXZGVlQaFQwMbGBi4uLmrP8fT0FCuCKlNZBNTs35BEZByMvs/+u+++ixEjRpRbpkGDBqWuDwgIQGFhIf799180btwYcrm8xIArqsdl9fOXSqWQSqXaB24GOCUbERERmbrevXuL/2/ZsiUCAgJQr149bNiwAXZ2dnp73Zr8G5KIjIPRJ/vu7u5wd3ev1HOTk5NhYWEBDw8PAEBgYCDef/99FBQUwNraGgAQHx+Pxo0bo1atWjqL2Vyomsyfu6UUE38iIiIiU+bi4oLnnnsOly5dwosvvoj8/HxkZGSo1e6npaWJFUFyuRzHjh1T20dFlUVERMbA6JvxayoxMRHLli3D6dOn8c8//2DdunWYPHky3njjDTGRHzJkCGxsbDB69GicO3cOP/30Ez777DNERUUZOHrjxSbzREREZE6ys7Nx+fJleHl5wd/fH9bW1khISBC3p6am4tq1awgMDATwuLLozJkzSE9PF8vEx8fD2dkZTZs2rfb4iYg0ZfQ1+5qSSqVYv3495s6di7y8PPj6+mLy5MlqibxMJsPvv/+OiIgI+Pv7o3bt2pg9ezan3aMKqcYsYAsHIiIi0/Lee++hb9++qFevHm7duoU5c+bA0tISgwcPhkwmw+jRoxEVFQVXV1c4OztjwoQJCAwMRIcOHQAAPXv2RNOmTTFs2DAsWbIECoUCM2fOREREBJvpm7Hua7sDAPaF7zNwJESVZzbJ/vPPPy9OkVKeli1b4o8//qiGiMhcqGYlUI1dwFkAiIiITMeNGzcwePBg3Lt3D+7u7ujcuTOOHDkidhNdunQpLCwsEBYWhry8PISEhGDVqlXi8y0tLbFjxw6MHz8egYGBcHBwQHh4OObPn2+oQyIi0ojZJPtE+lLarARM9omIiEzD+vXry91ua2uLlStXYuXKlWWWqVevHn799Vddh0ZEpFdm02efSF9UTfczHuVzVgIiIiIiIjIJrNknqgBnJSAiIiIiIlPDZJ9IA3KZLZN8IiIiIiIyGWzGT0RERERERGRmmOwTERERERERmRkm+0RERERERERmhsk+ERERERERkZlhsk9ERERERERkZpjsExEREREREZkZJvtEREREREREZobJPhEREREREZGZYbJPREREREREZGaY7BMRERERERGZGSb7REQ1QHR0NNq1awcnJyd4eHigf//+SE1NVSuTm5uLiIgIuLm5wdHREWFhYUhLSzNQxERERERUFUz2iYhqgAMHDiAiIgJHjhxBfHw8CgoK0LNnT+Tk5IhlJk+ejO3bt2Pjxo04cOAAbt26hYEDBxowaiIiIiKqLCtDB0BERPq3a9cutcexsbHw8PBAUlISunbtCqVSiW+++QZxcXF44YUXAABr1qyBn58fjhw5gg4dOhgibCIiIiKqJNbsExHVQEqlEgDg6uoKAEhKSkJBQQGCg4PFMk2aNEHdunWRmJhY5n7y8vKQmZmpthARERGR4THZJyKqYYqLizFp0iR06tQJzZs3BwAoFArY2NjAxcVFraynpycUCkWZ+4qOjoZMJhMXHx8ffYZORERERBpisk9khBTKXCRcSINCmWvoUMgMRURE4OzZs1i/fn2V9zVjxgwolUpxuX79ug4iJCIiIqKqYp99MloKZS7O3VKimbcMcpmtocOpNgplLsZ8dxwKZS7kMlt8PbxdjTp+0q/IyEjs2LEDBw8eRJ06dcT1crkc+fn5yMjIUKvdT0tLg1wuL3N/UqkUUqlUnyETERERUSWwZp+MkirhnbYpRUx8a4pzt5RQKHPhYmcj3vAgqipBEBAZGYktW7Zg79698PX1Vdvu7+8Pa2trJCQkiOtSU1Nx7do1BAYGVne4RERERFRFrNkno1RawltTardVLRlUNfvNvGWGDonMQEREBOLi4rBt2zY4OTmJ/fBlMhns7Owgk8kwevRoREVFwdXVFc7OzpgwYQICAwM5Ej8RERGRCWKyT0apJie8qqb7NbELA+nP6tWrAQBBQUFq69esWYMRI0YAAJYuXQoLCwuEhYUhLy8PISEhWLVqVTVHSkRERES6wGSfjFJNT3jlMtsad8ykX4IgVFjG1tYWK1euxMqVK6shIiIiIiLSJyb7NZCpDHzHhJeIiIiIiKhymOzXMBzpnYiIiIiIyPxxNP4ahiO9ExERERERmT8m+zWMqul+xqP8GjfwHRERERERUU3BZvw1TE0f+I6IiIiIiKgmYLJfA3HgOyIiIiIiIvPGZvxEeDxwYcKFNCiUuYYOhYiIiIiIqMpYs081HmcoICIiIiIic8OafarxOEMBERERERGZGyb7VONxhgIiIiIiIjI3bMZPNR5nKCAiIiIiInPDZJ8InKGAiIiIiIjMi8k041+0aBE6duwIe3t7uLi4lFrm2rVrCA0Nhb29PTw8PDBlyhQUFhaqldm/fz+ef/55SKVSPPvss4iNjdV/8ERERERERETVyGSS/fz8fLz66qsYP358qduLiooQGhqK/Px8HD58GGvXrkVsbCxmz54tlrly5QpCQ0PRvXt3JCcnY9KkSRgzZgx2795dXYdBREREREREpHcm04x/3rx5AFBmTfzvv/+O8+fPY8+ePfD09ETr1q2xYMECTJs2DXPnzoWNjQ1iYmLg6+uLTz75BADg5+eHQ4cOYenSpQgJCamuQ6EKqEbEZ/95IiIiIiKiyjGZmv2KJCYmokWLFvD09BTXhYSEIDMzE+fOnRPLBAcHqz0vJCQEiYmJZe43Ly8PmZmZagvpj2rO+2mbUjDmu+NQKHMNHRIREREREZHJMZtkX6FQqCX6AMTHCoWi3DKZmZl49OhRqfuNjo6GTCYTFx8fHz1ETyqc856IiIiIiKjqDJrsT58+HRKJpNzlr7/+MmSImDFjBpRKpbhcv37doPGYO855T0REREREVHUG7bP/7rvvYsSIEeWWadCggUb7ksvlOHbsmNq6tLQ0cZvqX9W6J8s4OzvDzs6u1P1KpVJIpVKNYqCq45z3REREREREVWfQZN/d3R3u7u462VdgYCAWLVqE9PR0eHh4AADi4+Ph7OyMpk2bimV+/fVXtefFx8cjMDBQJzGQbnDOeyIiIiIyBt3Xdhf/vy98nwEjIdKeyfTZv3btGpKTk3Ht2jUUFRUhOTkZycnJyM7OBgD07NkTTZs2xbBhw3D69Gns3r0bM2fOREREhFgz/9Zbb+Gff/7B1KlT8ddff2HVqlXYsGEDJk+ebMhDIyIiIiIiItIpk5l6b/bs2Vi7dq34uE2bNgCAffv2ISgoCJaWltixYwfGjx+PwMBAODg4IDw8HPPnzxef4+vri507d2Ly5Mn47LPPUKdOHXz99decdo+IiIiIiIjMiskk+7GxsYiNjS23TL169Uo0039aUFAQTp06pcPIiIiIiIiIiIyLyTTjJyIiIiIiMpTua7ur9eEnMnZM9omIiIiIiIjMjMk04yciIiIiIjK00mr3nxypX9Paf21H91ftl7MCkKZYs09ERGpWrlyJ+vXrw9bWFgEBATh27JihQyIiMgq8PlJZVE38tWnm/+Rz2EWA9IE1+0REJPrpp58QFRWFmJgYBAQEYNmyZQgJCUFqaio8PDwMHR4RkcHw+kjV4ena+4puAJRX2//kc9kaoGaSCIIgGDoIU5KZmQmZTAalUglnZ2dDh0NEJsiYryMBAQFo164dVqxYAQAoLi6Gj48PJkyYgOnTp1f4/MocG2syiMybNkkGr4/qmKxVDb9fNFPRjYLyylH1qcw1hDX7WlLdG8nMzDRwJERkqlTXD2O715qfn4+kpCTMmDFDXGdhYYHg4GAkJiaW+py8vDzk5eWJj5VKJQDtrpGFjworGTERmQJtrge8Pqp78vrI357a4/eLZrrEdNFpOV3bOWSn+P/QuNAq7+fJfTy576epypVXpjpV5vrIZF9LWVlZAAAfHx8DR0JEpi4rKwsymczQYYju3r2LoqIieHp6qq339PTEX3/9VepzoqOjMW/evBLreY0kIhXZeO2vc7w+llSZ95HIHOjqs1/afjTZt7Gde9pcH5nsa8nb2xvXr1+Hk5MTJBKJocNRk5mZCR8fH1y/ft3omr4ZCt+TkvielFTd74kgCMjKyoK3t7feX0vfZsyYgaioKPFxcXEx7t+/Dzc3N42ukebyeTSH4zCHYwB4HMZG2+Pg9fE/5vIZAMznWHgcxsdcjkWT46jM9ZHJvpYsLCxQp04dQ4dRLmdnZ5P+sOsD35OS+J6UVJ3viTHVWKnUrl0blpaWSEtLU1uflpYGuVxe6nOkUimkUqnaOhcXF61f21w+j+ZwHOZwDACPw9hocxy8Pqozl88AYD7HwuMwPuZyLBUdh7bXR069R0REAAAbGxv4+/sjISFBXFdcXIyEhAQEBgYaMDIiIsPi9ZGITBFr9omISBQVFYXw8HC0bdsW7du3x7Jly5CTk4ORI0caOjQiIoPi9ZGITA2TfTMilUoxZ86cEk3GajK+JyXxPSmJ78l/Xn/9ddy5cwezZ8+GQqFA69atsWvXrhKDUumKubz35nAc5nAMAI/D2JjLcQC8PlaFuRwLj8P4mMux6Os4JIKxzW1CRERERERERFXCPvtEREREREREZobJPhEREREREZGZYbJPREREREREZGaY7BMRERERERGZGSb7JmjRokXo2LEj7O3t4eLiUmqZa9euITQ0FPb29vDw8MCUKVNQWFioVmb//v14/vnnIZVK8eyzzyI2Nlb/wVejlStXon79+rC1tUVAQACOHTtm6JD06uDBg+jbty+8vb0hkUiwdetWte2CIGD27Nnw8vKCnZ0dgoODcfHiRbUy9+/fx9ChQ+Hs7AwXFxeMHj0a2dnZ1XgUuhMdHY127drByckJHh4e6N+/P1JTU9XK5ObmIiIiAm5ubnB0dERYWBjS0tLUymhyLlHlmPo5qslnzBR9+OGHkEgkmDRpkqFD0drNmzfxxhtvwM3NDXZ2dmjRogVOnDhh6LA0VlRUhFmzZsHX1xd2dnZo2LAhFixYAGMfS1kX3z81lbbXwY0bN6JJkyawtbVFixYt8Ouvv1ZTpBXT5lhiY2MhkUjUFltb22qMtnQVfZZLY4y/p7U9jv3795f4e0gkEigUiuoJuAyV/Z41tvOkMsehq3OEyb4Jys/Px6uvvorx48eXur2oqAihoaHIz8/H4cOHsXbtWsTGxmL27NlimStXriA0NBTdu3dHcnIyJk2ahDFjxmD37t3VdRh69dNPPyEqKgpz5szByZMn0apVK4SEhCA9Pd3QoelNTk4OWrVqhZUrV5a6fcmSJfj8888RExODo0ePwsHBASEhIcjNzRXLDB06FOfOnUN8fDx27NiBgwcPYty4cdV1CDp14MABRERE4MiRI4iPj0dBQQF69uyJnJwcsczkyZOxfft2bNy4EQcOHMCtW7cwcOBAcbsm5xJVjjmco5p8xkzN8ePH8cUXX6Bly5aGDkVrDx48QKdOnWBtbY3ffvsN58+fxyeffIJatWoZOjSNLV68GKtXr8aKFStw4cIFLF68GEuWLMHy5csNHVq5dPH9UxNpex08fPgwBg8ejNGjR+PUqVPo378/+vfvj7Nnz1Zz5CVV5pru7OyM27dvi8vVq1erMeLSVfRZfpqx/p7W9jhUUlNT1f4mHh4eeopQM5X5njXG86Syvxd0co4IZLLWrFkjyGSyEut//fVXwcLCQlAoFOK61atXC87OzkJeXp4gCIIwdepUoVmzZmrPe/3114WQkBC9xlxd2rdvL0RERIiPi4qKBG9vbyE6OtqAUVUfAMKWLVvEx8XFxYJcLhc++ugjcV1GRoYglUqFH3/8URAEQTh//rwAQDh+/LhY5rfffhMkEolw8+bNaotdX9LT0wUAwoEDBwRBeHz81tbWwsaNG8UyFy5cEAAIiYmJgiBodi5R5ZjjOfr0Z8zUZGVlCY0aNRLi4+OFbt26Ce+8846hQ9LKtGnThM6dOxs6jCoJDQ0VRo0apbZu4MCBwtChQw0UkfYq8/1TU2l7HXzttdeE0NBQtXUBAQHCm2++qdc4NaHtsZT1G9aYPP1ZLo0p/J7W5Dj27dsnABAePHhQLTFVlibfs8Z8nqhochy6OkdYs2+GEhMT0aJFC3h6eorrQkJCkJmZiXPnzollgoOD1Z4XEhKCxMTEao1VH/Lz85GUlKR2fBYWFggODjaL46uMK1euQKFQqL0nMpkMAQEB4nuSmJgIFxcXtG3bViwTHBwMCwsLHD16tNpj1jWlUgkAcHV1BQAkJSWhoKBA7T1p0qQJ6tatq/aeVHQukfbM9Rx9+jNmaiIiIhAaGlriu8FU/PLLL2jbti1effVVeHh4oE2bNvjqq68MHZZWOnbsiISEBPz9998AgNOnT+PQoUPo3bu3gSOrPE2+f2qiylwHjfW3W2Wv6dnZ2ahXrx58fHzQr18/k/xeNda/SWW1bt0aXl5eePHFF/Hnn38aOpwSNPmeNYW/iaa/F3RxjjDZN0MKhUItOQEgPlb1vSmrTGZmJh49elQ9gerJ3bt3UVRUVOrxGbrvkaGojru890ShUJRormVlZQVXV1eTf9+Ki4sxadIkdOrUCc2bNwfw+HhtbGxKjHvx9HtS0blE2jPHc7S0z5gpWb9+PU6ePIno6GhDh1Jp//zzD1avXo1GjRph9+7dGD9+PCZOnIi1a9caOjSNTZ8+HYMGDUKTJk1gbW2NNm3aYNKkSRg6dKihQ6s0Tb5/aqLKXAfL+k4y9PtYmWNp3Lgxvv32W2zbtg0//PADiouL0bFjR9y4caM6QtYZc/k97eXlhZiYGGzatAmbNm2Cj48PgoKCcPLkSUOHJtL0e9ZYzxMVTY9DV+eIVVUDJt2YPn06Fi9eXG6ZCxcuoEmTJtUUEZH5iIiIwNmzZ3Ho0CFDh0JmypQ/Y9evX8c777yD+Ph4oxggq7KKi4vRtm1bfPDBBwCANm3a4OzZs4iJiUF4eLiBo9PMhg0bsG7dOsTFxaFZs2ZiH2Bvb2+TOQYiTQQGBiIwMFB83LFjR/j5+eGLL77AggULDBhZzdS4cWM0btxYfNyxY0dcvnwZS5cuxffff2/AyP5jyt+zT9L0OHR1jjDZNxLvvvsuRowYUW6ZBg0aaLQvuVxeYgRU1Qjjcrlc/PfpUcfT0tLg7OwMOzs7DaM2TrVr14alpWWpx6c6/ppGddxpaWnw8vIS16elpaF169ZimacH0iksLMT9+/dN+n2LjIwUBxusU6eOuF4ulyM/Px8ZGRlqtftPfk40OZdIe+Z2jpb1GTMVSUlJSE9Px/PPPy+uKyoqwsGDB7FixQrk5eXB0tLSgBFqxsvLC02bNlVb5+fnh02bNhkoIu1NmTJFrN0HgBYtWuDq1auIjo422WRfk++fmqgy18GyfrsZ+rqpi2u6qiXLpUuX9BGi3pjz7+n27dsbTWKtzfessZ4nQNV+L1T2HGEzfiPh7u6OJk2alLvY2NhotK/AwECcOXNGLXGLj4+Hs7Oz+EMoMDAQCQkJas+Lj49Xu4NkqmxsbODv7692fMXFxUhISDCL46sMX19fyOVytfckMzMTR48eFd+TwMBAZGRkICkpSSyzd+9eFBcXIyAgoNpjripBEBAZGYktW7Zg79698PX1Vdvu7+8Pa2trtfckNTUV165dU3tPKjqXSHvmco5W9BkzFT169MCZM2eQnJwsLm3btsXQoUORnJxsEok+AHTq1KnEVEZ///036tWrZ6CItPfw4UNYWKj/NLO0tERxcbGBIqo6Tb5/aqLKXAeN9bebLq7pRUVFOHPmjNoNIVNgrH8TXUhOTjb436My37PG+DfRxe+FSp8jVR7ij6rd1atXhVOnTgnz5s0THB0dhVOnTgmnTp0SsrKyBEEQhMLCQqF58+ZCz549heTkZGHXrl2Cu7u7MGPGDHEf//zzj2Bvby9MmTJFuHDhgrBy5UrB0tJS2LVrl6EOS6fWr18vSKVSITY2Vjh//rwwbtw4wcXFRW1UdXOTlZUlfhYACJ9++qlw6tQp4erVq4IgCMKHH34ouLi4CNu2bRNSUlKEfv36Cb6+vsKjR4/EffTq1Uto06aNcPToUeHQoUNCo0aNhMGDBxvqkKpk/PjxgkwmE/bv3y/cvn1bXB4+fCiWeeutt4S6desKe/fuFU6cOCEEBgYKgYGB4nZNziWqHHM4RzX5jJkqUxyN/9ixY4KVlZWwaNEi4eLFi8K6desEe3t74YcffjB0aBoLDw8XnnnmGWHHjh3ClStXhM2bNwu1a9cWpk6daujQyqWL75+aqKLr4LBhw4Tp06eL5f/880/ByspK+Pjjj4ULFy4Ic+bMEaytrYUzZ84Y6hBE2h7LvHnzhN27dwuXL18WkpKShEGDBgm2trbCuXPnDHUIgiBU/FmePn26MGzYMLG8sf6e1vY4li5dKmzdulW4ePGicObMGeGdd94RLCwshD179hjqEARB0Ox71hTOk8och67OESb7Jig8PFwAUGLZt2+fWObff/8VevfuLdjZ2Qm1a9cW3n33XaGgoEBtP/v27RNat24t2NjYCA0aNBDWrFlTvQeiZ8uXLxfq1q0r2NjYCO3btxeOHDli6JD0SjVtytNLeHi4IAiPpz+aNWuW4OnpKUilUqFHjx5Camqq2j7u3bsnDB48WHB0dBScnZ2FkSNHijeRTE1p7wUAtc/5o0ePhLfffluoVauWYG9vLwwYMEC4ffu22n40OZeockz9HNXkM2aqTDHZFwRB2L59u9C8eXNBKpUKTZo0Eb788ktDh6SVzMxM4Z133hHq1q0r2NraCg0aNBDef/99o5/qUxffPzVVedfBbt26ie+hyoYNG4TnnntOsLGxEZo1aybs3LmzmiMumzbHMmnSJLGsp6en0KdPH+HkyZMGiFpdRZ/l8PBwoVu3biWeY2y/p7U9jsWLFwsNGzYUbG1tBVdXVyEoKEjYu3evYYJ/gibfs6ZwnlTmOHR1jkj+PwAiIiIiIiIiMhPss09ERERERERkZpjsExEREREREZkZJvtEREREREREZobJPhEREREREZGZYbJPREREREREZGaY7BMRERERERGZGSb7RERERERERGaGyT4RERERERGRmWGyTwQgNjYWLi4uFZaTSCTYunWrXmOZO3cuWrdurdfXICKqyIgRI9C/f/9qf11Nr8dERNWF10MyVRJBEARDB0E1Q1FREbp06QK5XI7NmzeL65VKJZo3b47hw4dj0aJFBont0aNHyMrKgoeHB4DHCffWrVuRnJysVk6hUKBWrVqQSqV6iyU7Oxt5eXlwc3PT22sQUc0mkUjK3T5nzhxMnjwZgiBU+w/Np6/HRET6xOshmTMrQwdANYelpSViY2PRunVrrFu3DkOHDgUATJgwAa6urpgzZ47BYrOzs4OdnV2F5eRyud5jcXR0hKOjo95fh4hqrtu3b4v//+mnnzB79mykpqaK6wx5HdL0ekxEpAu8HpI5YzN+qlbPPfccPvzwQ0yYMAG3b9/Gtm3bsH79enz33XewsbEp83n169fHggULMHjwYDg4OOCZZ57BypUr1cpcu3YN/fr1g6OjI5ydnfHaa68hLS1N3H769Gl0794dTk5OcHZ2hr+/P06cOAFAvZlUbGws5s2bh9OnT0MikUAikSA2NhZAyWb8Z86cwQsvvAA7Ozu4ublh3LhxyM7OFrermn19/PHH8PLygpubGyIiIlBQUFDmsT7djL8q+/j2229Rt25dODo64u2330ZRURGWLFkCuVwODw+PEi0pJBIJvvjiC7z00kuwt7eHn58fEhMTcenSJQQFBcHBwQEdO3bE5cuXy3xtIjJ+crlcXGQyGSQSido6R0fHEs1Wg4KCMGHCBEyaNAm1atWCp6cnvvrqK+Tk5GDkyJFwcnLCs88+i99++03ttc6ePYvevXvD0dERnp6eGDZsGO7evVtmbE83W1Vdz77//nvUr18fMpkMgwYNQlZWVpn7uHr1Kvr27YtatWrBwcEBzZo1w6+//lrp94uIzJe5Xw9V+9i9ezf8/Pzg6OiIXr16qd3k2L9/P9q3bw8HBwe4uLigU6dOuHr1qvZvJhkdJvtU7SZMmIBWrVph2LBhGDduHGbPno1WrVpV+LyPPvoIrVq1wqlTpzB9+nS88847iI+PBwAUFxejX79+uH//Pg4cOID4+Hj8888/eP3118XnDx06FHXq1MHx48eRlJSE6dOnw9rausTrvP7663j33f9r705DourCOID/UwgZk7QFKRCHitEZm8SwIg3FEhdEikylxCUSTdSasoIBQUMhDSmijChslCARoQWMVKJRQcdGEZdsKjNzAStSRNSS1PN+eOnSuKRt+jrv/wfz4Tz3LPf44cFz7zkzaXBzc0N/fz/6+/vN+vlmdHQUgYGBcHBwQENDA0pLS/HkyROkpKSY1dPr9ejs7IRer0dRUREKCwulhwcL9St9dHZ24vHjxygvL0dxcTEKCgoQEhKCvr4+VFdXIzc3F+np6Xj27JlZu6ysLMTExKC5uRmurq44cuQIEhMTodVq0djYCCHEjDkS0f9DUVER1q1bB6PRiNTUVCQlJSE8PBxeXl5oampCQEAAoqOjMTY2BgAYGhrC3r174eHhgcbGRpSXl+PDhw+IiIj4qXE7Ozvx4MEDlJWVoaysDNXV1cjJyZmzfnJyMsbHx1FTU4O2tjbk5uZyxxQR/VHLJR8CwNjYGPLy8nDnzh3U1NSgp6cHZ86cAQBMTEzgwIED8PX1RWtrKwwGAxISEuY93kDLhCBaAiaTSQAQarVafP36dd76zs7OIigoyCwWGRkpgoODhRBCVFZWCmtra9HT0yNdb29vFwCE0WgUQghhZ2cnCgsLZ+1fp9OJ1atXS+WMjAzh7u4+ox4Acf/+fSGEEDdv3hQODg5iZGREuv7o0SNhZWUl3r9/L4QQIjY2Vjg7O4uJiQmpTnh4uIiMjJxzrtPH/tU+ZDKZGB4elmKBgYFCLpeLyclJKebi4iIuXLhgNr/09HSpbDAYBABRUFAgxYqLi4WNjc2cYxPR8jI9/30TGxsr9u/fL5V9fX3Fnj17pPLExISwtbUV0dHRUqy/v18AEAaDQQghRFZWlggICDDrt7e3VwAQr169WtD9zJbPzp49K3bt2jXnnNRqtcjMzJzzOhHRbCwxH+p0OgFAvHnzRorl5+cLR0dHIYQQAwMDAoCoqqqasw9avvhmn5bE7du3IZPJ0NXVhb6+vgW12b1794yyyWQCAJhMJjg5OcHJyUm6rlKpYG9vL9U5ffo04uPj4e/vj5ycnN/eim4ymeDu7g5bW1sp5u3tjampKbOzXm5ubrC2tpbKGzZswMePH39qrF/pQy6Xw87OTio7OjpCpVLBysrKLDa9n23btpldBwC1Wm0W+/LlC4aHh39qDkS0/H2fH6ytrbF27doZ+QGAlFdaWlqg1+ulM6+rVq2Cq6srAPxUDp6ez+bLgSdOnEB2dja8vb2RkZGB1tbWBY9FRLQQyyUfAoBMJsPmzZtnbbNmzRrExcUhMDAQoaGhuHLlitkWf1reuNinRVdXV4fLly+jrKwMO3fuxLFjxyAW4UchMjMz0d7ejpCQEDx9+hQqlQr379//6+NOPyqwYsUKTE1N/fU+ZmuzkH6+r/NtC9dssZ+dAxEtf/Pllen5YWRkBKGhoWhubjb7dHR0wMfH57fG/VEOio+Px9u3bxEdHY22tjZ4enri6tWrCx6PiGg+yyUfztXm+/+9dTodDAYDvLy8UFJSAoVCgfr6+gXfE/13cbFPi2psbAxxcXFISkqCn58fCgoKYDQacePGjXnbTk869fX1UCqVAAClUone3l709vZK11+8eIGhoSGoVCopplAocOrUKVRWVuLgwYPQ6XSzjrVy5UpMTk7+8H6USiVaWlowOjoqxWpra2FlZQUXF5d550NEZOm2b9+O9vZ2yOVybNmyxezz/a6ov8HJyQnHjx/HvXv3kJaWhlu3bv3V8YiIfmQp8+FCeHh4QKvVoq6uDlu3bsXdu3eX+pboD+BinxaVVquFEEL6IhG5XI68vDycO3cO7969+2Hb2tpaXLx4Ea9fv0Z+fj5KS0tx8uRJAIC/vz/UajWioqLQ1NQEo9GImJgY+Pr6wtPTE58/f0ZKSgqqqqrQ3d2N2tpaNDQ0SA8LppPL5ejq6kJzczM+ffqE8fHxGXWioqJgY2OD2NhYPH/+HHq9HqmpqYiOjpa2bhER/Z8lJydjcHAQhw8fRkNDAzo7O1FRUYGjR4/O+0D1d2g0GlRUVKCrqwtNTU3Q6/Vz5nsiosWwVPlwPl1dXdBqtTAYDOju7kZlZSU6OjqYMy0EF/u0aKqrq5Gfnw+dTgeZTCbFExMT4eXlNe92/rS0NDQ2NsLDwwPZ2dm4dOkSAgMDAfy7Henhw4dwcHCAj48P/P39sWnTJpSUlAD49yzVwMAAYmJioFAoEBERgeDgYJw/f37WscLCwhAUFAQ/Pz+sX78excXFM+rIZDJUVFRgcHAQO3bswKFDh7Bv3z5cu3btd/5MREQWY+PGjaitrcXk5CQCAgKgVquh0Whgb29v9v0hf9rk5CSSk5OhVCoRFBQEhUKB69ev/7XxiIjms1T5cD4ymQwvX75EWFgYFAoFEhISkJycjMTExCW7J/pzVojFOCxN9Jvkcjk0Gg00Gs1S3woREREREdF/Ht/sExEREREREVkYLvaJiIiIiIiILAy38RMRERERERFZGL7ZJyIiIiIiIrIwXOwTERERERERWRgu9omIiIiIiIgsDBf7RERERERERBaGi30iIiIiIiIiC8PFPhEREREREZGF4WKfiIiIiIiIyMJwsU9ERERERERkYf4BCJxLRZVbSUUAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1200x400 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(1, 3, figsize=(12, 4))\n",
    "\n",
    "# Plot 1\n",
    "n = 2000 # use only the n first elements\n",
    "hits_u = hits.arrays()['PostPosition_X'][:n]\n",
    "hits_v = hits.arrays()['PostPosition_Y'][:n]\n",
    "ax[0].scatter(hits_u, hits_v, alpha=0.8, s=3)\n",
    "ax[0].set_title('Positions (2D) of hits in the camera')\n",
    "ax[0].set_xlabel('X position in mm')\n",
    "ax[0].set_ylabel('Y position in mm')\n",
    "\n",
    "# Plot 2\n",
    "times = hits.arrays()['GlobalTime']/1e9\n",
    "ax[1].hist(times, 100, facecolor='green', alpha=0.75)\n",
    "ax[1].set_title('Hits time (global)')\n",
    "ax[1].set_xlabel('Time in s')\n",
    "\n",
    "# Plot 3\n",
    "times = hits.arrays()['LocalTime']\n",
    "ax[2].hist(times, 100, facecolor='green', alpha=0.75)\n",
    "ax[2].set_title('Hits time (local)')\n",
    "ax[2].set_xlabel('Time in ns')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "cell_id": "00007-d230344f-c2fe-4569-be4d-1a75fdad9663",
    "output_cleared": false
   },
   "source": [
    "# Plot 2: edep, step length and track length"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:10.202919Z",
     "start_time": "2024-01-23T19:56:10.085840Z"
    },
    "cell_id": "00008-e2b8bf16-0053-4df6-85a9-00bdca67e379",
    "execution_millis": 2701,
    "execution_start": 1605009145833,
    "output_cleared": false,
    "source_hash": "fb76ead1"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'Length in mm')"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA+IAAAGJCAYAAAAOton7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABm6UlEQVR4nO3deVxV1f7/8fcBBBQFRIMjhWhqTjmFQ5QTSaKSZWl9LUsyrzZAOZSZXTOHijJT02vaqFSaZteszFRS1FKcKJwzNU1vCt4yISxBZf3+6Me+HkFlPAfl9Xw89uPhWftz9l5ryVmcD3vttW3GGCMAAAAAAOAUbq6uAAAAAAAAFQmJOAAAAAAATkQiDgAAAACAE5GIAwAAAADgRCTiAAAAAAA4EYk4AAAAAABORCIOAAAAAIATkYgDAAAAAOBEJOIAAAAAADgRiTiKbPXq1bLZbFq9erWrqwInGjt2rGw2m6urAQCWzp076/rrr3d1NQBcJh588EFVrVrV1dVwujp16ui2225zdTVwHhLxy4TNZivUVpjk+KWXXtLixYvLvM5z5sy5aF03bNhQ5nVwpXPb6uHhoYCAAIWFhWnIkCHatWuXq6tXKpz1swRcjrZv364+ffooNDRU3t7euvrqq3Xrrbdq+vTpDnGu/hzljdVbtmxxWR0u5siRIxo7dqxSU1NdXRUAxVCa32FdKe+CxK+//urqqhRo165dGjt2rA4ePOjqqqCQPFxdARTOBx984PD6/fffV2JiYr7yxo0bX/JYL730kvr06aNevXqVZhUvaPz48apbt26+8vr16zvl/K506623qn///jLGKCMjQ1u3blVCQoLeeOMNvfLKKxo+fLirq1hoo0eP1jPPPONQ5uyfJeBysX79ekVERKh27doaNGiQ7Ha7Dh8+rA0bNuj111/X448/bsXyObq4I0eOaNy4capTp45atmzp6uoAKKLS/A6LC9u1a5fGjRunzp07q06dOq6uDgqBRPwycf/99zu83rBhgxITE/OVl0fdu3dX69atXV0NnTx5Uj4+Pk4953XXXZfv/+jll19Wz5499eSTT6pRo0bq0aOHU+tUXB4eHvLwYMgACuPFF1+Un5+fNm/eLH9/f4d9x44dc02lAMAFivsd9s8//1SVKlXKsmqASzE1/Qpy8uRJPfnkkwoJCZGXl5caNmyoSZMmyRhjxdhsNp08eVIJCQnWVKAHH3xQkvTzzz/rscceU8OGDVW5cmXVqFFDd999d5lPcTl48KBsNpsmTZqkt956S/Xq1ZOXl5fatGmjzZs354v/4Ycf1KdPHwUEBMjb21utW7fW559/7hCTN9VyzZo1euyxxxQYGKhrrrnG2j9jxgxde+21qly5stq2batvvvlGnTt3VufOnSVJWVlZ8vHx0ZAhQ/Kd/z//+Y/c3d0VHx9frPbWqFFD8+fPl4eHh1588UWHfdnZ2Xr++edVv359eXl5KSQkRE8//bSys7Md4mw2m+Li4jR37lw1bNhQ3t7eCgsL09q1a/Od7/vvv1f37t3l6+urqlWrqkuXLvluCzh9+rTGjRunBg0ayNvbWzVq1FD79u2VmJhoxZx/j/jFfpaAim7//v1q2rRpviRckgIDA61/X+pz9Msvv+ihhx5SUFCQvLy81LRpU7333nsOx8tbt2PBggV69tlnZbfb5ePjo9tvv12HDx8utTYVpS4ff/yxXnzxRV1zzTXy9vZWly5dtG/fvnzHvNRYvHr1arVp00aSNGDAAKuP5syZ43CcXbt2KSIiQlWqVNHVV1+tiRMnllq7AZS9vPUeUlJS1LFjR1WpUkXPPvusJOmzzz5TdHS0goOD5eXlpXr16mnChAk6e/ZsvuNs3LhRPXr0UPXq1eXj46PmzZvr9ddfv+i5U1NTddVVV6lz587KysoqcVuK8j113bp1Gj58uK666ir5+Pjozjvv1H//+1+H2NzcXI0dO1bBwcGqUqWKIiIitGvXLtWpU8f6fTFnzhzdfffdkqSIiIgLTvf/9ttv1bZtW3l7e+vaa6/V+++/X+L2ovi4vHWFMMbo9ttvV1JSkgYOHKiWLVtq+fLlGjFihH755RdNmTJF0t/Tg/7xj3+obdu2Gjx4sCSpXr16kqTNmzdr/fr16tu3r6655hodPHhQM2fOVOfOnbVr165i/1UyIyMj3/00NptNNWrUcCibN2+e/vjjDz388MOy2WyaOHGi7rrrLv3000+qVKmSJGnnzp26+eabdfXVV+uZZ56Rj4+PPv74Y/Xq1Uv//ve/deeddzoc87HHHtNVV12lMWPG6OTJk5KkmTNnKi4uTh06dNCwYcN08OBB9erVS9WrV7eS9apVq+rOO+/UggULNHnyZLm7u1vH/Oijj2SMUb9+/YrVH5JUu3ZtderUSUlJScrMzJSvr69yc3N1++2369tvv9XgwYPVuHFjbd++XVOmTNGPP/6Y7x7SNWvWaMGCBXriiSfk5eWlN954Q926ddOmTZusxYt27typDh06yNfXV08//bQqVaqkN998U507d9aaNWvUrl07SX8n2fHx8dbPRmZmprZs2aLvvvtOt956a4FtuNjPElDRhYaGKjk5WTt27LjoYmIX+xylp6frxhtvtP7wdtVVV+mrr77SwIEDlZmZqaFDhzoc68UXX5TNZtPIkSN17NgxTZ06VZGRkUpNTVXlypVL1J6i1uXll1+Wm5ubnnrqKWVkZGjixInq16+fNm7caMUUZixu3Lixxo8frzFjxmjw4MHq0KGDJOmmm26yjvP777+rW7duuuuuu3TPPffok08+0ciRI9WsWTN17969RO0G4Dy//fabunfvrr59++r+++9XUFCQpL+TzKpVq2r48OGqWrWqVq1apTFjxigzM1Ovvvqq9f7ExETddtttqlWrloYMGSK73a7du3dryZIlBV5Ykf7+7hsVFaXWrVvrs88+K/FYWdTvqY8//riqV6+u559/XgcPHtTUqVMVFxenBQsWWDGjRo3SxIkT1bNnT0VFRWnr1q2KiorSqVOnrJiOHTvqiSee0LRp0/Tss89a0/zPne6/b98+9enTRwMHDlRMTIzee+89PfjggwoLC1PTpk1L1G4Uk8FlKTY21pz737d48WIjybzwwgsOcX369DE2m83s27fPKvPx8TExMTH5jvnnn3/mK0tOTjaSzPvvv2+VJSUlGUkmKSnponWcPXu2kVTg5uXlZcUdOHDASDI1atQwx48ft8o/++wzI8l88cUXVlmXLl1Ms2bNzKlTp6yy3Nxcc9NNN5kGDRrkO3f79u3NmTNnrPLs7GxTo0YN06ZNG3P69GmrfM6cOUaS6dSpk1W2fPlyI8l89dVXDu1q3ry5Q9yFSDKxsbEX3D9kyBAjyWzdutUYY8wHH3xg3NzczDfffOMQN2vWLCPJrFu3zuHYksyWLVussp9//tl4e3ubO++80yrr1auX8fT0NPv377fKjhw5YqpVq2Y6duxolbVo0cJER0dftD3PP/+8OX/IuNDPElDRrVixwri7uxt3d3cTHh5unn76abN8+XKTk5OTL/ZCn6OBAweaWrVqmV9//dWhvG/fvsbPz88as/PG5KuvvtpkZmZacR9//LGRZF5//fWL1jVvvNy8efMFY4pal8aNG5vs7Gwr7vXXXzeSzPbt240xRRuLN2/ebCSZ2bNn56tXp06d8v2Oys7ONna73fTu3fui7QbgGud/hzXmf5/lWbNm5Ysv6Pvpww8/bKpUqWJ9Hzxz5oypW7euCQ0NNb///rtDbG5urvXvmJgY4+PjY4wx5ttvvzW+vr4mOjra4XvlheR9D/rvf/97wZiifk+NjIx0qN+wYcOMu7u7OXHihDHGmLS0NOPh4WF69erlcJ6xY8caSQ6/OxYuXHjB7+ehoaFGklm7dq1VduzYMePl5WWefPLJS7YdZYOp6VeIpUuXyt3dXU888YRD+ZNPPiljjL766qtLHuPcvwKePn1av/32m+rXry9/f3999913xa7bjBkzlJiY6LAVVJ//+7//U/Xq1a3XeVc+fvrpJ0nS8ePHtWrVKt1zzz36448/9Ouvv+rXX3/Vb7/9pqioKO3du1e//PKLwzEHDRrkcDV7y5Yt+u233zRo0CCH+5379evncG5JioyMVHBwsObOnWuV7dixQ9u2bSuVe/PzHp/xxx9/SJIWLlyoxo0bq1GjRlbbfv31V91yyy2SpKSkJIf3h4eHKywszHpdu3Zt3XHHHVq+fLnOnj2rs2fPasWKFerVq5euvfZaK65WrVq677779O233yozM1OS5O/vr507d2rv3r0lbheAvxdqTE5O1u23366tW7dq4sSJioqK0tVXX51vimJBjDH697//rZ49e8oY4zAmREVFKSMjI9+43L9/f1WrVs163adPH9WqVUtLly4tUVuKU5cBAwbI09PTen3+eF6UsfhSqlat6jAme3p6qm3btta5AFwevLy8NGDAgHzl534/zfv+16FDB/3555/64YcfJP19G96BAwc0dOjQfLcEFfTo1aSkJEVFRalLly5atGiRvLy8Slz/4nxPHTx4sEP9OnTooLNnz+rnn3+WJK1cuVJnzpzRY4895vC+cxf8LKwmTZpYY7EkXXXVVWrYsCFjpQsxNf0K8fPPPys4ONjhS5j0vykpeR/oi/nrr78UHx+v2bNn65dffnG4tzwjI6PYdWvbtm2hFmurXbu2w+u8L2O///67pL+n1Bhj9Nxzz+m5554r8BjHjh3T1Vdfbb0+f7X2vH44f8V2Dw+PfCtMurm5qV+/fpo5c6a1YMjcuXPl7e1t3YdTEnn3IeX9n+3du1e7d+/WVVddVWD8+Qs8NWjQIF/Mddddpz///NO6v+jPP/9Uw4YN88U1btxYubm5Onz4sJo2barx48frjjvu0HXXXafrr79e3bp10wMPPKDmzZuXqI1ARdamTRstWrRIOTk52rp1qz799FNNmTJFffr0UWpqqpo0aXLB9/73v//ViRMn9NZbb+mtt94qMOZSY4LNZlP9+vVLvM5HcepyqfG8KGPxpVxzzTX5vmhXr15d27ZtK9JxALjW1Vdf7fAHvDw7d+7U6NGjtWrVKusCQp6876f79++XpIveCpTn1KlTio6OVlhYmD7++ONSW4i2ON9TiztWBgQEFPmPluefK+98eeeC85GIw/L4449r9uzZGjp0qMLDw+Xn5yebzaa+ffsqNze3zM9/7pXrc+X9QSCvDk899ZSioqIKjD1/oCrpvT79+/fXq6++qsWLF+vee+/VvHnzdNttt8nPz69Ex5X+vrru7u5u/bEgNzdXzZo10+TJkwuMDwkJKfE5L6Rjx47av3+/PvvsM61YsULvvPOOpkyZolmzZukf//hHmZ0XqAg8PT3Vpk0btWnTRtddd50GDBighQsX6vnnn7/ge/LGu/vvv18xMTEFxjjrD2XFqculxvPS5MxzASg7BX1nO3HihDp16iRfX1+NHz9e9erVk7e3t7777juNHDmyWN9Pvby81KNHD3322WdatmyZbrvtttKofrG+pzJWVmwk4leI0NBQff311/rjjz8crornTdkJDQ21ygqaoiNJn3zyiWJiYvTaa69ZZadOndKJEyfKptJFlDe9ulKlSoqMjCzWMfL6Yd++fYqIiLDKz5w5o4MHD+b7Mnn99derVatWmjt3rq655hodOnRI06dPL2YL/ufQoUNas2aNwsPDrf+vevXqaevWrerSpcsF/4/OVdA08h9//FFVqlSxrqpXqVJFe/bsyRf3ww8/yM3NzSG5DwgI0IABAzRgwABlZWWpY8eOGjt27EUT8cLUE8D/5M0OOnr0qFVW0OfoqquuUrVq1XT27NlCj3fnjwnGGO3bt6/ECXtx6nIpRRmLGWeAimv16tX67bfftGjRInXs2NEqP3DggENc3iKXO3bsuOQ4ZbPZNHfuXN1xxx26++679dVXX1lPaiiJ0vieer5zx8pzZ3n+9ttv+a5kM1ZefrhH/ArRo0cPnT17Vv/6178cyqdMmSKbzeawcqyPj0+BybW7u3u+v4pNnz69wMdDuEJgYKA6d+6sN9980+FLbJ7zH/dQkNatW6tGjRp6++23debMGat87ty5F5ya88ADD2jFihWaOnWqatSoUeJVeI8fP657771XZ8+e1T//+U+r/J577tEvv/yit99+O997/vrrL2vV9zzJyckO92UePnxYn332mbp27Sp3d3e5u7ura9eu+uyzzxympqanp2vevHlq3769fH19Jf09oJ+ratWqql+/fr7Hpp3vQj9LQEWXlJRU4FWGvPu1z71lpKDPkbu7u3r37q1///vf2rFjR77jFDTevf/++9aaE9Lff1w9evRoices4tTlUooyFvv4+EgSYw1QAeVdxT13PM3JydEbb7zhEHfDDTeobt26mjp1ar6xoqCx2NPTU4sWLVKbNm3Us2dPbdq0qcR1LY3vqefr0qWLPDw8NHPmTIfy87/vS4yVlyOuiF8hevbsqYiICP3zn//UwYMH1aJFC61YsUKfffaZhg4d6vBYqbCwMH399deaPHmygoODVbduXbVr10633XabPvjgA/n5+alJkyZKTk7W119/ne8xY0X11VdfWVfmz3XTTTc5LCJWGDNmzFD79u3VrFkzDRo0SNdee63S09OVnJys//znP9q6detF3+/p6amxY8fq8ccf1y233KJ77rlHBw8e1Jw5c1SvXr0C/5p433336emnn9ann36qRx991HqUWmH8+OOP+vDDD2WMUWZmprZu3aqFCxcqKytLkydPVrdu3azYBx54QB9//LEeeeQRJSUl6eabb9bZs2f1ww8/6OOPP9by5csd7rW//vrrFRUV5fD4MkkaN26cFfPCCy8oMTFR7du312OPPSYPDw+9+eabys7OdnjObpMmTdS5c2eFhYUpICBAW7Zs0SeffKK4uLiLtu9CP0tARff444/rzz//1J133qlGjRopJydH69ev14IFC1SnTh2HBYku9Dl6+eWXlZSUpHbt2mnQoEFq0qSJjh8/ru+++05ff/21jh8/7nDOgIAAtW/fXgMGDFB6erqmTp2q+vXra9CgQYWq83vvvadly5blKx8yZEiR63IpRRmL69WrJ39/f82aNUvVqlWTj4+P2rVrl28NEABXnptuuknVq1dXTEyMnnjiCdlsNn3wwQf5kms3NzfNnDlTPXv2VMuWLTVgwADVqlVLP/zwg3bu3Knly5fnO3blypW1ZMkS3XLLLerevbvWrFlTqHvMJ0+enO+Rvm5ubnr22WdL/D31fEFBQRoyZIhee+013X777erWrZu2bt2qr776SjVr1nQYK1u2bCl3d3e98sorysjIkJeXl2655RYFBgYW6ZxwIucu0o7SUtCjH/744w8zbNgwExwcbCpVqmQaNGhgXn31VYfHIhhjzA8//GA6duxoKleu7PDog99//90MGDDA1KxZ01StWtVERUWZH374wYSGhjo8HqE0Hl+mcx5Fk/f4sldffTXfMSSZ559/3qFs//79pn///sZut5tKlSqZq6++2tx2223mk08+yXfuCz2OZ9q0aSY0NNR4eXmZtm3bmnXr1pmwsDDTrVu3AuN79OhhJJn169dftM3n1z1vc3NzM/7+/qZVq1ZmyJAhZufOnQW+Jycnx7zyyiumadOmxsvLy1SvXt2EhYWZcePGmYyMDIdjx8bGmg8//NA0aNDAeHl5mVatWhX4f/Ldd9+ZqKgoU7VqVVOlShUTERGRrx0vvPCCadu2rfH39zeVK1c2jRo1Mi+++KLDo5YKenzZhX6WgIruq6++Mg899JBp1KiRqVq1qvH09DT169c3jz/+uElPT3eIvdjnKD093cTGxpqQkBBTqVIlY7fbTZcuXcxbb71lxeSNyR999JEZNWqUCQwMNJUrVzbR0dHm559/vmRdLzVWHz58uMh1WbhwocM58sb58x9BVtix+LPPPjNNmjQxHh4eDsfp1KmTadq0ab42xcTEmNDQ0Eu2HYDzXejxZQV9lo0xZt26debGG280lStXNsHBwdbjIAv6Lvrtt9+aW2+91VSrVs34+PiY5s2bm+nTp1v7z318WZ5ff/3VNGnSxNjtdrN3794L1jvve1BBm7u7uxVXku+pBX3HPnPmjHnuueeM3W43lStXNrfccovZvXu3qVGjhnnkkUcc3v/222+ba6+91ri7uzscJzQ0tMDH1Hbq1KlQj+RF2bAZwx36QG5urq666irdddddBU4Nv/POO7V9+3bt27fPBbXLz2azKTY2tsCpSQAqltWrVysiIkILFy5Unz59XF2dErnUWAwA+Hv6efXq1fXCCy843OaIywv3iKPCOXXqVL4pTe+//76OHz9e4GIdR48e1ZdffqkHHnjASTUEgCtfUcdiAKiI/vrrr3xlU6dOlSTGyssc94ijwtmwYYOGDRumu+++WzVq1NB3332nd999V9dff73D88EPHDigdevW6Z133lGlSpX08MMPu7DWAHBlKexYDAAV2YIFCzRnzhz16NFDVatW1bfffquPPvpIXbt21c033+zq6qEESMRR4dSpU0chISGaNm2ajh8/roCAAPXv318vv/yyPD09rbg1a9ZowIABql27thISEmS3211YawC4shR2LAaAiqx58+by8PDQxIkTlZmZaS3g9sILL7i6aigh7hEHAAAAAMCJuEccAAAAAAAnIhEHAABAmYqPj1ebNm1UrVo1BQYGqlevXtqzZ49DTOfOnWWz2Ry2Rx55xCHm0KFDio6OVpUqVRQYGKgRI0bozJkzDjGrV6/WDTfcIC8vL9WvX19z5swp6+YBQJFdsfeI5+bm6siRI6pWrZrDw+4BoDCMMfrjjz8UHBwsN7cr62+WjI8ASqqoY+SaNWsUGxurNm3a6MyZM3r22WfVtWtX7dq1Sz4+PlbcoEGDNH78eOt1lSpVrH+fPXtW0dHRstvtWr9+vY4ePar+/furUqVKeumllyT9vdBqdHS0HnnkEc2dO1crV67UP/7xD9WqVUtRUVGXrCfjI4CSKvT46LInmJexw4cPG0lsbGxsJdoOHz7s6uGs1DE+srGxldZW3DHy2LFjRpJZs2aNVdapUyczZMiQC75n6dKlxs3NzaSlpVllM2fONL6+viY7O9sYY8zTTz9tmjZt6vC+//u//zNRUVGFqhfjIxsbW2ltlxofr9gr4tWqVZMkHT58WL6+vi6uDYDLTWZmpkJCQqyx5ErC+AigpEo6RmZkZEiSAgICHMrnzp2rDz/8UHa7XT179tRzzz1nXRVPTk5Ws2bNFBQUZMVHRUXp0Ucf1c6dO9WqVSslJycrMjLS4ZhRUVEaOnRogfXIzs5Wdna29dr8/zWMGR8BFFdhx8crNhHPm07k6+vLQAqg2K7EqYmMjwBKS3HGyNzcXA0dOlQ333yzrr/+eqv8vvvuU2hoqIKDg7Vt2zaNHDlSe/bs0aJFiyRJaWlpDkm4JOt1WlraRWMyMzP1119/qXLlyg774uPjNW7cuHx1ZHwEUFKXGh+v2EQcAAAA5U9sbKx27Nihb7/91qF88ODB1r+bNWumWrVqqUuXLtq/f7/q1atXJnUZNWqUhg8fbr3Ou5IFAGXtylqBCAAAAOVWXFyclixZoqSkJF1zzTUXjW3Xrp0kad++fZIku92u9PR0h5i813a7/aIxvr6++a6GS5KXl5d19Zur4ACciUQcAAAAZcoYo7i4OH366adatWqV6tate8n3pKamSpJq1aolSQoPD9f27dt17NgxKyYxMVG+vr5q0qSJFbNy5UqH4yQmJio8PLyUWgIApYNEHAAAAGUqNjZWH374oebNm6dq1aopLS1NaWlp+uuvvyRJ+/fv14QJE5SSkqKDBw/q888/V//+/dWxY0c1b95cktS1a1c1adJEDzzwgLZu3arly5dr9OjRio2NlZeXlyTpkUce0U8//aSnn35aP/zwg9544w19/PHHGjZsmMvaDgAFIREHAABAmZo5c6YyMjLUuXNn1apVy9oWLFggSfL09NTXX3+trl27qlGjRnryySfVu3dvffHFF9Yx3N3dtWTJErm7uys8PFz333+/+vfv7/Dc8bp16+rLL79UYmKiWrRooddee03vvPNOoZ4hDgDOxGJtAAAAKFN5jwW7kJCQEK1Zs+aSxwkNDdXSpUsvGtO5c2d9//33RaofADgbV8QBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCJWTQfKgYiECOvfSTFJLqwJcGF5P6f8jAIAnIXfPbhScUUcAAAAAAAnIhEHAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQAAAAAAnIhEHABKydq1a9WzZ08FBwfLZrNp8eLFF4x95JFHZLPZNHXqVIfy48ePq1+/fvL19ZW/v78GDhyorKwsh5ht27apQ4cO8vb2VkhIiCZOnFgGrQEAAEBZIREHgFJy8uRJtWjRQjNmzLho3KeffqoNGzYoODg4375+/fpp586dSkxM1JIlS7R27VoNHjzY2p+ZmamuXbsqNDRUKSkpevXVVzV27Fi99dZbpd4eAAAAlA0PV1cAAK4U3bt3V/fu3S8a88svv+jxxx/X8uXLFR0d7bBv9+7dWrZsmTZv3qzWrVtLkqZPn64ePXpo0qRJCg4O1ty5c5WTk6P33ntPnp6eatq0qVJTUzV58mSHhB0AAADlF1fEAcBJcnNz9cADD2jEiBFq2rRpvv3Jycny9/e3knBJioyMlJubmzZu3GjFdOzYUZ6enlZMVFSU9uzZo99//73A82ZnZyszM9NhAwAAgOuQiAOAk7zyyivy8PDQE088UeD+tLQ0BQYGOpR5eHgoICBAaWlpVkxQUJBDTN7rvJjzxcfHy8/Pz9pCQkJK2hQAAACUAIk4ADhBSkqKXn/9dc2ZM0c2m82p5x41apQyMjKs7fDhw049PwAAAByRiAOAE3zzzTc6duyYateuLQ8PD3l4eOjnn3/Wk08+qTp16kiS7Ha7jh075vC+M2fO6Pjx47Lb7VZMenq6Q0ze67yY83l5ecnX19dhAwAAgOuQiAOAEzzwwAPatm2bUlNTrS04OFgjRozQ8uXLJUnh4eE6ceKEUlJSrPetWrVKubm5ateunRWzdu1anT592opJTExUw4YNVb16dec2CgAAAMXCqukAUEqysrK0b98+6/WBAweUmpqqgIAA1a5dWzVq1HCIr1Spkux2uxo2bChJaty4sbp166ZBgwZp1qxZOn36tOLi4tS3b1/rUWf33Xefxo0bp4EDB2rkyJHasWOHXn/9dU2ZMsV5DQUAAECJkIgDQCnZsmWLIiIirNfDhw+XJMXExGjOnDmFOsbcuXMVFxenLl26yM3NTb1799a0adOs/X5+flqxYoViY2MVFhammjVrasyYMTy6DAAA4DJCIg4ApaRz584yxhQ6/uDBg/nKAgICNG/evIu+r3nz5vrmm2+KWj0AAACUE9wjDgAAAACAE5GIAwAAAADgRCTiAAAAAAA4EYk4AAAAAABORCIOAAAAAIATkYgDAAAAAOBEJOIAAAAAADgRiTgAAAAAAE5EIg4AAAAAgBORiAMAAAAA4EQk4gAAAAAAOJGHqysAAAAAABcTkRBh/TspJsmFNQFKB1fEAQAAAABwIhJxAAAAAACciEQcAAAAAAAnIhEHAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwImKlIjHx8erTZs2qlatmgIDA9WrVy/t2bPHIebUqVOKjY1VjRo1VLVqVfXu3Vvp6ekOMYcOHVJ0dLSqVKmiwMBAjRgxQmfOnHGIWb16tW644QZ5eXmpfv36mjNnTvFaCAAAAABAOVKkRHzNmjWKjY3Vhg0blJiYqNOnT6tr1646efKkFTNs2DB98cUXWrhwodasWaMjR47orrvusvafPXtW0dHRysnJ0fr165WQkKA5c+ZozJgxVsyBAwcUHR2tiIgIpaamaujQofrHP/6h5cuXl0KTAQAAAABwHY+iBC9btszh9Zw5cxQYGKiUlBR17NhRGRkZevfddzVv3jzdcsstkqTZs2ercePG2rBhg2688UatWLFCu3bt0tdff62goCC1bNlSEyZM0MiRIzV27Fh5enpq1qxZqlu3rl577TVJUuPGjfXtt99qypQpioqKKrBu2dnZys7Otl5nZmYWqSMAAAAAAHCGEt0jnpGRIUkKCAiQJKWkpOj06dOKjIy0Yho1aqTatWsrOTlZkpScnKxmzZopKCjIiomKilJmZqZ27txpxZx7jLyYvGMUJD4+Xn5+ftYWEhJSkqYBAAAAAFAmip2I5+bmaujQobr55pt1/fXXS5LS0tLk6ekpf39/h9igoCClpaVZMecm4Xn78/ZdLCYzM1N//fVXgfUZNWqUMjIyrO3w4cPFbRoAAAAAAGWm2Il4bGysduzYofnz55dmfYrNy8tLvr6+DhsAAABcjwV/AcBRsRLxuLg4LVmyRElJSbrmmmuscrvdrpycHJ04ccIhPj09XXa73Yo5f1DNe32pGF9fX1WuXLk4VQYAAICLsOAvADgqUiJujFFcXJw+/fRTrVq1SnXr1nXYHxYWpkqVKmnlypVW2Z49e3To0CGFh4dLksLDw7V9+3YdO3bMiklMTJSvr6+aNGlixZx7jLyYvGMAAADg8rFs2TI9+OCDatq0qVq0aKE5c+bo0KFDSklJkSRrwd/JkyfrlltuUVhYmGbPnq3169drw4YNkmQt+Pvhhx+qZcuW6t69uyZMmKAZM2YoJydHkhwW/G3cuLHi4uLUp08fTZkyxWVtB4CCFCkRj42N1Ycffqh58+apWrVqSktLU1pamnXftp+fnwYOHKjhw4crKSlJKSkpGjBggMLDw3XjjTdKkrp27aomTZrogQce0NatW7V8+XKNHj1asbGx8vLykiQ98sgj+umnn/T000/rhx9+0BtvvKGPP/5Yw4YNK+XmAwAAwNnKy4K/2dnZyszMdNgAwBmKlIjPnDlTGRkZ6ty5s2rVqmVtCxYssGKmTJmi2267Tb1791bHjh1lt9u1aNEia7+7u7uWLFkid3d3hYeH6/7771f//v01fvx4K6Zu3br68ssvlZiYqBYtWui1117TO++8c8FHlwFAebB27Vr17NlTwcHBstlsWrx4sbXv9OnTGjlypJo1ayYfHx8FBwerf//+OnLkiMMxjh8/rn79+snX11f+/v4aOHCgsrKyHGK2bdumDh06yNvbWyEhIZo4caIzmgcApaI8LfjLU3cAuEqRniNujLlkjLe3t2bMmKEZM2ZcMCY0NFRLly696HE6d+6s77//vijVAwCXOnnypFq0aKGHHnrI4b5GSfrzzz/13Xff6bnnnlOLFi30+++/a8iQIbr99tu1ZcsWK65fv346evSodQ/lgAEDNHjwYM2bN0+SlJmZqa5duyoyMlKzZs3S9u3b9dBDD8nf31+DBw92ansBoDjyFvz99ttvXV0VjRo1SsOHD7deZ2ZmkowDcIoiJeIAgAvr3r27unfvXuA+Pz8/JSYmOpT961//Utu2bXXo0CHVrl1bu3fv1rJly7R582a1bt1akjR9+nT16NFDkyZNUnBwsObOnaucnBy999578vT0VNOmTZWamqrJkyeTiAMo9/IW/F27du0FF/w996r4+Qv+btq0yeF4JV3w18vLy7o1EgCcqdiPLwMAlExGRoZsNpv1pTM5OVn+/v5WEi5JkZGRcnNz08aNG62Yjh07ytPT04qJiorSnj179Pvvvxd4Hu6BBOBqLPgLAI5IxAHABU6dOqWRI0fq3nvvla+vr6S/720MDAx0iPPw8FBAQECR7pE8H/dAAnA1FvwFAEck4gDgZKdPn9Y999wjY4xmzpxZ5ucbNWqUMjIyrO3w4cNlfk4AOBcL/gKAI+4Rh1NEJERY/06KSXJhTQDXykvCf/75Z61atcq6Gi79fW/juVMuJenMmTM6fvz4Je9/zNtXEO6BBOBqLPgLAI64Ig4ATpKXhO/du1dff/21atSo4bA/PDxcJ06cUEpKilW2atUq5ebmql27dlbM2rVrdfr0aSsmMTFRDRs2VPXq1Z3TEAAAAJQIiTgAlJKsrCylpqYqNTVVknTgwAGlpqbq0KFDOn36tPr06aMtW7Zo7ty5Onv2rHWPZE5OjiSpcePG6tatmwYNGqRNmzZp3bp1iouLU9++fRUcHCxJuu++++Tp6amBAwdq586dWrBggV5//XWHx+8AAACgfGNqOgCUki1btigi4n+3YeQlxzExMRo7dqw+//xzSVLLli0d3peUlKTOnTtLkubOnau4uDh16dJFbm5u6t27t6ZNm2bF+vn5acWKFYqNjVVYWJhq1qypMWPG8OgyAACAywiJOACUks6dO1/0PsjC3CMZEBCgefPmXTSmefPm+uabb4pcPwAAAJQPTE0HAAAAAMCJuCIOAAAAoNw492k7wJWKK+IAAAAAADgRiTgAAAAAAE5EIg4AAAAAgBORiAMAAAAA4EQk4gAAAAAAOBGJOAAAAAAATkQiDgAAAACAE5GIAwAAAADgRCTiAAAAAAA4EYk4AAAAAABORCIOAAAAAIATkYgDAAAAAOBEJOIAAAAAADgRiTgAoEgiEiIUkRDh6moAAABctkjEAQAAAABwIhJxAAAAAACciEQcAAAAAAAnIhEHAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQAAAAAAnIhEHAAAAAAAJyIRBwAAAADAiUjEAQAAAABwIhJxACgla9euVc+ePRUcHCybzabFixc77DfGaMyYMapVq5YqV66syMhI7d271yHm+PHj6tevn3x9feXv76+BAwcqKyvLIWbbtm3q0KGDvL29FRISookTJ5Z10wAAAFCKSMQBoJScPHlSLVq00IwZMwrcP3HiRE2bNk2zZs3Sxo0b5ePjo6ioKJ06dcqK6devn3bu3KnExEQtWbJEa9eu1eDBg639mZmZ6tq1q0JDQ5WSkqJXX31VY8eO1VtvvVXm7QMAAEDp8HB1BQDgStG9e3d17969wH3GGE2dOlWjR4/WHXfcIUl6//33FRQUpMWLF6tv377avXu3li1bps2bN6t169aSpOnTp6tHjx6aNGmSgoODNXfuXOXk5Oi9996Tp6enmjZtqtTUVE2ePNkhYQcAAED5xRVxAHCCAwcOKC0tTZGRkVaZn5+f2rVrp+TkZElScnKy/P39rSRckiIjI+Xm5qaNGzdaMR07dpSnp6cVExUVpT179uj3338v8NzZ2dnKzMx02AAAAOA6JOIA4ARpaWmSpKCgIIfyoKAga19aWpoCAwMd9nt4eCggIMAhpqBjnHuO88XHx8vPz8/aQkJCSt4gAAAAFBuJOABc4UaNGqWMjAxrO3z4sKurBAAAUKGRiAOAE9jtdklSenq6Q3l6erq1z26369ixYw77z5w5o+PHjzvEFHSMc89xPi8vL/n6+jpsAAAAcB0ScQBwgrp168put2vlypVWWWZmpjZu3Kjw8HBJUnh4uE6cOKGUlBQrZtWqVcrNzVW7du2smLVr1+r06dNWTGJioho2bKjq1as7qTUAAAAoCRJxACglWVlZSk1NVWpqqqS/F2hLTU3VoUOHZLPZNHToUL3wwgv6/PPPtX37dvXv31/BwcHq1auXJKlx48bq1q2bBg0apE2bNmndunWKi4tT3759FRwcLEm677775OnpqYEDB2rnzp1asGCBXn/9dQ0fPtxFrQYAAEBR8fgyACglW7ZsUUREhPU6LzmOiYnRnDlz9PTTT+vkyZMaPHiwTpw4ofbt22vZsmXy9va23jN37lzFxcWpS5cucnNzU+/evTVt2jRrv5+fn1asWKHY2FiFhYWpZs2aGjNmDI8uAwAAuIyQiANAKencubOMMRfcb7PZNH78eI0fP/6CMQEBAZo3b95Fz9O8eXN98803xa4nAAAAXIup6QAAAAAAOBGJOAAAAAAATlTkRHzt2rXq2bOngoODZbPZtHjxYof9Dz74oGw2m8PWrVs3h5jjx4+rX79+8vX1lb+/vwYOHKisrCyHmG3btqlDhw7y9vZWSEiIJk6cWPTWAQAAAABQzhQ5ET958qRatGihGTNmXDCmW7duOnr0qLV99NFHDvv79eunnTt3KjExUUuWLNHatWsdFhrKzMxU165dFRoaqpSUFL366qsaO3as3nrrraJWFwAAAACAcqXIi7V1795d3bt3v2iMl5eX7HZ7gft2796tZcuWafPmzWrdurUkafr06erRo4cmTZqk4OBgzZ07Vzk5OXrvvffk6emppk2bKjU1VZMnT2ZlYAAAAADAZa1M7hFfvXq1AgMD1bBhQz366KP67bffrH3Jycny9/e3knBJioyMlJubmzZu3GjFdOzYUZ6enlZMVFSU9uzZo99//73Ac2ZnZyszM9NhAwAAAACgvCn1RLxbt256//33tXLlSr3yyitas2aNunfvrrNnz0qS0tLSFBgY6PAeDw8PBQQEKC0tzYoJCgpyiMl7nRdzvvj4ePn5+VlbSEhIaTcNAAAAxcAaQwDgqNSfI963b1/r382aNVPz5s1Vr149rV69Wl26dCnt01lGjRql4cOHW68zMzNJxgEAAMqBvDWGHnroId11110FxnTr1k2zZ8+2Xnt5eTns79evn44eParExESdPn1aAwYM0ODBgzVv3jxJ/1tjKDIyUrNmzdL27dv10EMPyd/fn1sbAZQ7pZ6In+/aa69VzZo1tW/fPnXp0kV2u13Hjh1ziDlz5oyOHz9u3Vdut9uVnp7uEJP3+kL3nnt5eeUbsAEAAOB6rDEEAI7K/Dni//nPf/Tbb7+pVq1akqTw8HCdOHFCKSkpVsyqVauUm5urdu3aWTFr167V6dOnrZjExEQ1bNhQ1atXL+sqAwAAwMlYYwhARVLkRDwrK0upqalKTU2VJB04cECpqak6dOiQsrKyNGLECG3YsEEHDx7UypUrdccdd6h+/fqKioqSJDVu3FjdunXToEGDtGnTJq1bt05xcXHq27evgoODJUn33XefPD09NXDgQO3cuVMLFizQ66+/7jD1HAAAAFcG1hgCUNEUeWr6li1bFBERYb3OS45jYmI0c+ZMbdu2TQkJCTpx4oSCg4PVtWtXTZgwwWHa+Ny5cxUXF6cuXbrIzc1NvXv31rRp06z9fn5+WrFihWJjYxUWFqaaNWtqzJgxTCsCAAC4ArHGEICKpsiJeOfOnWWMueD+5cuXX/IYAQEB1sIaF9K8eXN98803Ra0eAAAALnOsMQTgSlfm94gDAAAARcEaQwCudCTiAAAAKFOsMQQAjkjEAQAAUKa2bNmiVq1aqVWrVpL+XmOoVatWGjNmjNzd3bVt2zbdfvvtuu666zRw4ECFhYXpm2++ybfGUKNGjdSlSxf16NFD7du311tvvWXtz1tj6MCBAwoLC9OTTz7JGkMAyq0yf444AAAAKjbWGAIAR1wRBwAAAADAiUjEAQAAAABwIhJxAAAAAACciEQcAAAAAAAnIhEHAAAAAMCJSMQBAAAAAHAiEnEAAAAAl42IhAhFJES4uhpAiZCIAwAAAADgRCTiAAAAAAA4EYk4AAAAAABORCIOAAAAAIATkYgDAAAAAOBEJOIA4CRnz57Vc889p7p166py5cqqV6+eJkyYIGOMFWOM0ZgxY1SrVi1VrlxZkZGR2rt3r8Nxjh8/rn79+snX11f+/v4aOHCgsrKynN0cAAAAFBOJOAA4ySuvvKKZM2fqX//6l3bv3q1XXnlFEydO1PTp062YiRMnatq0aZo1a5Y2btwoHx8fRUVF6dSpU1ZMv379tHPnTiUmJmrJkiVau3atBg8e7IomAQAAoBg8XF0BAKgo1q9frzvuuEPR0dGSpDp16uijjz7Spk2bJP19NXzq1KkaPXq07rjjDknS+++/r6CgIC1evFh9+/bV7t27tWzZMm3evFmtW7eWJE2fPl09evTQpEmTFBwc7JrGAQAAoNC4Ig4ATnLTTTdp5cqV+vHHHyVJW7du1bfffqvu3btLkg4cOKC0tDRFRkZa7/Hz81O7du2UnJwsSUpOTpa/v7+VhEtSZGSk3NzctHHjxgLPm52drczMTIcNAAAArsMVcQBwkmeeeUaZmZlq1KiR3N3ddfbsWb344ovq16+fJCktLU2SFBQU5PC+oKAga19aWpoCAwMd9nt4eCggIMCKOV98fLzGjRtX2s1RREKE9e+kmKRSPz4AAMCViiviAOAkH3/8sebOnat58+bpu+++U0JCgiZNmqSEhIQyPe+oUaOUkZFhbYcPHy7T8wEAAODiuCIOAE4yYsQIPfPMM+rbt68kqVmzZvr5558VHx+vmJgY2e12SVJ6erpq1aplvS89PV0tW7aUJNntdh07dszhuGfOnNHx48et95/Py8tLXl5eZdAiAAAAFAdXxAHASf7880+5uTkOu+7u7srNzZUk1a1bV3a7XStXrrT2Z2ZmauPGjQoPD5ckhYeH68SJE0pJSbFiVq1apdzcXLVr184JrQAAAEBJcUUcAJykZ8+eevHFF1W7dm01bdpU33//vSZPnqyHHnpIkmSz2TR06FC98MILatCggerWravnnntOwcHB6tWrlySpcePG6tatmwYNGqRZs2bp9OnTiouLU9++fVkxHQAA4DJBIg4ATjJ9+nQ999xzeuyxx3Ts2DEFBwfr4Ycf1pgxY6yYp59+WidPntTgwYN14sQJtW/fXsuWLZO3t7cVM3fuXMXFxalLly5yc3NT7969NW3aNFc0CQAAAMVAIg4ATlKtWjVNnTpVU6dOvWCMzWbT+PHjNX78+AvGBAQEaN68eWVQQwAAADgD94gDAAAAAOBEJOIAAAAAADgRiTgAAAAAAE5EIg4AAAAAgBORiAMAAAAA4EQk4gAAAEAJRCREKCIhwtXVAHAZIREHAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQAAAAAAnIhEHAAAAAAAJyIRBwAAAADAiUjEAQAAAABwIhJxAECJRSREKCIhwtXVAAAAuCyQiAMAAAAA4EQk4gAAAAAAOBGJOAAAAAAATkQiDgAAAACAE5GIAwAAAADgRCTiAAAAAAA4kYerK1DRnPt4n6SYJBfWBAAAAADgClwRBwAAAADAiYqciK9du1Y9e/ZUcHCwbDabFi9e7LDfGKMxY8aoVq1aqly5siIjI7V3716HmOPHj6tfv37y9fWVv7+/Bg4cqKysLIeYbdu2qUOHDvL29lZISIgmTpxY9NYBAJwqIiHC2gAAAFCwIifiJ0+eVIsWLTRjxowC90+cOFHTpk3TrFmztHHjRvn4+CgqKkqnTp2yYvr166edO3cqMTFRS5Ys0dq1azV48GBrf2Zmprp27arQ0FClpKTo1Vdf1dixY/XWW28Vo4kAAAAAAJQfRU7Eu3fvrhdeeEF33nlnvn3GGE2dOlWjR4/WHXfcoebNm+v999/XkSNHrCvnu3fv1rJly/TOO++oXbt2at++vaZPn6758+fryJEjkqS5c+cqJydH7733npo2baq+ffvqiSee0OTJk0vWWgAAADgdMyoBwFGp3iN+4MABpaWlKTIy0irz8/NTu3btlJycLElKTk6Wv7+/WrdubcVERkbKzc1NGzdutGI6duwoT09PKyYqKkp79uzR77//XuC5s7OzlZmZ6bABAADA9ZhRCQCOSnXV9LS0NElSUFCQQ3lQUJC1Ly0tTYGBgY6V8PBQQECAQ0zdunXzHSNvX/Xq1fOdOz4+XuPGjSudhgAAAKDUdO/eXd27dy9w3/kzKiXp/fffV1BQkBYvXqy+fftaMyo3b95sXcyZPn26evTooUmTJik4ONhhRqWnp6eaNm2q1NRUTZ482SFhB4Dy4IpZNX3UqFHKyMiwtsOHD7u6SgCQzy+//KL7779fNWrUUOXKldWsWTNt2bLF2l9a0zMB4HLBjEoAFVGpJuJ2u12SlJ6e7lCenp5u7bPb7Tp27JjD/jNnzuj48eMOMQUd49xznM/Ly0u+vr4OGwCUJ7///rtuvvlmVapUSV999ZV27dql1157zWGWT2lMzwSAy0lpzqgs6BjnnuN88fHx8vPzs7aQkJCSNwgACqFUE/G6devKbrdr5cqVVllmZqY2btyo8PBwSVJ4eLhOnDihlJQUK2bVqlXKzc1Vu3btrJi1a9fq9OnTVkxiYqIaNmxY4LR0ALgcvPLKKwoJCdHs2bPVtm1b1a1bV127dlW9evUkld6Cl+fjig8AFIwZlQBcpciJeFZWllJTU5Wamirp7+lEqampOnTokGw2m4YOHaoXXnhBn3/+ubZv367+/fsrODhYvXr1kiQ1btxY3bp106BBg7Rp0yatW7dOcXFx6tu3r4KDgyVJ9913nzw9PTVw4EDt3LlTCxYs0Ouvv67hw4eXWsMBwNk+//xztW7dWnfffbcCAwPVqlUrvf3229b+0pqeeT6u+AAoz5hRCaAiKnIivmXLFrVq1UqtWrWSJA0fPlytWrXSmDFjJElPP/20Hn/8cQ0ePFht2rRRVlaWli1bJm9vb+sYc+fOVaNGjdSlSxf16NFD7du3d1jR0s/PTytWrNCBAwcUFhamJ598UmPGjGHqJYDL2k8//aSZM2eqQYMGWr58uR599FE98cQTSkhIkFR60zPPxxUfAOUZMyoBVERFXjW9c+fOMsZccL/NZtP48eM1fvz4C8YEBARo3rx5Fz1P8+bN9c033xS1egBQbuXm5qp169Z66aWXJEmtWrXSjh07NGvWLMXExJTZeb28vOTl5VVmxweAS8nKytK+ffus13kzKgMCAlS7dm1rRmWDBg1Ut25dPffccxecUTlr1iydPn26wBmV48aN08CBAzVy5Ejt2LFDr7/+uqZMmeKKJgPARV0xq6YDQHlXq1YtNWnSxKGscePGOnTokKTSm54JAOUNMyoBwFGpPkccAHBhN998s/bs2eNQ9uOPPyo0NFSS4/TMli1bSvrf9MxHH31UkuP0zLCwMEn5p2cCQHnDjEoAcEQiDgBOMmzYMN1000166aWXdM8992jTpk166623rCs65y54WZLpmQAAACjfSMQBwEnatGmjTz/9VKNGjdL48eNVt25dTZ06Vf369bNinn76aZ08eVKDBw/WiRMn1L59+wKnZ8bFxalLly5yc3NT7969NW3aNFc0CQAAAMVAIg4ATnTbbbfptttuu+D+0pqeCQAAgPKLxdoAAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQAAAAAAnIhEHAAAAAAAJyIRBwAAAADAiUjEAQAAAABwIhJxAAAAAACciEQcAAAAAAAnIhEHAJSJiIQIRSREuLoaAAAA5Q6JOAAAAAAATkQiDgAAAACAE5GIAwAAAADgRCTiAAAAAAA4kYerKwAAuLKdu2BbUkySC2sCAABQPnBFHAAAAAAAJyIRBwAAAADAiUjEAQAAAABwIhJxAAAAAACciEQcAAAAAAAnIhEHAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQCA00QkRCgiIcLV1QAAAHApEnEAcJGXX35ZNptNQ4cOtcpOnTql2NhY1ahRQ1WrVlXv3r2Vnp7u8L5Dhw4pOjpaVapUUWBgoEaMGKEzZ844ufYAAAAoLhJxAHCBzZs3680331Tz5s0dyocNG6YvvvhCCxcu1Jo1a3TkyBHddddd1v6zZ88qOjpaOTk5Wr9+vRISEjRnzhyNGTPG2U0AAABAMZGIA4CTZWVlqV+/fnr77bdVvXp1qzwjI0PvvvuuJk+erFtuuUVhYWGaPXu21q9frw0bNkiSVqxYoV27dunDDz9Uy5Yt1b17d02YMEEzZsxQTk6Oq5oEAACAIiARBwAni42NVXR0tCIjIx3KU1JSdPr0aYfyRo0aqXbt2kpOTpYkJScnq1mzZgoKCrJioqKilJmZqZ07dxZ4vuzsbGVmZjpsAAAAcB0PV1cAACqS+fPn67vvvtPmzZvz7UtLS5Onp6f8/f0dyoOCgpSWlmbFnJuE5+3P21eQ+Ph4jRs3rhRqDwAAgNJAIg4ATnL48GENGTJEiYmJ8vb2dtp5R40apeHDh1uvMzMzFRIS4rTzAwBQGDxVAxUJiTgAOElKSoqOHTumG264wSo7e/as1q5dq3/9619avny5cnJydOLECYer4unp6bLb7ZIku92uTZs2ORw3b1X1vJjzeXl5ycvLq5RbUzLnftlKiklyYU0AAACcj3vEAcBJunTpou3btys1NdXaWrdurX79+ln/rlSpklauXGm9Z8+ePTp06JDCw8MlSeHh4dq+fbuOHTtmxSQmJsrX11dNmjRxepsAAABQdFwRBwAnqVatmq6//nqHMh8fH9WoUcMqHzhwoIYPH66AgAD5+vrq8ccfV3h4uG688UZJUteuXdWkSRM98MADmjhxotLS0jR69GjFxsaWu6veAAAAKBiJOACUI1OmTJGbm5t69+6t7OxsRUVF6Y033rD2u7u7a8mSJXr00UcVHh4uHx8fxcTEaPz48S6sdcnkTVNnijoAAKgoSMQBwIVWr17t8Nrb21szZszQjBkzLvie0NBQLV26tIxrBgAAgLLCPeIAAAAAADgRiTgAAAAAAE5EIg4AAAAAgBORiAMAAAAA4EQk4gAAAAAAOBGJOAAAAAAATlTqifjYsWNls9kctkaNGln7T506pdjYWNWoUUNVq1ZV7969lZ6e7nCMQ4cOKTo6WlWqVFFgYKBGjBihM2fOlHZVAQAAAABwujJ5jnjTpk319ddf/+8kHv87zbBhw/Tll19q4cKF8vPzU1xcnO666y6tW7dOknT27FlFR0fLbrdr/fr1Onr0qPr3769KlSrppZdeKovqAgDKgYiEiHxlSTFJLqgJAABA2SqTqekeHh6y2+3WVrNmTUlSRkaG3n33XU2ePFm33HKLwsLCNHv2bK1fv14bNmyQJK1YsUK7du3Shx9+qJYtW6p79+6aMGGCZsyYoZycnLKoLgAAAFyIGZUAKpoyScT37t2r4OBgXXvtterXr58OHTokSUpJSdHp06cVGRlpxTZq1Ei1a9dWcnKyJCk5OVnNmjVTUFCQFRMVFaXMzEzt3LnzgufMzs5WZmamwwYAAIDLQ9OmTXX06FFr+/bbb619w4YN0xdffKGFCxdqzZo1OnLkiO666y5rf96MypycHK1fv14JCQmaM2eOxowZ44qmAMAllfrU9Hbt2mnOnDlq2LChjh49qnHjxqlDhw7asWOH0tLS5OnpKX9/f4f3BAUFKS0tTZKUlpbmkITn7c/bdyHx8fEaN25c6TYGAAAATpE3o/J8eTMq582bp1tuuUWSNHv2bDVu3FgbNmzQjTfeaM2o/PrrrxUUFKSWLVtqwoQJGjlypMaOHStPT88Cz5mdna3s7GzrNRdyADhLqV8R7969u+6++241b95cUVFRWrp0qU6cOKGPP/64tE/lYNSoUcrIyLC2w4cPl+n5AAAAUHpcMaMyPj5efn5+1hYSElJGrQMAR2X++DJ/f39dd9112rdvn+x2u3JycnTixAmHmPT0dOsvoHa7Pd89P3mvC/oraR4vLy/5+vo6bAAAACj/8mZULlu2TDNnztSBAwfUoUMH/fHHH2U6o5ILOQBcpcwT8aysLO3fv1+1atVSWFiYKlWqpJUrV1r79+zZo0OHDik8PFySFB4eru3bt+vYsWNWTGJionx9fdWkSZOyri4AAACczFUzKrmQA8BVSj0Rf+qpp7RmzRodPHhQ69ev15133il3d3fde++98vPz08CBAzV8+HAlJSUpJSVFAwYMUHh4uG688UZJUteuXdWkSRM98MAD2rp1q5YvX67Ro0crNjZWXl5epV1dAAAAlDPOmlEJAK5S6on4f/7zH917771q2LCh7rnnHtWoUUMbNmzQVVddJUmaMmWKbrvtNvXu3VsdO3aU3W7XokWLrPe7u7tryZIlcnd3V3h4uO6//371799f48ePL+2qAgAAoBxiRiWAK12pr5o+f/78i+739vbWjBkzNGPGjAvGhIaGaunSpaVdNQAAAJRDTz31lHr27KnQ0FAdOXJEzz//fIEzKgMCAuTr66vHH3/8gjMqJ06cqLS0NGZUVgARCRHWv5NiklxYE6DoSj0RBwAAAIoib0blb7/9pquuukrt27fPN6PSzc1NvXv3VnZ2tqKiovTGG29Y78+bUfnoo48qPDxcPj4+iomJYUYlgHKLRBwAAAAuxYxKABUNiTgAoNxi2iEAALgSlfnjy4ArQURChENCAAAAAADFRSIOALgs8AcxAABwpSARBwAAAADAiUjEAQAAAABwIhZrA4shAQAAAIATcUUcAAAAAAAn4oo4AOCyUtAsHmb2AACAywlXxAEAAAAAcCIScQAAAAAAnIhEHACcJD4+Xm3atFG1atUUGBioXr16ac+ePQ4xp06dUmxsrGrUqKGqVauqd+/eSk9Pd4g5dOiQoqOjVaVKFQUGBmrEiBE6c+aMM5tSbvBscQAAcDkiEQcAJ1mzZo1iY2O1YcMGJSYm6vTp0+ratatOnjxpxQwbNkxffPGFFi5cqDVr1ujIkSO66667rP1nz55VdHS0cnJytH79eiUkJGjOnDkaM2aMK5oEAACAYmCxNgBwkmXLljm8njNnjgIDA5WSkqKOHTsqIyND7777rubNm6dbbrlFkjR79mw1btxYGzZs0I033qgVK1Zo165d+vrrrxUUFKSWLVtqwoQJGjlypMaOHStPT09XNA0AAABFQCIOAC6SkZEhSQoICJAkpaSk6PTp04qMjLRiGjVqpNq1ays5OVk33nijkpOT1axZMwUFBVkxUVFRevTRR7Vz5061atUq33mys7OVnZ1tvc7MzCyrJpULBU1VZyV1AABQnjA1HQBcIDc3V0OHDtXNN9+s66+/XpKUlpYmT09P+fv7O8QGBQUpLS3Nijk3Cc/bn7evIPHx8fLz87O2kJCQUm4NAAAAioJEHABcIDY2Vjt27ND8+fPL/FyjRo1SRkaGtR0+fLjMzwkAAIALY2o6ADhZXFyclixZorVr1+qaa66xyu12u3JycnTixAmHq+Lp6emy2+1WzKZNmxyOl7eqel7M+by8vOTl5VXKrQAAAEBxcUUcAJzEGKO4uDh9+umnWrVqlerWreuwPywsTJUqVdLKlSutsj179ujQoUMKDw+XJIWHh2v79u06duyYFZOYmChfX181adLEOQ0BAABAiXBFHACcJDY2VvPmzdNnn32matWqWfd0+/n5qXLlyvLz89PAgQM1fPhwBQQEyNfXV48//rjCw8N14403SpK6du2qJk2a6IEHHtDEiROVlpam0aNHKzY2lqveAAAAlwkScQBwkpkzZ0qSOnfu7FA+e/ZsPfjgg5KkKVOmyM3NTb1791Z2draioqL0xhtvWLHu7u5asmSJHn30UYWHh8vHx0cxMTEaP368s5oBAACAEiIRryDyHufDI3wA1zHGXDLG29tbM2bM0IwZMy4YExoaqqVLl5Zm1SqMcx9txngIAABchUQcAHDFK+jZ4vyBEgAAuAqLtQEAAAAA4EQk4gAAAAAAOBFT0wEAFRr3jQMAAGfjijgAAP9fREJEgfeTAwAAlCYScQAAAAAAnIhEHAAAAAAAJ+IecQAAzlPQ9PRz7x/n0WcAAKAkSMQBACgE7h0HAAClhUQcAIBiYsV1AABQHNwjDgAAAACAE5GIAwBQCnj0GQAAKCwScQAAAAAAnIh7xCswrtwAQOnjvnEAAHApJOLn4MtT+ccjgwAAAABc7piaDgAAAACAE5GIAwBQRljADQAAFISp6RdQGlOgL8ep7hf7wljabbgc+wcAAAAASopEvAwUlMxyReR/6AsAAAAAFRmJ+CVcKmksjSu5V/KVYZJuAAAAAHBEIl5CrOJduq7kP0oAAAAAgMRibVcUFgUCAAAAgPKPK+IudLF7ycvqanBJEvVLXa0uSd35AwIAAACKi1mquNyQiJeS0k4kCztFuzwmsKVVJwZUAAAAAFciEvHLQEEJaWGSXWck6VfKOQAAAADAWUjEr0AkxwAAAABQfpGIX0ZIfgt2sX652L3sF9oPAACAyxPf83C5IBHHZamwf5QoywXmAAAAAKA4SMRRYbh6RkF5T/rLe/2AyxlXaADA+Qr67pc3Bhf2e2FJngbEeI+LKdeJ+IwZM/Tqq68qLS1NLVq00PTp09W2bVtXVwtXoOJcYc/DIFu+XakJEOMjAFwYYyQupKgXZvjuh7JSbhPxBQsWaPjw4Zo1a5batWunqVOnKioqSnv27FFgYKCrqwdYSmNAB4qC8REALowxEmWtoD/yX+r73cWuknMFvWKyGWOMqytRkHbt2qlNmzb617/+JUnKzc1VSEiIHn/8cT3zzDOXfH9mZqb8/PyUkZEhX1/fQp2TBAnlgasG4Sv1l0Bxr4gXZwxxFleMjxJjpDMU9JjKgr7kXWmfU1x+rtQxsqTjI5/NouN3S+FdahHii8XBeQo7jpTLK+I5OTlKSUnRqFGjrDI3NzdFRkYqOTm5wPdkZ2crOzvbep2RkSHp744orDN/nSlmjYHSk/czGz0vOt++L+/70vp33v5zywpS2Li8n/+ifGZKQ2HrV1znfq6L0ra82PL2t0pXjY8SY6QzdJjVodhleQr6LBU0nhT2vQUdp6CxqCjHKaqijhPn1qmsxpaK7koZI0t7fHT279ArAb9bCu9iY39x4kpb3nhb2N85FzvGuccpre+6zlLo8dGUQ7/88ouRZNavX+9QPmLECNO2bdsC3/P8888bSWxsbGyluh0+fNgZw16hMT6ysbGVp+1yHyMZH9nY2Mpqu9T4WC6viBfHqFGjNHz4cOt1bm6ujh8/rho1ashms13y/ZmZmQoJCdHhw4fL3RSr8ox+Kzr6rOhc0WfGGP3xxx8KDg52yvnKUknHR4mf29JEX5Yu+rN0FbY/r5QxkvGxdNEX/0NfOKpI/VHY8bFcJuI1a9aUu7u70tPTHcrT09Nlt9sLfI+Xl5e8vLwcyvz9/Yt8bl9f3yv+h6Ms0G9FR58VnbP7zM/Pz2nnKixXjo8SP7elib4sXfRn6SpMf14JYyTjY9mgL/6HvnBUUfqjMOOjmxPqUWSenp4KCwvTypUrrbLc3FytXLlS4eHhLqwZALgW4yMAXBhjJIDLRbm8Ii5Jw4cPV0xMjFq3bq22bdtq6tSpOnnypAYMGODqqgGASzE+AsCFMUYCuByU20T8//7v//Tf//5XY8aMUVpamlq2bKlly5YpKCioTM7n5eWl559/Pt/0JFwc/VZ09FnR0WeOnD0+SvwflCb6snTRn6XrSuhPvkO6Dn3xP/SFI/ojv3L7HHEAAAAAAK5E5fIecQAAAAAArlQk4gAAAAAAOBGJOAAAAAAATkQiDgAAAACAE1WoRHzGjBmqU6eOvL291a5dO23atOmi8QsXLlSjRo3k7e2tZs2aaenSpU6qaflSlH7buXOnevfurTp16shms2nq1KnOq2g5UpQ+e/vtt9WhQwdVr15d1atXV2Rk5CV/Nq9ERemzRYsWqXXr1vL395ePj49atmypDz74wIm1rViKOnbib2vXrlXPnj0VHBwsm82mxYsXO+w3xmjMmDGqVauWKleurMjISO3du9c1lS3n4uPj1aZNG1WrVk2BgYHq1auX9uzZ4xBz6tQpxcbGqkaNGqpatap69+6t9PR0F9W4fJs5c6aaN28uX19f+fr6Kjw8XF999ZW1n74svIowPpbW5+/QoUOKjo5WlSpVFBgYqBEjRujMmTPObEqpe/nll2Wz2TR06FCrrCL1xS+//KL7779fNWrUUOXKldWsWTNt2bLF2l+Y33PHjx9Xv3795OvrK39/fw0cOFBZWVnOboprmApi/vz5xtPT07z33ntm586dZtCgQcbf39+kp6cXGL9u3Trj7u5uJk6caHbt2mVGjx5tKlWqZLZv3+7kmrtWUftt06ZN5qmnnjIfffSRsdvtZsqUKc6tcDlQ1D677777zIwZM8z3339vdu/ebR588EHj5+dn/vOf/zi55q5T1D5LSkoyixYtMrt27TL79u0zU6dONe7u7mbZsmVOrvmVr6j/N/ifpUuXmn/+859m0aJFRpL59NNPHfa//PLLxs/PzyxevNhs3brV3H777aZu3brmr7/+ck2Fy7GoqCgze/Zss2PHDpOammp69OhhateubbKysqyYRx55xISEhJiVK1eaLVu2mBtvvNHcdNNNLqx1+fX555+bL7/80vz4449mz5495tlnnzWVKlUyO3bsMMbQl4VVUcbH0vj8nTlzxlx//fUmMjLSfP/992bp0qWmZs2aZtSoUa5oUqnYtGmTqVOnjmnevLkZMmSIVV5R+uL48eMmNDTUPPjgg2bjxo3mp59+MsuXLzf79u2zYgrze65bt26mRYsWZsOGDeabb74x9evXN/fee68rmuR0FSYRb9u2rYmNjbVenz171gQHB5v4+PgC4++55x4THR3tUNauXTvz8MMPl2k9y5ui9tu5QkNDK2QiXpI+M+bvAbpatWomISGhrKpY7pS0z4wxplWrVmb06NFlUb0KrTT+b2DyJeK5ubnGbrebV1991So7ceKE8fLyMh999JELanh5OXbsmJFk1qxZY4z5u+8qVapkFi5caMXs3r3bSDLJycmuquZlpXr16uadd96hL4ugoo6Pxfn8LV261Li5uZm0tDQrZubMmcbX19dkZ2c7twGl4I8//jANGjQwiYmJplOnTlYiXpH6YuTIkaZ9+/YX3F+Y33O7du0ykszmzZutmK+++srYbDbzyy+/lF3ly4kKMTU9JydHKSkpioyMtMrc3NwUGRmp5OTkAt+TnJzsEC9JUVFRF4y/EhWn3yq60uizP//8U6dPn1ZAQEBZVbNcKWmfGWO0cuVK7dmzRx07dizLqlY4jAFl58CBA0pLS3PoWz8/P7Vr146+LYSMjAxJssbJlJQUnT592qE/GzVqpNq1a9Ofl3D27FnNnz9fJ0+eVHh4OH1ZSBV5fCzO5y85OVnNmjVTUFCQFRMVFaXMzEzt3LnTibUvHbGxsYqOjs6XK1Skvvj888/VunVr3X333QoMDFSrVq309ttvW/sL83suOTlZ/v7+at26tRUTGRkpNzc3bdy40XmNcZEKkYj/+uuvOnv2rMMPvCQFBQUpLS2twPekpaUVKf5KVJx+q+hKo89Gjhyp4ODgfIP7laq4fZaRkaGqVavK09NT0dHRmj59um699dayrm6FwhhQdvL6j74tutzcXA0dOlQ333yzrr/+ekl/96enp6f8/f0dYunPC9u+fbuqVq0qLy8vPfLII/r000/VpEkT+rKQKur4WNzP34W+V+ftu5zMnz9f3333neLj4/Ptq0h98dNPP2nmzJlq0KCBli9frkcffVRPPPGEEhISJBXu91xaWpoCAwMd9nt4eCggIOCy6ovi8nB1BQD8z8svv6z58+dr9erV8vb2dnV1yrVq1aopNTVVWVlZWrlypYYPH65rr71WnTt3dnXVAJSh2NhY7dixQ99++62rq3JZa9iwoVJTU5WRkaFPPvlEMTExWrNmjaurhXKuon/+Dh8+rCFDhigxMbHCf0/Lzc1V69at9dJLL0mSWrVqpR07dmjWrFmKiYlxce0uDxXiinjNmjXl7u6eb8XC9PR02e32At9jt9uLFH8lKk6/VXQl6bNJkybp5Zdf1ooVK9S8efOyrGa5Utw+c3NzU/369dWyZUs9+eST6tOnT4F/nUbxMQaUnbz+o2+LJi4uTkuWLFFSUpKuueYaq9xutysnJ0cnTpxwiKc/L8zT01P169dXWFiY4uPj1aJFC73++uv0ZSFVxPGxJJ+/C32vztt3uUhJSdGxY8d0ww03yMPDQx4eHlqzZo2mTZsmDw8PBQUFVZi+qFWrlpo0aeJQ1rhxYx06dEhS4X7P2e12HTt2zGH/mTNndPz48cuqL4qrQiTinp6eCgsL08qVK62y3NxcrVy5UuHh4QW+Jzw83CFekhITEy8YfyUqTr9VdMXts4kTJ2rChAlatmyZw30yFUFp/Zzl5uYqOzu7LKpYYTEGlJ26devKbrc79G1mZqY2btxI3xbAGKO4uDh9+umnWrVqlerWreuwPywsTJUqVXLozz179ujQoUP0ZyHljaH0ZeFUpPGxND5/4eHh2r59u0PSlZiYKF9f33zJXHnWpUsXbd++XampqdbWunVr9evXz/p3RemLm2++Od9j7H788UeFhoZKKtzvufDwcJ04cUIpKSlWzKpVq5Sbm6t27do5oRUu5uLF4pxm/vz5xsvLy8yZM8fs2rXLDB482Pj7+1srFj7wwAPmmWeeseLXrVtnPDw8zKRJk8zu3bvN888/X2EfX1aUfsvOzjbff/+9+f77702tWrXMU089Zb7//nuzd+9eVzXB6YraZy+//LLx9PQ0n3zyiTl69Ki1/fHHH65qgtMVtc9eeukls2LFCrN//36za9cuM2nSJOPh4WHefvttVzXhinWp/xtc2B9//GGNh5LM5MmTzffff29+/vlnY8zfn31/f3/z2WefmW3btpk77riDx5ddwKOPPmr8/PzM6tWrHcbJP//804p55JFHTO3atc2qVavMli1bTHh4uAkPD3dhrcuvZ555xqxZs8YcOHDAbNu2zTzzzDPGZrOZFStWGGPoy8KqKONjaXz+8h7Z1bVrV5OammqWLVtmrrrqqsvukV0FOXfVdGMqTl9s2rTJeHh4mBdffNHs3bvXzJ0711SpUsV8+OGHVkxhfs9169bNtGrVymzcuNF8++23pkGDBjy+7Eo0ffp0U7t2bePp6Wnatm1rNmzYYO3r1KmTiYmJcYj/+OOPzXXXXWc8PT1N06ZNzZdffunkGpcPRem3AwcOGEn5tk6dOjm/4i5UlD4LDQ0tsM+ef/5551fchYrSZ//85z9N/fr1jbe3t6levboJDw838+fPd0GtK4aL/d/gwpKSkgr8bOf9LOfm5prnnnvOBAUFGS8vL9OlSxezZ88e11a6nCqoHyWZ2bNnWzF//fWXeeyxx0z16tVNlSpVzJ133mmOHj3qukqXYw899JAJDQ01np6e5qqrrjJdunSxknBj6MuiqAjjY2l9/g4ePGi6d+9uKleubGrWrGmefPJJc/r0aSe3pvSdn4hXpL744osvzPXXX2+8vLxMo0aNzFtvveWwvzC/53777Tdz7733mqpVqxpfX18zYMCACnMxymaMMc659g4AAAAAACrEPeIAAAAAAJQXJOIAAAAAADgRiTgAAAAAAE5EIg4AAAAAgBORiAMAAAAA4EQk4gAAAAAAOBGJOAAAAAAATkQiDgAAAACAE5GIA+fo3Lmzhg4d6upqAECZe/DBB9WrV68ivWfs2LFq2bJlmdQHAMoLxkc4A4k4iu3BBx+UzWbLt3Xr1s3VVSu2RYsWacKECSU6Rl6/PPLII/n2xcbGymaz6cEHHyzUsVJSUmSz2bRhw4YC93fp0kV33XVXSaoLoAwV58tcaTt48KBsNptSU1NLfKynnnpKK1euLHmlAFR4jI+o6EjEUSLdunXT0aNHHbaPPvqoTM+Zk5NTZscOCAhQtWrVSnyckJAQzZ8/X3/99ZdVdurUKc2bN0+1a9cu9HHCwsLUokULvffee/n2HTx4UElJSRo4cGCJ6wsAhVG1alXVqFHD1dUAgHKH8RFFRSKOEvHy8pLdbnfYqlevbu232Wx65513dOedd6pKlSpq0KCBPv/8c4dj7NixQ927d1fVqlUVFBSkBx54QL/++qu1v3PnzoqLi9PQoUNVs2ZNRUVFSZI+//xzNWjQQN7e3oqIiFBCQoJsNptOnDihkydPytfXV5988onDuRYvXiwfHx/98ccfBbbn/KnpderU0UsvvaSHHnpI1apVU+3atfXWW29dsl9uuOEGhYSEaNGiRVbZokWLVLt2bbVq1cohNjc3V/Hx8apbt64qV66sFi1aONR74MCBWrBggf7880+H982ZM0e1atW6rGcgABVdYca/J554Qk8//bQCAgJkt9s1duxYh2P88MMPat++vby9vdWkSRN9/fXXstlsWrx4sSSpbt26kqRWrVrJZrOpc+fODu+fNGmSatWqpRo1aig2NlanT5++YH3Pn3qZd0WrOMd47733VLt2bVWtWlWPPfaYzp49q4kTJ8putyswMFAvvviiw/tsNpvefPNN3XbbbapSpYoaN26s5ORk7du3T507d5aPj49uuukm7d+//yI9DuBywfjI+HilIxFHmRs3bpzuuecebdu2TT169FC/fv10/PhxSdKJEyd0yy23qFWrVtqyZYuWLVum9PR03XPPPQ7HSEhIkKenp9atW6dZs2bpwIED6tOnj3r16qWtW7fq4Ycf1j//+U8r3sfHR3379tXs2bMdjjN79mz16dOnSFe9X3vtNbVu3Vrff/+9HnvsMT366KPas2fPJd/30EMPOZz/vffe04ABA/LFxcfH6/3339esWbO0c+dODRs2TPfff7/WrFkjSerXr5+ys7MdknNjjBISEvTggw/K3d290G0BUH4UZfzz8fHRxo0bNXHiRI0fP16JiYmSpLNnz6pXr16qUqWKNm7cqLfeesthLJSkTZs2SZK+/vprHT161OEPhElJSdq/f7+SkpKUkJCgOXPmaM6cOUVqR3GOsX//fn311VdatmyZPvroI7377ruKjo7Wf/7zH61Zs0avvPKKRo8erY0bNzq8b8KECerfv79SU1PVqFEj3XfffXr44Yc1atQobdmyRcYYxcXFFan+AMofxkfGxwrBAMUUExNj3N3djY+Pj8P24osvWjGSzOjRo63XWVlZRpL56quvjDHGTJgwwXTt2tXhuIcPHzaSzJ49e4wxxnTq1Mm0atXKIWbkyJHm+uuvdyj75z//aSSZ33//3RhjzMaNG427u7s5cuSIMcaY9PR04+HhYVavXn3BNnXq1MkMGTLEeh0aGmruv/9+63Vubq4JDAw0M2fOvGi/3HHHHebYsWPGy8vLHDx40Bw8eNB4e3ub//73v+aOO+4wMTExxhhjTp06ZapUqWLWr1/vcIyBAweae++913rdt29f06lTJ+v1ypUrjSSzd+/eC9YDgOvljQcFKez41759e4eYNm3amJEjRxpjjPnqq6+Mh4eHOXr0qLU/MTHRSDKffvqpMcaYAwcOGEnm+++/z1e30NBQc+bMGavs7rvvNv/3f/93wfY8//zzpkWLFiU+RpUqVUxmZqZVFhUVZerUqWPOnj1rlTVs2NDEx8dbr8//fZKcnGwkmXfffdcq++ijj4y3t/cFzw2g/GB8LPgYjI8Vh4cLcn9cQSIiIjRz5kyHsoCAAIfXzZs3t/7t4+MjX19fHTt2TJK0detWJSUlqWrVqvmOvX//fl133XWS/r5X+lx79uxRmzZtHMratm2b73XTpk2VkJCgZ555Rh9++KFCQ0PVsWPHIrXx3PrbbDbZ7Xar/hdz1VVXKTo6WnPmzJExRtHR0apZs6ZDzL59+/Tnn3/q1ltvdSjPyclxmML+0EMPKSoqSvv371e9evX03nvvqVOnTqpfv36R2gKg/Cjs+HfuGCRJtWrVssagPXv2KCQkRHa73dp//lh4MU2bNnWYVVOrVi1t3769SO0ozjHq1KnjMDMpKChI7u7ucnNzcyg7f6w9ty+CgoIkSc2aNXMoO3XqlDIzM+Xr61ukdgAoPxgfGR8rAhJxlIiPj88lk8FKlSo5vLbZbMrNzZUkZWVlqWfPnnrllVfyva9WrVoO5ymOf/zjH5oxY4aeeeYZzZ49WwMGDJDNZivSMS5W/0t56KGHrGlAM2bMyLc/KytLkvTll1/q6quvdtjn5eVl/btLly6qXbu25syZoxEjRmjRokV68803i9QOAOVLYce/koxBl1Iaxy7OMQp6T2GOc25M3lheUFlp9Q8A12B8dHwP4+OViUQcLnXDDTfo3//+t+rUqSMPj8L/ODZs2FBLly51KNu8eXO+uPvvv19PP/20pk2bpl27dikmJqbEdS6Kbt26KScnRzabzVpk7lxNmjSRl5eXDh06pE6dOl3wOG5ubhowYIDeffddXX311fL09FSfPn3KsuoAylhxx79zNWzYUIcPH1Z6erp1BeT8sdDT01PS3/dLAsDlgPERFQGLtaFEsrOzlZaW5rCdu6LlpcTGxur48eO69957tXnzZu3fv1/Lly/XgAEDLjooPvzww/rhhx80cuRI/fjjj/r444+txS/OveJdvXp13XXXXRoxYoS6du2qa665pthtLQ53d3ft3r1bu3btKnBRtWrVqumpp57SsGHDlJCQoP379+u7777T9OnTlZCQ4BA7YMAA/fLLL3r22Wd17733qnLlys5qBoASyMjIUGpqqsN2+PDhYo9/57r11ltVr149xcTEaNu2bVq3bp1Gjx4t6X9jYWBgoCpXrmwtdpSRkVFmbQWAomB8REVGIo4SWbZsmWrVquWwtW/fvtDvDw4O1rp163T27Fl17dpVzZo109ChQ+Xv7+9wL8z56tatq08++USLFi1S8+bNNXPmTGslzHOndEt/P/4rJydHDz30UPEaWUK+vr4XvRdnwoQJeu655xQfH6/GjRurW7du+vLLL61HauSpXbu2IiMj9fvvv7usLQCKbvXq1WrVqpXDNm7cuGKPf+dyd3fX4sWLlZWVpTZt2ugf//iHNRZ6e3tLkjw8PDRt2jS9+eabCg4O1h133FFmbQWAomB8REVmM8YYV1cCKA0vvviiZs2apcOHDzuUf/DBBxo2bJiOHDliTUECgCvVunXr1L59e+3bt0/16tVzdXUAoNxgfER5wj3iuGy98cYbatOmjWrUqKF169bp1VdfdXg+4p9//qmjR4/q5Zdf1sMPP0wSDuCK9Omnn6pq1apq0KCB9u3bpyFDhujmm2/mSyaACo/xEeUZiTguW3v37tULL7yg48ePq3bt2nryySc1atQoa//EiRP14osvqmPHjg7lAHAl+eOPPzRy5EgdOnRINWvWVGRkpF577TVXVwsAXI7xEeUZU9MBAAAAAHAiFmsDAAAAAMCJSMQBAAAAAHAiEnEAAAAAAJyIRBwAAAAAACciEQcAAAAAwIlIxAEAAAAAcCIScQAAAAAAnIhEHAAAAAAAJ/p/4+viNjS5H1cAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1200x400 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig2, ax2 = plt.subplots(1, 3, figsize=(12, 4))\n",
    "\n",
    "# Plot 1\n",
    "edep = hits.arrays()['TotalEnergyDeposit']\n",
    "ax2[0].hist(edep, 100, facecolor='green', alpha=0.75)\n",
    "ax2[0].set_title('Total Energy Deposit')\n",
    "ax2[0].set_xlabel('Energy in MeV')\n",
    "\n",
    "# Plot 1\n",
    "stepLength = hits.arrays()['StepLength']\n",
    "n, bins, patches = ax2[1].hist(stepLength, 100, facecolor='green', alpha=0.75)\n",
    "ax2[1].set_title('Step Length')\n",
    "ax2[1].set_xlabel('Length in mm')\n",
    "\n",
    "# Plot 2\n",
    "trackLength = hits.arrays()['TrackLength']\n",
    "n, bins, patches = ax2[2].hist(trackLength, 100, facecolor='green', alpha=0.75)\n",
    "ax2[2].set_title('Track Length')\n",
    "ax2[2].set_xlabel('Length in mm')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:10.203427Z",
     "start_time": "2024-01-23T19:56:10.194894Z"
    },
    "cell_id": "00009-e878ae6a-d878-4821-8fde-7698918ba06a",
    "execution_millis": 7,
    "execution_start": 1605009148550,
    "output_cleared": false,
    "source_hash": "10778042"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of Singles : 9\n",
      "Singles: \n",
      "PostPosition_X  PostPosition_Y  PostPosition_Z  TotalEnergyDeposit  PreStepUniqueVolumeID  GlobalTime  LocalTime  StepLength  TrackLength  "
     ]
    }
   ],
   "source": [
    "singles = f['Singles_spect_crystal']\n",
    "print('Number of Singles :', len(singles))\n",
    "print('Singles: ')\n",
    "for i in singles.keys():\n",
    "    print(i, ' ', end='')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "cell_id": "00010-6b9ee439-23f5-4451-be31-751890e44398",
    "output_cleared": false
   },
   "source": [
    "# Plot 3: energy  of the singles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:10.280827Z",
     "start_time": "2024-01-23T19:56:10.201275Z"
    },
    "cell_id": "00011-e7c5fc9c-0e27-4bb6-8821-438cf7653ff5",
    "execution_millis": 577,
    "execution_start": 1605009148571,
    "output_cleared": false,
    "source_hash": "401d2dd6"
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'Energy in keV')"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAHHCAYAAAChjmJTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABA3UlEQVR4nO3de1xVdb7/8fcGZXvBDSJyG5VIy0upGZXuM2WSDKhkmjplmmJijobNeBkzTv3UciYb7WIX00fHSay0NCebSVNDBW0ULSny0uioB8NSwFRAUVHh+/tjHqzjDrxgKCx8PR+P9Yi9vp+11vf7dQPv1l5r4TDGGAEAANiIV3V3AAAAoLIIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMACUlpYmh8OhtLS06u4KqoDD4dDUqVOruxvAVUWAAaqJw+G4rOVyQsULL7ygTz75pEr7N2zYsAv2qV69elV6rOo0depUj7E1aNBALVq0UO/evTV//nwVFxdXdxd/sU2bNmnq1KnKz8+v7q4AVaZOdXcAuF699957Hq/fffddpaSklFvftm3bS+7rhRde0IABA9S3b9+q7KKcTqfmzZtXbr23t3eVHqcmmDNnjnx9fVVcXKwff/xRq1ev1vDhwzVr1iwtX75czZs3r+4uXrZTp06pTp3/+/G+adMmPffccxo2bJj8/f2rr2NAFSLAANXk0Ucf9Xi9efNmpaSklFtfnerUqVNj+lNUVKSGDRtetf0PGDBAgYGB1uvJkydr4cKFGjp0qH77299q8+bNV+3YVa02nSEDLoSPkIAarKioSBMmTFDz5s3ldDrVunVrvfTSSzr/j8g7HA4VFRVpwYIF1scgw4YNkyR9//33euKJJ9S6dWvVr19fTZo00W9/+1vt37+/yvqYnJwsh8OhjRs3avz48WratKkaNmyoBx98UIcPHy5Xv3LlSt1zzz1q2LChGjVqpLi4OO3cudOjZtiwYfL19dW+ffvUq1cvNWrUSIMHD5b0n7MLv//97xUYGKhGjRrpgQce0I8//uhx3UdqaqocDoeWLVtW7viLFi2Sw+FQenr6Jcc2ePBgjRgxQlu2bFFKSopH25YtW9SjRw/5+fmpQYMGuvfee7Vx40aPmrKPp3bt2qWHHnpILpdLTZo00R/+8AedPn3ao/bcuXOaNm2aWrZsKafTqRtuuEH//d//Xe4jrK1btyo2NlaBgYGqX7++IiIiNHz4cI+a8+di6tSpmjhxoiQpIiLCeo9U5XsAqA6cgQFqKGOMHnjgAaWmpiohIUG33XabVq9erYkTJ+rHH3/Uq6++Kuk/H0WNGDFCd911l0aOHClJatmypSTpq6++0qZNmzRw4EA1a9ZM+/fv15w5c9StWzd99913atCgwSX78dNPP5Vb5+PjI5fL5bHuySefVOPGjTVlyhTt379fs2bN0pgxY7R48WKr5r333lN8fLxiY2P1l7/8RSdPntScOXN0991365tvvtENN9xg1Z47d06xsbG6++679dJLL1l9HTZsmJYsWaIhQ4aoS5cuWr9+veLi4jz60q1bNzVv3lwLFy7Ugw8+6NG2cOFCtWzZUm63+5Jjl6QhQ4bo7bff1ueff67f/OY3kqR169apZ8+eioyM1JQpU+Tl5aX58+frvvvu0xdffKG77rrLYx8PPfSQbrjhBk2fPl2bN2/W66+/rmPHjundd9+1akaMGKEFCxZowIABmjBhgrZs2aLp06frX//6lxXE8vLyFBMTo6ZNm+rpp5+Wv7+/9u/fr48//viC/e/Xr5/+/e9/64MPPtCrr75qnWVq2rTpZY0fqLEMgBohMTHRnP8t+cknnxhJ5k9/+pNH3YABA4zD4TB79+611jVs2NDEx8eX2+fJkyfLrUtPTzeSzLvvvmutS01NNZJMamqqtS4+Pt5IqnCJjY216ubPn28kmejoaFNaWmqtHzdunPH29jb5+fnGGGOOHz9u/P39zeOPP+7Rn5ycHOPn5+exvuzYTz/9tEdtRkaGkWTGjh3rsX7YsGFGkpkyZYq1LikpyTidTuv4xhiTl5dn6tSp41E3ZcoUI8kcPny43FwZY8yxY8eMJPPggw8aY4wpLS01N910k4mNjfUY78mTJ01ERIT5zW9+U27fDzzwgMc+n3jiCSPJfPvtt8YYYzIzM40kM2LECI+6P/7xj0aSWbdunTHGmGXLlhlJ5quvvqqwr2V+PhczZ840kkxWVtZFtwPshI+QgBrqs88+k7e3t37/+997rJ8wYYKMMVq5cuUl91G/fn3r67Nnz+rIkSNq1aqV/P399fXXX19y+3r16iklJaXc8uKLL5arHTlypBwOh/X6nnvuUUlJib7//ntJUkpKivLz8/XII4/op59+shZvb2917txZqamp5fY5evRoj9erVq2SJD3xxBMe65988sly2w4dOlTFxcVaunSptW7x4sU6d+5cpa7r8fX1lSQdP35ckpSZmak9e/Zo0KBBOnLkiDWOoqIide/eXRs2bFBpaanHPhITEyvs72effebx3/Hjx3vUTZgwQZK0YsUKSbIuwF2+fLnOnj172WMAaiM+QgJqqO+//15hYWFq1KiRx/qyu5LKgsHFnDp1StOnT9f8+fP1448/elw7U1BQcMntvb29FR0dfVn9bdGihcfrxo0bS5KOHTsmSdqzZ48k6b777qtw+59/JFWnTh01a9bMY933338vLy8vRUREeKxv1apVuf21adNGd955pxYuXKiEhARJ//n4qEuXLhXWX8iJEyckyfp3KBtHfHz8BbcpKCiwxi9JN910k0d7y5Yt5eXlZV2HUjaun/crJCRE/v7+1r/1vffeq/79++u5557Tq6++qm7duqlv374aNGiQnE7nZY8JqA0IMEAt9uSTT2r+/PkaO3as3G63/Pz85HA4NHDgwHJnCX6pC91aXRaayo733nvvKSQkpFzd+bf9Sv+5hdvL65edJB46dKj+8Ic/6IcfflBxcbE2b96sN998s1L72LFjh6T/C0ll45g5c6Zuu+22CrcpO2tzIeefqbqc9ee3L126VJs3b9ann35q3er98ssva/PmzZc8LlCbEGCAGio8PFxr1qzR8ePHPc7C7Nq1y2ovc6FffEuXLlV8fLxefvlla93p06er5YFmZRcWBwUFXfZZnZ8LDw9XaWmpsrKyPM5q7N27t8L6gQMHavz48frggw906tQp1a1bVw8//HCljln2XJ7Y2FhJ/zcOl8t12ePYs2ePx1mjvXv3qrS01LpouWxce/bs8XjuT25urvLz8z3+rSWpS5cu6tKli/785z9r0aJFGjx4sD788EONGDGiwuNfKhgBdsQ1MEAN1atXL5WUlJQ7Y/Dqq6/K4XCoZ8+e1rqGDRtWGEq8vb09PjaSpDfeeEMlJSVXpc8XExsbK5fLpRdeeKHC6zcquuW6on1I0ltvveWx/o033qiwPjAwUD179tT777+vhQsXqkePHh7PermURYsWad68eXK73erevbskKTIyUi1bttRLL71kfbx0qXHMnj27wv6W/Rv26tVLkjRr1iyPuldeeUWSrLusjh07Vu7fs+ws0MWeGFz2/ByexIvahDMwQA3Vu3dvRUVF6ZlnntH+/fvVsWNHff755/r73/+usWPHWmcCpP/8Ul2zZo1eeeUVhYWFKSIiQp07d9b999+v9957T35+fmrXrp3S09O1Zs0aNWnS5LL6cO7cOb3//vsVtj344IOVerCcy+XSnDlzNGTIEN1+++0aOHCgmjZtquzsbK1YsUK//vWvL/nxTmRkpPr3769Zs2bpyJEj1m3U//73vyVVfKZh6NChGjBggCRp2rRpF9z30qVL5evrqzNnzlhP4t24caM6duyojz76yKrz8vLSvHnz1LNnT91yyy167LHH9Ktf/Uo//vijUlNT5XK59Omnn3rsOysrSw888IB69Oih9PR0vf/++xo0aJA6duwoSerYsaPi4+P19ttvKz8/X/fee6++/PJLLViwQH379lVUVJQkacGCBXrrrbf04IMPqmXLljp+/Lj+53/+Ry6XywpBF5o3SXrmmWc0cOBA1a1bV717976qDwYErrpqvQcKgOXnt1Eb859bj8eNG2fCwsJM3bp1zU033WRmzpzpcfuuMcbs2rXLdO3a1dSvX99Ism6pPnbsmHnsscdMYGCg8fX1NbGxsWbXrl0mPDzc47bryt5GrfNuyS27jfrnt/ZWtM+y9bGxscbPz8/Uq1fPtGzZ0gwbNsxs3brV49gNGzascJ6KiopMYmKiCQgIML6+vqZv375m9+7dRpJ58cUXy9UXFxebxo0bGz8/P3Pq1Kly7WW3Opct9erVM82aNTP333+/eeedd8zp06cr7Mc333xj+vXrZ5o0aWKcTqcJDw83Dz30kFm7dm25fX/33XdmwIABplGjRqZx48ZmzJgx5fpy9uxZ89xzz5mIiAhTt25d07x5c5OUlORx/K+//to88sgjpkWLFsbpdJqgoCBz//33e8ydMeVvozbGmGnTpplf/epXxsvLi1uqUSs4jPnZ+UgAsJnMzEx16tRJ77//vvXE3jLnzp1TWFiYevfurb/+9a/XtF9Tp07Vc889p8OHD1fqoysAl8Y1MABs5dSpU+XWzZo1S15eXuratWu5tk8++USHDx/W0KFDr0X3AFwjXAMDwFZmzJihjIwMRUVFqU6dOlq5cqVWrlypkSNHevzF6C1btmjbtm2aNm2aOnXqpHvvvbcaew2gqhFgANjKf/3XfyklJUXTpk3TiRMn1KJFC02dOlXPPPOMR92cOXP0/vvv67bbblNycnL1dBbAVcM1MAAAwHa4BgYAANgOAQYAANhOrb0GprS0VAcPHlSjRo14jDYAADZhjNHx48cVFhZ20b+HVmsDzMGDBz3uSAAAAPZx4MCBcn+R/ny1NsCU/fG7AwcOyOVyVXNvAADA5SgsLFTz5s09/ohtRWptgCn72MjlchFgAACwmUtd/sFFvAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMAAAwHYIMABgU1ELohS1IKq6uwFUCwIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwnUoFmDlz5qhDhw5yuVxyuVxyu91auXKl1d6tWzc5HA6PZdSoUR77yM7OVlxcnBo0aKCgoCBNnDhR586d86hJS0vT7bffLqfTqVatWik5OfnKRwgAAGqdOpUpbtasmV588UXddNNNMsZowYIF6tOnj7755hvdcsstkqTHH39czz//vLVNgwYNrK9LSkoUFxenkJAQbdq0SYcOHdLQoUNVt25dvfDCC5KkrKwsxcXFadSoUVq4cKHWrl2rESNGKDQ0VLGxsVUxZgAAYHMOY4z5JTsICAjQzJkzlZCQoG7duum2227TrFmzKqxduXKl7r//fh08eFDBwcGSpLlz52rSpEk6fPiwfHx8NGnSJK1YsUI7duywths4cKDy8/O1atWqy+5XYWGh/Pz8VFBQIJfL9UuGCAA1UtSCKElSanxqNfcEqDqX+/v7iq+BKSkp0YcffqiioiK53W5r/cKFCxUYGKhbb71VSUlJOnnypNWWnp6u9u3bW+FFkmJjY1VYWKidO3daNdHR0R7Hio2NVXp6+kX7U1xcrMLCQo8FAADUTpX6CEmStm/fLrfbrdOnT8vX11fLli1Tu3btJEmDBg1SeHi4wsLCtG3bNk2aNEm7d+/Wxx9/LEnKycnxCC+SrNc5OTkXrSksLNSpU6dUv379Cvs1ffp0Pffcc5UdDgAAsKFKB5jWrVsrMzNTBQUFWrp0qeLj47V+/Xq1a9dOI0eOtOrat2+v0NBQde/eXfv27VPLli2rtOM/l5SUpPHjx1uvCwsL1bx586t6TAAAUD0q/RGSj4+PWrVqpcjISE2fPl0dO3bUa6+9VmFt586dJUl79+6VJIWEhCg3N9ejpux1SEjIRWtcLtcFz75IktPptO6OKlsAAEDt9IufA1NaWqri4uIK2zIzMyVJoaGhkiS3263t27crLy/PqklJSZHL5bI+hnK73Vq7dq3HflJSUjyuswEAANe3Sn2ElJSUpJ49e6pFixY6fvy4Fi1apLS0NK1evVr79u3TokWL1KtXLzVp0kTbtm3TuHHj1LVrV3Xo0EGSFBMTo3bt2mnIkCGaMWOGcnJy9OyzzyoxMVFOp1OSNGrUKL355pt66qmnNHz4cK1bt05LlizRihUrqn70AADAlioVYPLy8jR06FAdOnRIfn5+6tChg1avXq3f/OY3OnDggNasWaNZs2apqKhIzZs3V//+/fXss89a23t7e2v58uUaPXq03G63GjZsqPj4eI/nxkRERGjFihUaN26cXnvtNTVr1kzz5s3jGTAAAMDyi58DU1PxHBgAtR3PgUFtdNWfAwMAAFBdCDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2CDAAAMB2KhVg5syZow4dOsjlcsnlcsntdmvlypVW++nTp5WYmKgmTZrI19dX/fv3V25ursc+srOzFRcXpwYNGigoKEgTJ07UuXPnPGrS0tJ0++23y+l0qlWrVkpOTr7yEQIAgFqnUgGmWbNmevHFF5WRkaGtW7fqvvvuU58+fbRz505J0rhx4/Tpp5/qo48+0vr163Xw4EH169fP2r6kpERxcXE6c+aMNm3apAULFig5OVmTJ0+2arKyshQXF6eoqChlZmZq7NixGjFihFavXl1FQwYAAHbnMMaYX7KDgIAAzZw5UwMGDFDTpk21aNEiDRgwQJK0a9cutW3bVunp6erSpYtWrlyp+++/XwcPHlRwcLAkae7cuZo0aZIOHz4sHx8fTZo0SStWrNCOHTusYwwcOFD5+flatWrVZfersLBQfn5+KigokMvl+iVDBIAaKWpBlCQpNT61mnsCVJ3L/f19xdfAlJSU6MMPP1RRUZHcbrcyMjJ09uxZRUdHWzVt2rRRixYtlJ6eLklKT09X+/btrfAiSbGxsSosLLTO4qSnp3vso6ymbB8XUlxcrMLCQo8FAADUTpUOMNu3b5evr6+cTqdGjRqlZcuWqV27dsrJyZGPj4/8/f096oODg5WTkyNJysnJ8QgvZe1lbRerKSws1KlTpy7Yr+nTp8vPz89amjdvXtmhAQAAm6h0gGndurUyMzO1ZcsWjR49WvHx8fruu++uRt8qJSkpSQUFBdZy4MCB6u4SAAC4SupUdgMfHx+1atVKkhQZGamvvvpKr732mh5++GGdOXNG+fn5HmdhcnNzFRISIkkKCQnRl19+6bG/sruUzq/5+Z1Lubm5crlcql+//gX75XQ65XQ6KzscAABgQ7/4OTClpaUqLi5WZGSk6tatq7Vr11ptu3fvVnZ2ttxutyTJ7XZr+/btysvLs2pSUlLkcrnUrl07q+b8fZTVlO0DAACgUmdgkpKS1LNnT7Vo0ULHjx/XokWLlJaWptWrV8vPz08JCQkaP368AgIC5HK59OSTT8rtdqtLly6SpJiYGLVr105DhgzRjBkzlJOTo2effVaJiYnW2ZNRo0bpzTff1FNPPaXhw4dr3bp1WrJkiVasWFH1owcAALZUqQCTl5enoUOH6tChQ/Lz81OHDh20evVq/eY3v5Ekvfrqq/Ly8lL//v1VXFys2NhYvfXWW9b23t7eWr58uUaPHi23262GDRsqPj5ezz//vFUTERGhFStWaNy4cXrttdfUrFkzzZs3T7GxsVU0ZAAAYHe/+DkwNRXPgQFQ2/EcGNRGV/05MAAAANWFAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGyHAAMAAGynUgFm+vTpuvPOO9WoUSMFBQWpb9++2r17t0dNt27d5HA4PJZRo0Z51GRnZysuLk4NGjRQUFCQJk6cqHPnznnUpKWl6fbbb5fT6VSrVq2UnJx8ZSMEAAC1TqUCzPr165WYmKjNmzcrJSVFZ8+eVUxMjIqKijzqHn/8cR06dMhaZsyYYbWVlJQoLi5OZ86c0aZNm7RgwQIlJydr8uTJVk1WVpbi4uIUFRWlzMxMjR07ViNGjNDq1at/4XABAEBtUKcyxatWrfJ4nZycrKCgIGVkZKhr167W+gYNGigkJKTCfXz++ef67rvvtGbNGgUHB+u2227TtGnTNGnSJE2dOlU+Pj6aO3euIiIi9PLLL0uS2rZtq3/+85969dVXFRsbW9kxAgCAWuYXXQNTUFAgSQoICPBYv3DhQgUGBurWW29VUlKSTp48abWlp6erffv2Cg4OttbFxsaqsLBQO3futGqio6M99hkbG6v09PRf0l0AAFBLVOoMzPlKS0s1duxY/frXv9att95qrR80aJDCw8MVFhambdu2adKkSdq9e7c+/vhjSVJOTo5HeJFkvc7JybloTWFhoU6dOqX69euX609xcbGKi4ut14WFhVc6NAAAUMNdcYBJTEzUjh079M9//tNj/ciRI62v27dvr9DQUHXv3l379u1Ty5Ytr7ynlzB9+nQ999xzV23/AACg5riij5DGjBmj5cuXKzU1Vc2aNbtobefOnSVJe/fulSSFhIQoNzfXo6bsddl1MxeqcblcFZ59kaSkpCQVFBRYy4EDByo/MAAAYAuVCjDGGI0ZM0bLli3TunXrFBERccltMjMzJUmhoaGSJLfbre3btysvL8+qSUlJkcvlUrt27ayatWvXeuwnJSVFbrf7gsdxOp1yuVweCwAAqJ0qFWASExP1/vvva9GiRWrUqJFycnKUk5OjU6dOSZL27dunadOmKSMjQ/v379c//vEPDR06VF27dlWHDh0kSTExMWrXrp2GDBmib7/9VqtXr9azzz6rxMREOZ1OSdKoUaP0v//7v3rqqae0a9cuvfXWW1qyZInGjRtXxcMHAAB2VKkAM2fOHBUUFKhbt24KDQ21lsWLF0uSfHx8tGbNGsXExKhNmzaaMGGC+vfvr08//dTah7e3t5YvXy5vb2+53W49+uijGjp0qJ5//nmrJiIiQitWrFBKSoo6duyol19+WfPmzeMWagAAIElyGGNMdXfiaigsLJSfn58KCgr4OAlArRS1IEqSlBqfWs09AarO5f7+5m8hAQAA2yHAAAAA2yHAAAAA2yHAAAAA2yHAADYRtSDKumgTAK53BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7BBgAAGA7lQow06dP15133qlGjRopKChIffv21e7duz1qTp8+rcTERDVp0kS+vr7q37+/cnNzPWqys7MVFxenBg0aKCgoSBMnTtS5c+c8atLS0nT77bfL6XSqVatWSk5OvrIRAgCAWqdSAWb9+vVKTEzU5s2blZKSorNnzyomJkZFRUVWzbhx4/Tpp5/qo48+0vr163Xw4EH169fPai8pKVFcXJzOnDmjTZs2acGCBUpOTtbkyZOtmqysLMXFxSkqKkqZmZkaO3asRowYodWrV1fBkAEAgN05jDHmSjc+fPiwgoKCtH79enXt2lUFBQVq2rSpFi1apAEDBkiSdu3apbZt2yo9PV1dunTRypUrdf/99+vgwYMKDg6WJM2dO1eTJk3S4cOH5ePjo0mTJmnFihXasWOHdayBAwcqPz9fq1atuqy+FRYWys/PTwUFBXK5XFc6RKDGiFoQJUlKjU+t5p6gpuA9gdrocn9//6JrYAoKCiRJAQEBkqSMjAydPXtW0dHRVk2bNm3UokULpaenS5LS09PVvn17K7xIUmxsrAoLC7Vz506r5vx9lNWU7aMixcXFKiws9FgAAEDtdMUBprS0VGPHjtWvf/1r3XrrrZKknJwc+fj4yN/f36M2ODhYOTk5Vs354aWsvaztYjWFhYU6depUhf2ZPn26/Pz8rKV58+ZXOjQAAFDDXXGASUxM1I4dO/Thhx9WZX+uWFJSkgoKCqzlwIED1d0lAABwldS5ko3GjBmj5cuXa8OGDWrWrJm1PiQkRGfOnFF+fr7HWZjc3FyFhIRYNV9++aXH/sruUjq/5ud3LuXm5srlcql+/foV9snpdMrpdF7JcAAAgM1U6gyMMUZjxozRsmXLtG7dOkVERHi0R0ZGqm7dulq7dq21bvfu3crOzpbb7ZYkud1ubd++XXl5eVZNSkqKXC6X2rVrZ9Wcv4+ymrJ9AACA61ulzsAkJiZq0aJF+vvf/65GjRpZ16z4+fmpfv368vPzU0JCgsaPH6+AgAC5XC49+eSTcrvd6tKliyQpJiZG7dq105AhQzRjxgzl5OTo2WefVWJionUGZdSoUXrzzTf11FNPafjw4Vq3bp2WLFmiFStWVPHwAQCAHVXqDMycOXNUUFCgbt26KTQ01FoWL15s1bz66qu6//771b9/f3Xt2lUhISH6+OOPrXZvb28tX75c3t7ecrvdevTRRzV06FA9//zzVk1ERIRWrFihlJQUdezYUS+//LLmzZun2NjYKhgyAACwu1/0HJiajOfAoLbhmR/4Od4TqI2uyXNgAAAAqgMBBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2E6lA8yGDRvUu3dvhYWFyeFw6JNPPvFoHzZsmBwOh8fSo0cPj5qjR49q8ODBcrlc8vf3V0JCgk6cOOFRs23bNt1zzz2qV6+emjdvrhkzZlR+dAAAoFaqdIApKipSx44dNXv27AvW9OjRQ4cOHbKWDz74wKN98ODB2rlzp1JSUrR8+XJt2LBBI0eOtNoLCwsVExOj8PBwZWRkaObMmZo6darefvvtynYXAADUQnUqu0HPnj3Vs2fPi9Y4nU6FhIRU2Pavf/1Lq1at0ldffaU77rhDkvTGG2+oV69eeumllxQWFqaFCxfqzJkzeuedd+Tj46NbbrlFmZmZeuWVVzyCDgAAuD5dlWtg0tLSFBQUpNatW2v06NE6cuSI1Zaeni5/f38rvEhSdHS0vLy8tGXLFquma9eu8vHxsWpiY2O1e/duHTt2rMJjFhcXq7Cw0GMBAAC1U5UHmB49eujdd9/V2rVr9Ze//EXr169Xz549VVJSIknKyclRUFCQxzZ16tRRQECAcnJyrJrg4GCPmrLXZTU/N336dPn5+VlL8+bNq3poAACghqj0R0iXMnDgQOvr9u3bq0OHDmrZsqXS0tLUvXv3qj6cJSkpSePHj7deFxYWEmIAAKilrvpt1DfeeKMCAwO1d+9eSVJISIjy8vI8as6dO6ejR49a182EhIQoNzfXo6bs9YWurXE6nXK5XB4LAACona56gPnhhx905MgRhYaGSpLcbrfy8/OVkZFh1axbt06lpaXq3LmzVbNhwwadPXvWqklJSVHr1q3VuHHjq91lAABQw1U6wJw4cUKZmZnKzMyUJGVlZSkzM1PZ2dk6ceKEJk6cqM2bN2v//v1au3at+vTpo1atWik2NlaS1LZtW/Xo0UOPP/64vvzyS23cuFFjxozRwIEDFRYWJkkaNGiQfHx8lJCQoJ07d2rx4sV67bXXPD4iAgAA169KB5itW7eqU6dO6tSpkyRp/Pjx6tSpkyZPnixvb29t27ZNDzzwgG6++WYlJCQoMjJSX3zxhZxOp7WPhQsXqk2bNurevbt69eqlu+++2+MZL35+fvr888+VlZWlyMhITZgwQZMnT+YWagAAIOkKLuLt1q2bjDEXbF+9evUl9xEQEKBFixZdtKZDhw764osvKts9AABwHeBvIQEAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANshwAAAANupdIDZsGGDevfurbCwMDkcDn3yySce7cYYTZ48WaGhoapfv76io6O1Z88ej5qjR49q8ODBcrlc8vf3V0JCgk6cOOFRs23bNt1zzz2qV6+emjdvrhkzZlR+dAAAoFaqdIApKipSx44dNXv27ArbZ8yYoddff11z587Vli1b1LBhQ8XGxur06dNWzeDBg7Vz506lpKRo+fLl2rBhg0aOHGm1FxYWKiYmRuHh4crIyNDMmTM1depUvf3221cwRAAAUNvUqewGPXv2VM+ePStsM8Zo1qxZevbZZ9WnTx9J0rvvvqvg4GB98sknGjhwoP71r39p1apV+uqrr3THHXdIkt544w316tVLL730ksLCwrRw4UKdOXNG77zzjnx8fHTLLbcoMzNTr7zyikfQAQAA16cqvQYmKytLOTk5io6Ottb5+fmpc+fOSk9PlySlp6fL39/fCi+SFB0dLS8vL23ZssWq6dq1q3x8fKya2NhY7d69W8eOHavw2MXFxSosLPRYAABA7VSlASYnJ0eSFBwc7LE+ODjYasvJyVFQUJBHe506dRQQEOBRU9E+zj/Gz02fPl1+fn7W0rx5818+IAAAUCPVmruQkpKSVFBQYC0HDhyo7i4BAICrpEoDTEhIiCQpNzfXY31ubq7VFhISory8PI/2c+fO6ejRox41Fe3j/GP8nNPplMvl8lgAAEDtVKUBJiIiQiEhIVq7dq21rrCwUFu2bJHb7ZYkud1u5efnKyMjw6pZt26dSktL1blzZ6tmw4YNOnv2rFWTkpKi1q1bq3HjxlXZZQAAYEOVDjAnTpxQZmamMjMzJf3nwt3MzExlZ2fL4XBo7Nix+tOf/qR//OMf2r59u4YOHaqwsDD17dtXktS2bVv16NFDjz/+uL788ktt3LhRY8aM0cCBAxUWFiZJGjRokHx8fJSQkKCdO3dq8eLFeu211zR+/PgqGzgAALCvSt9GvXXrVkVFRVmvy0JFfHy8kpOT9dRTT6moqEgjR45Ufn6+7r77bq1atUr16tWztlm4cKHGjBmj7t27y8vLS/3799frr79utfv5+enzzz9XYmKiIiMjFRgYqMmTJ3MLNQAAkCQ5jDGmujtxNRQWFsrPz08FBQVcD4NaIWrBf/7HITU+tZp7gpqC9wRqo8v9/V1r7kICAADXDwIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwHQIMAACwnSoPMFOnTpXD4fBY2rRpY7WfPn1aiYmJatKkiXx9fdW/f3/l5uZ67CM7O1txcXFq0KCBgoKCNHHiRJ07d66quwoAAGyqztXY6S233KI1a9b830Hq/N9hxo0bpxUrVuijjz6Sn5+fxowZo379+mnjxo2SpJKSEsXFxSkkJESbNm3SoUOHNHToUNWtW1cvvPDC1eguAACwmasSYOrUqaOQkJBy6wsKCvTXv/5VixYt0n333SdJmj9/vtq2bavNmzerS5cu+vzzz/Xdd99pzZo1Cg4O1m233aZp06Zp0qRJmjp1qnx8fK5GlwEAgI1clWtg9uzZo7CwMN14440aPHiwsrOzJUkZGRk6e/asoqOjrdo2bdqoRYsWSk9PlySlp6erffv2Cg4OtmpiY2NVWFionTt3XvCYxcXFKiws9FgAAEDtVOUBpnPnzkpOTtaqVas0Z84cZWVl6Z577tHx48eVk5MjHx8f+fv7e2wTHBysnJwcSVJOTo5HeClrL2u7kOnTp8vPz89amjdvXrUDAwAANUaVf4TUs2dP6+sOHTqoc+fOCg8P15IlS1S/fv2qPpwlKSlJ48ePt14XFhYSYgAAqKWu+m3U/v7+uvnmm7V3716FhITozJkzys/P96jJzc21rpkJCQkpd1dS2euKrqsp43Q65XK5PBYAAFA7XfUAc+LECe3bt0+hoaGKjIxU3bp1tXbtWqt99+7dys7OltvtliS53W5t375deXl5Vk1KSopcLpfatWt3tbsLAABsoMo/QvrjH/+o3r17Kzw8XAcPHtSUKVPk7e2tRx55RH5+fkpISND48eMVEBAgl8ulJ598Um63W126dJEkxcTEqF27dhoyZIhmzJihnJwcPfvss0pMTJTT6azq7gIAABuq8gDzww8/6JFHHtGRI0fUtGlT3X333dq8ebOaNm0qSXr11Vfl5eWl/v37q7i4WLGxsXrrrbes7b29vbV8+XKNHj1abrdbDRs2VHx8vJ5//vmq7ioAALCpKg8wH3744UXb69Wrp9mzZ2v27NkXrAkPD9dnn31W1V0DAAC1BH8LCQAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgAA2A4BBgBQ40UtiFLUgqjq7gZqEAIMbI8fbABw/anRAWb27Nm64YYbVK9ePXXu3FlffvlldXcJqJEIcQCulZry86bGBpjFixdr/PjxmjJlir7++mt17NhRsbGxysvLq+6uATVWTfnBAgBXW53q7sCFvPLKK3r88cf12GOPSZLmzp2rFStW6J133tHTTz9dzb2D3ZX9kk+NT62ytqvVj8tdX9ma813OeH6+z6qaAwC4EjUywJw5c0YZGRlKSkqy1nl5eSk6Olrp6enV2LPyLvQL6Pwf9hdqO3/9hX7h8Euiav18ni/2i/5K265ETTlrcqUBraqD3bVEMAPsqUYGmJ9++kklJSUKDg72WB8cHKxdu3ZVuE1xcbGKi4ut1wUFBZKkwsLCKu9f3KK4cuvumXvPBesv1HaxbSpTcy2sGLTC+vrn4y9rO399Ret+yXEvZz81Za7s6Odzd6Xvzdrwb1DdY/j5e/5i33tlrsbPuZrm3Klzkq6PsdZ0V/vfomy/xpiLF5oa6McffzSSzKZNmzzWT5w40dx1110VbjNlyhQjiYWFhYWFhaUWLAcOHLhoVqiRZ2ACAwPl7e2t3Nxcj/W5ubkKCQmpcJukpCSNHz/eel1aWqqjR4+qSZMmcjgcVda3wsJCNW/eXAcOHJDL5aqy/dod83JhzE3FmJeKMS8VY14qVhvnxRij48ePKyws7KJ1NTLA+Pj4KDIyUmvXrlXfvn0l/SeQrF27VmPGjKlwG6fTKafT6bHO39//qvXR5XLVmjdLVWJeLoy5qRjzUjHmpWLMS8Vq27z4+fldsqZGBhhJGj9+vOLj43XHHXforrvu0qxZs1RUVGTdlQQAAK5fNTbAPPzwwzp8+LAmT56snJwc3XbbbVq1alW5C3sBAMD1p8YGGEkaM2bMBT8yqi5Op1NTpkwp93HV9Y55uTDmpmLMS8WYl4oxLxW7nufFYcyl7lMCAACoWWrsnxIAAAC4EAIMAACwHQIMAACwHQIMAACwHQJMJc2ePVs33HCD6tWrp86dO+vLL7+s7i5dU1OnTpXD4fBY2rRpY7WfPn1aiYmJatKkiXx9fdW/f/9yT1SuDTZs2KDevXsrLCxMDodDn3zyiUe7MUaTJ09WaGio6tevr+joaO3Zs8ej5ujRoxo8eLBcLpf8/f2VkJCgEydOXMNRVL1LzcuwYcPKvX969OjhUVMb52X69Om688471ahRIwUFBalv377avXu3R83lfO9kZ2crLi5ODRo0UFBQkCZOnKhz585dy6FUqcuZl27dupV7z4waNcqjprbNy5w5c9ShQwfr4XRut1srV6602q/H90pFCDCVsHjxYo0fP15TpkzR119/rY4dOyo2NlZ5eXnV3bVr6pZbbtGhQ4es5Z///KfVNm7cOH366af66KOPtH79eh08eFD9+vWrxt5eHUVFRerYsaNmz55dYfuMGTP0+uuva+7cudqyZYsaNmyo2NhYnT592qoZPHiwdu7cqZSUFC1fvlwbNmzQyJEjr9UQropLzYsk9ejRw+P988EHH3i018Z5Wb9+vRITE7V582alpKTo7NmziomJUVFRkVVzqe+dkpISxcXF6cyZM9q0aZMWLFig5ORkTZ48uTqGVCUuZ14k6fHHH/d4z8yYMcNqq43z0qxZM7344ovKyMjQ1q1bdd9996lPnz7auXOnpOvzvVKhKvnri9eJu+66yyQmJlqvS0pKTFhYmJk+fXo19uramjJliunYsWOFbfn5+aZu3brmo48+stb961//MpJMenr6NerhtSfJLFu2zHpdWlpqQkJCzMyZM611+fn5xul0mg8++MAYY8x3331nJJmvvvrKqlm5cqVxOBzmxx9/vGZ9v5p+Pi/GGBMfH2/69OlzwW2uh3kxxpi8vDwjyaxfv94Yc3nfO5999pnx8vIyOTk5Vs2cOXOMy+UyxcXF13YAV8nP58UYY+69917zhz/84YLbXA/zYowxjRs3NvPmzeO9ch7OwFymM2fOKCMjQ9HR0dY6Ly8vRUdHKz09vRp7du3t2bNHYWFhuvHGGzV48GBlZ2dLkjIyMnT27FmPOWrTpo1atGhxXc1RVlaWcnJyPObBz89PnTt3tuYhPT1d/v7+uuOOO6ya6OhoeXl5acuWLde8z9dSWlqagoKC1Lp1a40ePVpHjhyx2q6XeSkoKJAkBQQESLq875309HS1b9/e42nksbGxKiwstP7P3O5+Pi9lFi5cqMDAQN16661KSkrSyZMnrbbaPi8lJSX68MMPVVRUJLfbzXvlPDX6Sbw1yU8//aSSkpJyf8ogODhYu3btqqZeXXudO3dWcnKyWrdurUOHDum5557TPffcox07dignJ0c+Pj7l/ohmcHCwcnJyqqfD1aBsrBW9V8racnJyFBQU5NFep04dBQQE1Oq56tGjh/r166eIiAjt27dP//3f/62ePXsqPT1d3t7e18W8lJaWauzYsfr1r3+tW2+9VZIu63snJyenwvdUWZvdVTQvkjRo0CCFh4crLCxM27Zt06RJk7R79259/PHHkmrvvGzfvl1ut1unT5+Wr6+vli1bpnbt2ikzM/O6f6+UIcCgUnr27Gl93aFDB3Xu3Fnh4eFasmSJ6tevX409gx0MHDjQ+rp9+/bq0KGDWrZsqbS0NHXv3r0ae3btJCYmaseOHR7XjuHC83L+9U/t27dXaGiounfvrn379qlly5bXupvXTOvWrZWZmamCggItXbpU8fHxWr9+fXV3q0bhI6TLFBgYKG9v73JXeufm5iokJKSaelX9/P39dfPNN2vv3r0KCQnRmTNnlJ+f71Fzvc1R2Vgv9l4JCQkpd/H3uXPndPTo0etqrm688UYFBgZq7969kmr/vIwZM0bLly9XamqqmjVrZq2/nO+dkJCQCt9TZW12dqF5qUjnzp0lyeM9UxvnxcfHR61atVJkZKSmT5+ujh076rXXXrvu3yvnI8BcJh8fH0VGRmrt2rXWutLSUq1du1Zut7sae1a9Tpw4oX379ik0NFSRkZGqW7euxxzt3r1b2dnZ19UcRUREKCQkxGMeCgsLtWXLFmse3G638vPzlZGRYdWsW7dOpaWl1g/o68EPP/ygI0eOKDQ0VFLtnRdjjMaMGaNly5Zp3bp1ioiI8Gi/nO8dt9ut7du3ewS8lJQUuVwutWvX7toMpIpdal4qkpmZKUke75naNi8VKS0tVXFx8XX7XqlQdV9FbCcffvihcTqdJjk52Xz33Xdm5MiRxt/f3+NK79puwoQJJi0tzWRlZZmNGzea6OhoExgYaPLy8owxxowaNcq0aNHCrFu3zmzdutW43W7jdruruddV7/jx4+abb74x33zzjZFkXnnlFfPNN9+Y77//3hhjzIsvvmj8/f3N3//+d7Nt2zbTp08fExERYU6dOmXto0ePHqZTp05my5Yt5p///Ke56aabzCOPPFJdQ6oSF5uX48ePmz/+8Y8mPT3dZGVlmTVr1pjbb7/d3HTTTeb06dPWPmrjvIwePdr4+fmZtLQ0c+jQIWs5efKkVXOp751z586ZW2+91cTExJjMzEyzatUq07RpU5OUlFQdQ6oSl5qXvXv3mueff95s3brVZGVlmb///e/mxhtvNF27drX2URvn5emnnzbr1683WVlZZtu2bebpp582DofDfP7558aY6/O9UhECTCW98cYbpkWLFsbHx8fcddddZvPmzdXdpWvq4YcfNqGhocbHx8f86le/Mg8//LDZu3ev1X7q1CnzxBNPmMaNG5sGDRqYBx980Bw6dKgae3x1pKamGknllvj4eGPMf26l/n//7/+Z4OBg43Q6Tffu3c3u3bs99nHkyBHzyCOPGF9fX+Nyucxjjz1mjh8/Xg2jqToXm5eTJ0+amJgY07RpU1O3bl0THh5uHn/88XL/A1Ab56WiOZFk5s+fb9VczvfO/v37Tc+ePU39+vVNYGCgmTBhgjl79uw1Hk3VudS8ZGdnm65du5qAgADjdDpNq1atzMSJE01BQYHHfmrbvAwfPtyEh4cbHx8f07RpU9O9e3crvBhzfb5XKuIwxphrd74HAADgl+MaGAAAYDsEGAAAYDsEGAAAYDsEGAAAYDsEGAAAYDsEGAAAYDsEGAAAYDsEGAC1Srdu3TR27Nhq3weAq4sAA1zHhg0bJofDUW7p0aNHdXftin388ceaNm1atR3/ySefVNu2bStsy87Olre3t/7xj39c414BtQ8BBrjO9ejRQ4cOHfJYPvjgg6t6zDNnzly1fQcEBKhRo0ZXbf+XkpCQoF27dmnTpk3l2pKTkxUUFKRevXpVQ8+A2oUAA1znnE6nQkJCPJbGjRtb7Q6HQ/PmzdODDz6oBg0a6Kabbip3BmHHjh3q2bOnfH19FRwcrCFDhuinn36y2rt166YxY8Zo7NixCgwMVGxsrCTpH//4h2666SbVq1dPUVFRWrBggRwOh/Lz81VUVCSXy6WlS5d6HOuTTz5Rw4YNdfz48QrH8/OPf2644Qa98MILGj58uBo1aqQWLVro7bffrtQcrVixQn5+flq4cKEk6cCBA3rooYfk7++vgIAA9enTR/v375ck3Xbbbbr99tv1zjvveOzDGKPk5GTFx8erTp06lTo+gPIIMAAu6bnnntNDDz2kbdu2qVevXho8eLCOHj0qScrPz9d9992nTp06aevWrVq1apVyc3P10EMPeexjwYIF8vHx0caNGzV37lxlZWVpwIAB6tu3r7799lv97ne/0zPPPGPVN2zYUAMHDtT8+fM99jN//nwNGDCgUmdZXn75Zd1xxx365ptv9MQTT2j06NHavXv3ZW27aNEiPfLII1q4cKEGDx6ss2fPKjY2Vo0aNdIXX3yhjRs3ytfXVz169LDOLCUkJGjJkiUqKiqy9pOWlqasrCwNHz78svsN4CKq+Y9JAqhG8fHxxtvb2zRs2NBj+fOf/2zVSDLPPvus9frEiRNGklm5cqUxxphp06aZmJgYj/0eOHDASLL+Ave9995rOnXq5FEzadIkc+utt3qse+aZZ4wkc+zYMWOMMVu2bDHe3t7m4MGDxhhjcnNzTZ06dUxaWtoFx3TvvfeaP/zhD9br8PBw8+ijj1qvS0tLTVBQkJkzZ84l9/Hmm28aPz8/j+O99957pnXr1qa0tNRaV1xcbOrXr29Wr15tjDHm2LFjpl69eh5/bXrIkCHm7rvvvuAxAVQO5zGB61xUVJTmzJnjsS4gIMDjdYcOHayvGzZsKJfLpby8PEnSt99+q9TUVPn6+pbb9759+3TzzTdLkiIjIz3adu/erTvvvNNj3V133VXu9S233KIFCxbo6aef1vvvv6/w8HB17dq1UmM8v/8Oh0MhISFW/y9k6dKlysvL08aNGz36+e2332rv3r3lzgCdPn1a+/btkyT5+/urX79+eueddzRs2DAVFhbqb3/7m2bPnl2pfgO4MAIMcJ1r2LChWrVqddGaunXrerx2OBwqLS2VJJ04cUK9e/fWX/7yl3LbhYaGehznSowYMUKzZ8/W008/rfnz5+uxxx6Tw+Go1D4u1v8L6dSpk77++mu98847uuOOO6xjnjhxQpGRkdb1MOdr2rSp9XVCQoK6d++uvXv3KjU1Vd7e3vrtb39bqX4DuDACDIBf5Pbbb9ff/vY33XDDDZW6OLV169b67LPPPNZ99dVX5eoeffRRPfXUU3r99df13XffKT4+/hf3+XK0bNlSL7/8srp16yZvb2+9+eabkv4z3sWLFysoKEgul+uC20dFRSkiIkLz589XamqqBg4ceMUhDkB5XMQLXOeKi4uVk5PjsZx/B9GlJCYm6ujRo3rkkUf01Vdfad++fVq9erUee+wxlZSUXHC73/3ud9q1a5cmTZqkf//731qyZImSk5MlyeMMS+PGjdWvXz9NnDhRMTExatas2RWPtbJuvvlmpaam6m9/+5t1Z9PgwYMVGBioPn366IsvvlBWVpbS0tL0+9//Xj/88IO1rcPh0PDhwzVnzhylp6crISHhmvUbuB4QYIDr3KpVqxQaGuqx3H333Ze9fVhYmDZu3KiSkhLFxMSoffv2Gjt2rPz9/eXldeEfMREREVq6dKk+/vhjdejQQXPmzLHuQnI6nR61CQkJOnPmTLXcwdO6dWutW7dOH3zwgSZMmKAGDRpow4YNatGihfr166e2bdsqISFBp0+fLndGZtiwYSooKNAtt9yizp07X/O+A7WZwxhjqrsTACBJf/7znzV37lwdOHDAY/17772ncePG6eDBg/Lx8amm3gGoSbgGBkC1eeutt3TnnXeqSZMm2rhxo2bOnKkxY8ZY7SdPntShQ4f04osv6ne/+x3hBYCFj5AAVJs9e/aoT58+ateunaZNm6YJEyZo6tSpVvuMGTPUpk0bhYSEKCkpqfo6CqDG4SMkAABgO5yBAQAAtkOAAQAAtkOAAQAAtkOAAQAAtkOAAQAAtkOAAQAAtkOAAQAAtkOAAQAAtkOAAQAAtvP/AQnIeeydpYUbAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig3, ax3 = plt.subplots()\n",
    "energy = singles.arrays()['TotalEnergyDeposit']/0.001\n",
    "n, bins, patches = ax3.hist(energy, 200, facecolor='green', alpha=0.75)\n",
    "ax3.set_title('TotalEnergyDeposit')\n",
    "ax3.set_xlabel('Energy in keV')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-01-23T19:56:10.281059Z",
     "start_time": "2024-01-23T19:56:10.273679Z"
    }
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "deepnote_execution_queue": [],
  "deepnote_notebook_id": "2fdaf907-cc5d-4f6c-96dc-43e4d6245abf",
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
